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ABSTRACT 


This document is the Programmer's Guide for the GNAT Computer Program 
developed under MSC/TRW Task 705-2, "Apollo Cyrogenic Storage System (CSS) 
Analysis," Subtask 2. Detailed logic flowcharts and compiled program 
listings are provided for all program elements. 
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1 . 0 INTRODUCTION 


The General Numerical Analysis of Transport Computet Program (GNAT) was 
developed to describe the thermal stratification and fluid mixing occurring 
in the supercritical oxygen storage system of the Apollo Command and Service 
Module. These phenomena are governed in two dimensions by a system of four 

non-linear partial differential equations which specify the conservation 
of mass, momantum, and energy. These equations are formulated in rectilinear 
Eulerian coordinates and are solved simultaneously using an explicit finite 
difference technique which Is described In Reference 1. 
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2.0 PROGRAM DESCRIPTION 


The two-dimensional volume is divided into a grid of cubic nodes which 
are individually identified by subscripts (I, J) which run in the x-and 
y-directions respectively. The indicies I and J must have values between 
1 and 20 inclusive. The tank geometric configuration is defined by 
specifying the beginning (NG) and ending (NS) values of the J-subscrlpt for 

each value of I such that: 

NG(I) < J < NS (I) ; 1< I < 20 

The same limits apply to the I-subscript at a given value of J so that: 

NG (J) < I < NS(J); 1 < J < 20 

The program variables were defined in a self-consistent set of 
engineering units so that no conversion constants were required to formulate 
the governing equations. All program variables are expressed in units which 
are combinations of one or more of the following units: lbf, slug, ft, sec, 

°F. Input data, expressed in the commonly used units shown in Tables 3.2 and 3.3, are 
converted immediately upon entry into the program. Output data is converted 

to the commonly used units shown in Table 3.4 for ease in interpretation. All 
computations performed in the program use the consistent set of units shown in 
Table 2.1. 

The state of the fluid at any given time is specified by the density, 

,iai energy, and the x- and y- components of momentum at each node point I, J. 

I these values, the corresponding pressure, temperature, and velocities 

are computed. 

The program is started by specifying the initial state at rime t Q . This 
initial state may be specified by card input in which case zero velocity 
conditions are assumed, or it may { input from a previously generated data 
tape which contains the Information for a developed fluid state obatined from 


an earlier computer run. 



The program time is advanced in a stepwiBe manner from the Initial time 

by successively integrating the four governing equations at each node (I,J) 
within the fluid volume. This procedure is illustrated in the logic 
flow chart of the main program (Figure 2.1). 
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TABLE 2.1 


PROGRAM UNITS 


Variable 

Description 

Units 

P 

Pressure 

(lb f /ft 2 ) 

H 

Temperature, T 

(°R) 

R 

Density, p 

(slugs/ft 3 ) 

U, V 

Velocity 

(f t/sec) 

E 

Specific Internal Energy 

(£ t“lb f /slug) 

RU, RV 

Momentum 

(siug-f c/sec/ft 

RE 

Internal Energy 

(ft-lb f /ft 3 ) 

T, DT 

Time 

(sec) 

K 

Thermal Conductivity, k 

,ft~lbf . 

l ft-sec-°R ; 

vsc 

Absolute Viscosity, p 

(lb f - r >ec/f t 2 ) 



FIGURE 2.1 

FLOWCHART OF MAIN PROGRAM 
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3.0 STORAGE ALLOCATION 

The program requires the storage of the fluid state variables for 
each of the possible 400 node poirts. All node centered variables are 
dimensioned ?0 x 20. The values of density, energy, and x~ and y- 
momentum at the node faces are obtained by interpolation between nodes 
and are stored in arrays which are dimensioned 21 x 21 and are identified 
by the suffix x or y. The associated subscripts (1,J) refer to the left 
or lower face of node I,J. 

Moat of the progi*am variables are ultimately equlvalenced to the 
master Input/Output (1/0) variable A which is dimensioned 3700 and which 
is placed in the COMMON block DATA, The single equivalencing variable A 
was defined to simplify reading and writing of the state data and other 
necessary control data on magnetic tape to provide a program restart 
capability, 

The program data are combined in logical groupings depending upon 
their function. For example, the array CTL contains twenty locations which 
are used for program initiation and control. Table 3.1 show the relation- 
ship of these data blocks to the A-array. 

Tables 3.2 through 3.8 show the breakdown of the data blocks to the 
basic program variables. In addition to the data blocks equlvalenced 
to the A-array, other COMMON blocks are used to transfer data between the 
program elements. Table 3.9 shows a cross-reference between all data blocks 
in COMMON storage and the program elements in which the data blocks are used. 
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TABLE 3.1 

STORAGE ALLOCATION: A(3700) 


TYPE: Real 

PURPOSE: Provide a common I/O variable to which most program 
variables are equivalenced. 


A(l) - 

A(20) 

CTL 

(Table 3.2) 

A(21 - 

A(30) 

PROP 

(Table 3.3) 

A(31) - 

■ A (40) 

LABEL 

(Table 3.4) 

A(41) - 

■ A(90) 

LIMITS 

(Table 3.5) 

A(91) - 

■ A(100) 

AVE 

(Table 3,6) 

A(101) 

- A(2100) 

STATE 

(Table 3.7) 

A(2101) 

- A(3700) 

RATES 

(Table 3.8) 
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TABLE 3.2 

STORAGE ALLOCATION: CTL(20) 

TYPE: Real 

PURPOSE: Input of program control variables 


1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 


CTL(I) 


Descrip 1 1 o n (Input Units) 


TO Program start time (sec) 

DT Integration time step (sec) 

TSTOP Program stop time (sec) 

DTPR Data output time interval (sec) 

WDOT Tank mass flowrate (lbm/hr) 

DQHEAT Tank heater input (B/hr) 

2 

DQBC(l) Boundary heat flux-left wall (B/ft -hr) 

DQBC(2) Boundary heat flux-right wall (B/ft^-hr) 

2 

DQBC(3) Boundary heat flux-bottom wall (B/ft -hr) 

2 

DQBC(4) Boundary heat flux-top wall (B/ft -hr) 

GX Acceleration, x-component (g's) 

GY Acceleration, y-component (g’s) 

SCALE Scale factor 



I/O Tape control described in Table 3.2a 


IT APE (6) 

STOPFG Program stop flag. Stop if > 1. 
Equlvalenced to A(l) 
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TABLE 3.2a 

STORAGE ALLOCATION: ITAPE(6) 

TYPE: Integer 

PURPOSE: I/O Tape control variables. 


1 

ITAPE (I) 

Description 

1 

IUNITR 

Tape unit for input state data. 

2 

IFILER 

File of input data on IUNITR 

3 

IRECR 

Record of input data in IFILER 

A 

IUNITW 

Tape unit to output state data. 

5 

IFILEW 

File of output data on IUNITW 

6 

IRECW 

Record of output data In IFILEW 


ITAPE is equivalenced to CTL(IA) 


15 



TABLE 3.3 

STORAGE ALLOCATION: PR0P(1O) 


TYPE: 

PURPOSE: 

Real 
Input of 

property data 


I 


PR0P(I) 

Description (Input Units) 

1 

2 

3 


k 

Thermal conductivity (B/f t-hr-°R) 

4 


vsc 

Absolute viscosity (Poise) 

5 

6 


L 

Node dimension (ft) 

7 

8 


PO 

Initial pressure (psl) 

9 


HO 

Initial temperature (°R) 

10 


PPLAG 

Heater flag: l.=on, 0,=off. 



Equlvalenced to A(21). 
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TABLE 3.4 

STORAGE ALLOCATION: LABEL (10) 


TYPE: 

Real 


PURPOSE: 

Output array. 


I 


LABEL (I) 

1 


T 

2 


DTI 

3 


TSCALE 

4 


PBAR1 

5 


RBAR1 

6 


HMIN 

7 


HBAR1 

8 


HMAX 

9 


PCOL 

10 


WT 


Description (Output Units) 
Time (minutes) 

Program time step (sec) 

Scaled problem time (min) 

Average tank pressure (psi) 

3 

Average tank density (lbm/ ft ) 
Minimum node temperature (°R) 
Average temperature (°R) 

Maximum node temperature (°R) 
Collapse pressure (psi) 

Weight of fluid contained (lbm) 


Equivalenced to A (31) 
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TABLE 3.5 

STORAGE ALLOCATION: LIMITS (50) 


TYPE: Integer 

PURPOSE: I, J Indexing limits for geometry definition 


1 

1 

+ 

20 

21 

+ 

40 

41 

42 

43 


LIMITS (I) 

Description 

NG (1) 

Node starting index array 


J1=NG(I) or I1=NG (J) 

NG(20) 

NS (1) 

Node ending index array 

+ 

J2=NS(I) or I2=NS (J) 

NS (20) 

NO 

Minimum NG 

NE 

Maximum NS 

NODES 

Total number of nodes 


Equivalenced to A (41) 
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TABLE 3.6 

STORAGE ALLOCATION; AVE(IO) 

TYPE: Double precision 

PURPOSE: Double precision storage of average state variables. 


I 

AVE(l) 

Description 

1 

PBAR 

Average pressure (lbf/ft ) 

2 

+ 


3 

RBAR 

Average density (slugs /ft ) 

4 

i 


5 

REBAR 

3 

Average intefnal energy (ft-lbf/ft ) 


6 + 

7 

8 

9 

10 

Equlvalenced to A (91) 
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TABLE 3,7 

STORAGE ALLOCATION: STATE (2000) 

TYPE: Real 

PURPOSE: Contains the state variables P, R, U, V. 


I 

Rates (I) 

Description 

1 

P(l,l) 

2 

Relative pressure (lbf/ft ) 

+ 

i 


400 

P(20,20) 


401 

R(l,l) 

3 

Relative density (slugs/ft ) 

+ 



800 

R(20,20) 


801 

0 ( 1 , 1 ) 

x-component velocity (ft/ sec) 




1200 

U (20,20) 


1201 

V(l,lj 

y-component velocity (ft/sec) 

i 

+ 


1600 

V ( 20 , 20 ) 



Equivalence:! to A (101) 
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TABLE 3.8 

STORAGE ALLOCATION: RATES (1600) 


TYPE: Real 

PURPOSl • Contains the additional state variables RE, RU, RV 


l 

Rates (I) 

Description 

1 

— 


+ 



400 



401 

1 

RE (1,1) 

Relative internal energy (ft-lbf/ft 

T 

800 

RE (20, 20) 


801 

RU(1, 1) 

3 

x-component momentum (ft/sec/ft ) 

+ 

4 - 


1200 

RU (20,20) 


1201 

RV (1,1) 

3 

y-component momentum (ft/sec/ft ) 

4 

+ 


160C 

RV (20,20) 



Bquivalenced to A (2101) 
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DATA BLOCKS IN COMMON 


TABLE 3.9 STORAGE CROSS-REFERENCE 
PROGRAM ELEMENTS 
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RDTAPE 


4.0 SUBPROGRAM DESCRIPTION 


All program logic and the evaluation and integration of the governing 
equations are performed in the main program. However, extensive use is 
made of subroutine subprograms to perform support and peripheral functions. 

A summary of these subprograms and a brief description of the function served 
is given in Table 4.1. 

The subprograms are described in more detail in the following subsections. 

The operation of each subprogram is shown In the accompanying flowcharts, 
(Figures 4.1 through 4.16). 
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TABLE 4.1 
SUBPROGRAM SUMMARY 


Sub 

■Brass.™ 

Function 

1 . 

MEANA 

Face-centered interpolation for momentum 

2. 

MEANB 

Face-centered interpolation for mass and energy 

3. 

DIFF 

Differences for momentum 

4. 

DIFF2 

Differences for mass and energy 

5. 

HEATER 

Heater input simulation 

6 . 

BC0UT 

Withdrawal port boundary conditions 

7. 

TEMP 

Temperature calculation from internal energy 

8. 

BETA 

Internal energy calculation from temperature 

9. 

PRESS 

Ideal gas equation of state 

10. 

RTPRES 

Stewart's equation of state 

11. 

BULK 

Average tank quantities 

12. 

OUTPUT 

Data output control 

13. 

DISPLY 

Display of one state variable 

14. 

CR2TAP 

Input data display 

15. 

RDTAPE 

Data tape input 

16. 

TAPEI0 

Tape Input/output package 

17. 

OPTD 

Interpolation of Weber's oxygen density data 
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4.1 Subroutine MEANA 


CALLING SEQUENCE: 
INPUT DATA: 

OUTPUT DATA: 
USAGE: 

FUNCTION: 


CALL MEANA 

RU, RV, U, V, NO, NF, NG, NS 
RUX, RUY, RVX, RVY, UX, VY 

This subroutine is called from MAIN once each program 
time step. 

Linear interpolation between node centers is used 
to define x- and y~ components of velocity at node 
faces. Quadratic interpolation between three node 
centers is used to obtain the values of x- and y- 
momentum at each of the node faces. The signs of 
UX and VY are used to bias the interpolation in 
the upstream direction. A linear Interpolation is 
used at wall nodes where a third upstream point is 
not available. 
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MEANA CONT 



COMPUTE RUY, RVY 
AT WALL POSITION 
-centered 






4 . 2 Subroutine MEANB 


CALLING SEQUENCE: 
INPUT DATA: 
OUTPUT DATA: 
USAGE: 

FUNCTION: 


CALL MEANB 
R, RE, U, V 

RX, RY, REX, REY, UX, VY 

This subroutine is called from MAIN once each 
program time step. 

Linear interpolation between node centers is used 
to define x- and y- components of velocity at the 
node faces. Quadratic interpolation between three 
node centers is used to obtain the values of density 
and internal energy at each of the node faces. The 
signs of UX and VY are used to bias the interpola- 
tion in the upstream direction. A linear interpolation 
is used at wall nodes where a third upstream point is 
not available. 
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FIGURE 4.2 


FLOWCHART OF MEANB 



COMPUTE 

RX(I,J), REX(I, J) 
- centered 
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4 . 3 Subroutine D IFF 


CALLING SEQUENCE: 
INPUT DATA: 

OUTPUT DATA: 

USAGE: 

FUNCTION: 


CALL DIFF (I , J) 

P, U, V, RUX, RUY, RVX , RVY, UX, VY, L, I, J 
DRUUDX , DRVUDX, DRVDY , DRVUDY , DPDX, DPDY, D2UDX2, 
D2VDY2, D2UDXY , D2VDXY 

This subroutine is called from MAIN for each node 
(I,J) at each program time step. 

Difference terms used in the momentum equations are 
computed in this subroutine. 











DIFF CONT 1 
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RETURN 


DIKE OONT 2 


230 

COMPUTE AT BOTTOM HALT, 
DPDX, D2UDX2, D2VDX2 , 
D2UDXY, D2VDXY 
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DPDY, D2UDY2 , D2VDY2 











4.4 Subroutine DIFF2 


CALLING SEQUENCE: 
INPUT DATA: 

OUTPUT DATA: 

USAGE: 

FUNCTION: 


CALL DIFF2 (I , J) 

RX, BY, REX, REY, UX, VY, H, K, 001, DQ2, D03, 

DQ4, L, I, J 

DUDX, DVDY, DRUDX, DRVDY, DREUDX , DREVOY, D2HDX2, 
D2HDY2 

This subroutine Is called from MAIN for each node 
(I,J) at each program time step. 

Difference terms used in the continuity and energy 
equations are computed in this subroutine. Heat 
leak boundary conditions are imposed at exterior 
node faces as prescribed temperature gradients. 
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FIGURE 4.4 FLOWCHART OF DIFF2 
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4 . 5 Subroutine HEATER 


CALLING SEQUENCE: 
INPUT DATA: 

OUTPUT DATA: 
USAGE: 

FUNCTION: 


CALL HEATER (DQHEAT, L3) 

DQHEAT , L3, DRE 

Heater node location is hard-coded. 

DRE 

This subroutine is called from MAIN once each 
program time step whenever DQHEAT > 0. 

The heating rate of the specified nodes is increased 
to represent heater input. 
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FIGURE 4.5 


FLOWCHART OF HEATER 
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4.6 Subroutine BC0UT 


CALLING SEQUENCE: 
INPUT DATA: 
OUTPUT DATA: 
USAGE: 

FUNCTION: 


CALL BC0UT (N, WD0T) 

N, WD0T, R, RBAR, DRU, DRV, DRE, RE, REBAR, L 
DRU, DRV, DR, DRE 

This subroutine is called twice from MAIN during 
ennh program time step: The first time with 

N=l; the second time with N=2. 

The outlet port location is hard-coded. During 
the first call, the fluid velocity at the node 
wall, UWALL, is computed such that the prescribed 
flowrate, WD0T will take place. This velocity is 
used to establish the convection rate of x- and y- 
momentum from the node. 

During the second call, UWALL is U3ed to establish 
the convection rate of mass and internal energy 
from the node. 
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4.7 Function TEMP 


CALLING SEQUENCE: 
INPUT DATA: 

OUTPUT DATA: 

USAGE: 

FUNCTION: 

NON-STANDARD UNITS: 
REFERENCE: 


T(I,J)=TEMP(ERHO, RHO) 

ERHO, RHO, data tables E, T 
TEMP 

This function is called from MAIN for each node 
(I,J) at each time step. 

Linear interpolation in a real-gas data table is 
used to compute a node temperature corresponding 
to the interval energy of the node. 

E(B/lbm) 

The data for these tables was taken from 
Reference 2 and assumes the nominal tank pressure 
of 60 atmospheres (882 psi) . 
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FIGURE 4.7 



FLOWCHART OF TEMP 











4.8 Function BETA 


CALLING SEQUENCE: 
INPUT DATA: 

OUTi I DATA: 
USAGE: 

FUNCTION: 


REFERENCE: 


EO a BETA(TO)*778. 156*32. 2 
TO, data tables E, T 
BETA 

This function is called once at the beginning 
of MAIN during initialization. 

Linear interpolation in the real-gas data table 
is used to obtain the initial specific internal 
energy corresponding to a given initial fluid 
temperature. 

The data for these tables was taken from 
Reference 2 and assumes the nominal tank pressure 
of 60 atmospheres (882 psi) . 
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4 , 9 Subroutine PRESS 


CALLING SEQUENCE: 
INPUT DATA: 

OUTPUT DATA: 
USAGE: 

FUNCTION: 


CALL PRESS (P, R, HH) 

R, HH, PBAR, RBAR, Z, RCONST 
P 

Not used for real-gas problems 

The relative pressure (p-p) is computed in double 
precision from the ideal gas equation with a 
compressibility factor. 
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FIGURE 4.9 FLOWCHART OF PRES 
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4.10 Function RTPRES 


CALLING SEQUENCE: 
INPUT DAT a: 

OUTPUT DATA: 
USAGE: 

FUNCTION: 


P a RTPRES (RR, TO) 

RHO, TO 
RTPRES 

This function is called from MAIN once for each 
node (I,J) at each time step. 

A node pressure is computed using Stewart's 
equation of state (Reference 3) given density 
and temperature. 
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FIGURE 4.10 


FLOWCHART OF A7rP.ES 









4.11 Subroutine BULK 


CALLING SEQUENCE: 
INPUT DATA: 
OUTPUT DATA: 

USACE: 

FUNCTION: 


SUBPROGRAM 
REFERENCES : 


CALL BULK 

P, R, RE, H, PBAR, REBAR, L, NG, NS, NO, NF, N0DES 
P, R, RE, PBAR, RBAR, REBAR, UMIN, HBAR1 , UMAX, 

PC0L, WT 

This subroutine is called periodically during the 
execution of MAIN, normally just prior to a call 
of the OUTPUT subroutine. 

The average tank quantities PBAR, RBAR, and REBAR 
are computed in this subroutine. The relative 
values, P, R, RE are revised according to the 
hange In the average values. The equilibrium 
temperature and the corresponding potential collapse 
pressure are also computed in this subroutine. 

TEMP, RTPRES 
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4.12 Subroutine OUTPUT 


CALLING SEQUENCE: 
INPUT DATA: 

OUTPUT DATA: 
USAGE: 

FUNCTION: 


SUBPROGRAM 

REFERENCES: 


CALL OUTPUT 'T, DT, DTCR, REF1, REF 2, REF3) 

P, R, H, U, V, PBAR, RBAR, T, REF1, REF2, REF3 
Heading LABEL 

This subroutine is called from MIN every DTPR 
seconds. 

The printing of all output data Is controlled by 
this subroutine. The units of PBAR and RBAR 
are adjusted for display and the output heading 
and the data in array LABEL is printed. Also, 
the printing of P, R, H, U, and V by subroutine 
DISPLY is controlled. If an output tape is being 
gen' rated (IUNITW > 0) this subroutine calls the 
tape I/O package to write a record. 

TAPEI0 , DISPLY 
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FIGURE 4.12 


FLOWCHART OF OUTPUT 
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4 . 13 Subroutine DISPLY 


CALLING SEQUENCE: 
INPUT DATA: 
OUTPUT DATA: 
USAGE : 

FUNCTION: 


CALL DISPLY (X, XBAR, CONST) 

X, XBAR, CONST, NG, NS, NO, NF, TIME 
TIME, Y 

This subroutine is called five times during the 
execution of a CALL to OUTPUT. 

The variables P, R, H, U, and V are displayed by 
successive calls to this subroutine from OUTPUT. 
Prior to printing, the reference value XBAR is 
subtracted from each element in X and the units 
are adjusted by the multiplier. CONST: Y(I,J)= 
(X(I,J) - XBAR)*CONST. The rows of data are 
displayed in inverted order to coincide with the 
physical description. 
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FIGURE 4.13 
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4.14 Subroutine CR2TAP 


CALLING SEQUENCE: 
INPUT DATA: 

OUTPUT DATA: 

USAGE : 

FUNCTION : 


CALL CR2TAB (ENDECK, 15, LINES) 

Card images, 15, LINES 

Card images on computer internal unit (15) . 

This subroutine is called from MAIN whenever card 
input is expected. 

The purpose of this subroutine is to document all 
the data input to the program by card. Images of 
the input data cards are stored on the computer 
internal unit (15) and are output on the line 
printer. Card reading is initiated by a CALL of 
this subroutine and is terminated upon reading a 
card beginning with the alphameric word stored in 
ENDECK. This word was selected to be the NAMELIST 
terminator word bgENDb. Data is input through 
NAMELIST from the Internal unit 15. 
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FIGURE 4.14 FLOWCHART OF CIR2TAP 
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4.15 Subroutine RDTAPE 


CALLING SEQUENCE: 
INPUT DATA: 

OUTPUT DATA: 
USAGE: 

FUNCTION: 


SUBPROGRAM 

REFERENCES: 


CALL RDTAPE 

IUNITR, IFILER, IRECR, previously generated tape 
data, 

A-array (data from input tape) . 

This subroutine is called once at the beginning 
of MAIN if an input data tape is used to initialize 
the program. 

The series of operations performed by successive 
calls to the TAPEI0 package is collected here. 

The input tape read parameters supercede these 
input from tape. 

TAPEIfc 
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FIGURE A. 15 


FLOWCHART OF RDTAPE 
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h J 6 SUBROUTINE TAPE 10 


IDENTIFICATION 

Name/Title - TAPEI0 (Tape Input/Output:) 

Programmer /Date - John Prewitt, September 1970 

Organization/Installation - TRW for EP5-MSC 

Source Language - FORTRAN V 


PURPOSE 


Subroutine TAPEI0 performs unformatted tape input /output functions with 
either NTRAN or FORTRAN I/O packages. 


USAGE 


• Calling Sequence 


CALL 

TAPEPS 

(ITYPE, 

IUNIT % 

IFILE, IREC) 


CALL 

TAPERD 

(ITYPE, 

IUNIT, 

1W0RDS , A, 

LSTAT) 

CALL 

TAPEWR 

(ITYPE, 

IUNIT, 

IW0RDS, A, 

LSTAT) 

CALL 

TAPERW 

(ITYPE, 

IUNIT) 



CALL 

TAPE0F 

(ITYPE , 

IUNIT) 



CALL 

TAPEPR 

(ITYPE, 

IUNIT) 



CALL 

TAFECK 

(JSTAT) 





Arguments : 


• Function 
Position tape 
Read tape 
Write tape 
Rewind tape 
Write end-of-file 
Print status of tape 
position 

Check status of NTRAN 
read /write 


Parameter 


Name 

In/Out 

Dimension 

ITYPE 

In 

1 

IUNIT 

In 

1 

IFILE 

In 

1 

IREC 

In 

1 

IW0RDS 

In 

1 


Type 

I 

I 

I 

I 

I 


Description 
Type of data tape 

0 = NTRAN 

1 = FORTRAN 

Physical unit for tape assignment 
Number of files to be skipped 
Number of blocks or physical 
records to be skipped 
Number of data words to be 
transmitted 
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Parameter 

Name In/Out Dimension Type Descriptio n 

A In/Out IW0RDS I Storage area for data wends 

LSTAT Out 1 I Status word fr?: NTRAN read/write 

-1 = Transmission i.3 not complete 
-2 = End-of-filL* for read, 

end-of-tape or drim-file for 
write 

-3 ■- Device error 
■4 e. Transmission abort 
lVwRDS = Number of data words trans- 
mitted when transmission is 
complete 

Status word for FORTRAN read/write 
IW0RDS = Number of data words transmitted 
when transmission is complete 

JSTAT In/Out 1 I Status word for NTRAN raad/wrlte 

• Data In/Out 

None 

• Error Messages 

If an error occxirs subroutine KILLER is called and a walkback is 
generated, 

• Storage 

Coding occupies 774g (508^) locations. Internal data occupies 
142g (98 ^q) locations. 

METHOD 

• Model 

The logic for each of the functions in TAPEI0 is basically the same. 

First the status of the tapes position is updated, next a test is made 
on the type of data tape (NTRAN or FORTRAN) being processed and then the 
appropriate function is called. 

This routine works on all tape or cape simulated I/O devices (see 
restrictions) . 

• Symbol Definition 

None 
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References 

"UNIVAC 1108 FORTRAN V Programmer Reference Manual'', UP-3569 Rev, 1, 

P. N. Bertstresser, T. ^W. Rimkus, "Compute 1 * Systems Bulletin No, 64, 

File and Physical Record Skipping on the Unlvac 1108" , TRW I0C 
P513. 20-37, 20 August 1959, 

H. W. Bryan, "Computer Systems Bulletin No, 7.1, System Supplied Editing 
Routines", TRW IOC 5513,20-60, 2 February 1970. 

RESTRICTIONS 
« Analytic 

IUNIT can have any value between 1-29 except for 5, 6, and 17. 

• Operational 

The tape positioning function does not work for high speed drum (FH232). 

• Other subprograms required 

KILLER 

NTRAN 

QQFIL 

QQREC 

ACCURACY 

Not applicable 

VALIDITY 

A main driver was developed to check all possible combinations of tape 
Jnput/output functions. 

CODING INFORMATION 

• Special Program Constants 

None needed 
« Timing 

The speed is dependent on the amount of information to be transmitted 
and the I/O device used. 
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FIGURE 4.16 
FLOWCHART OF TAPE 10 


DETAILED FLOW CHART 


SUBROUTINE TAPE 10 




6? 












SUBIK/UTINIi TAPEI0 




9/70 


TAPEI0-6 






SUBROUTINE TAPEIH 



TATEI0-7 





SUBROUTINE TAPE 10 



TAPEI0-8 



jRB ROUTINE l'AI'EIfJ 



9/70 


TAPEX0-9 
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SUBROUTINE TAPlilO 



t 

( RETURN 


[ TAPF.CK 


YES 



TEST 

FOR NTRAN 
READ/WRITE ERROR 
JSTAT:<0 


r YES CALL KILLER 
^RROR) ST0P 



9/70 TAPEI(5-10 
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4 . 17 Function OPTD 

CALLING SEQUENCE: ROOPTD (PI, HO) /32.2 

INPUT DATA: PI, HO 

OUTPUT DATA: OPTD (oxygen density) 

USAGE: When initial conditions are generated by the program, this 

function is called once from MAIN to establish the initial bulk density 
and then is also called at each node point at time t=o to establish the 
initial individual node densities. 

FUNCTION: This subprogram employs linear interpolation of data tobies 

obtained from Reference 2 to obtain the density of oxygen at a given pressure 
and temperature, 

REFERENCE: Further documentation of this subprogram is contained under 

Subroutine THERM in Reference 4. 
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APPENDIX I 

COMPUTER PROGRAM NOMENCLATURE 


The following nomenclature list contains the program variable names 
appearing in the storage assignment lists generated for each program 
element. The x's at the right of the page designate in which program 
elements the variable names appeared. The following key Identifies the 
element numbers with the names of the program elements. 


Number 

Name 

0 

MAIN 

1 

MEANA 

2 

MEANB 

3 

DIFF 

4 

DIFF2 

5 

HEATER 

6 

BC0UT 

7 

TEMP 

8 

BETA 

9 

PRES 

10 

RTPRES 

11 

BULK 

12 

OUTPUT 

13 

DISPLY 

14 

R.DTAPE 

15 

— 
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NOMENCLATURE 


VAK jyp UlM UESCHIFTIOim (UMTS) 


A R 3700 MASTER EQUIVALFNCING ARRAY FOR T APF I/C 

A VE R 1U ARRAY CONTAINING AVERAGE STATE VALUES 

BETA R .1 RETURNED SPEC, INT. ENERGY (E) FROM FCTN beta 

C R 1 SPECIFIC HE AT t CP* (UNUSED) (FT -LBF/SLl'G— R 1 

CSTL2 R 1 CONSTANT 1/L? (FT-? ) 

CSTiL R l CONSTANT 1/L (FT-1) 

CSUL? R 1 CONSTANT 1/L?{FT“2) 

CST2L R 1 CONSTANT 1/2L (FT-1) 

CST2L2 R 1 CONSTANT 1/2L2 (FT~2) 

CST4L P 1 CONSTANT 1/4L (FT-1) 

CST4L? P 1 CONSTANT 1/4L2 (FJ-2) 

CTL R 20 PROGRAM CONTROl DATA ARRAY 
UPtiAK R 1 CHANGE IN PBAK 

DPUX R 1 DERIVATIVE OF P WRT V (LPF/FT2/FT) 

DPDY R 1 DERIVATIVE OF P WRT Y (LPF/FT2/FT) 

OGbC R 4 HEAT LErtK BOUNDARY CONDITION INPUT ARRAY 

DOH R 1 INTERNAL HEAT GENERATION ( FT-LPF/SFC/NODF ) 

DOHEAT R 1 INTERNAL HEAT GENERATION ( FT-LPF/SFC/NCDF ) 

D01 R 1 HEAT FLUX - LEFT BOUNDARY (FT-LPF/FT2-SFC > 

QG2 R 1 HEAT FLyX - RIGHT BOUNDARY 

DOS R 1 HEA 1 FLUX - POTT. BOUNDARY 

0Q4 P 1 HEAT FLUX - TOP POUNPARY 

DR R 2UX20 NODE DE.nSiTy RATE (SLUG/ FT3/SEC) 

DRBAR R 1 CHANGE IN RDAK 

URE R 20X20 NODE INTERNAL ENERGY KATE ( FT-LBF/FT3/SEC J* 
DREBAR R 1 CHANGE IN RERAP 

DREUDX R 1 DERIVATIVE OF PE U wRT X ( FT-LPF /PT3/SFC ) 

DRtVUY 9 1 DERIVATIVE OF PE V «Rf Y ( FT-LPF/FT3/SFC ) 

URU R 20X20 NODE X-MOMENTUN RATE ( SLUG— FT /FT3— S£C2 ) 

URUDX R 1 DERIVATIVE OF PU WRT X ( SLUG— FT/FT 3- SEC/ FT ) 

URUUUX R 1 DERIVATIVE OF RU U W^T X ( SLUG-FT/FT3- C EC 3 

DRUVDY R 1 DERIVATIVE OF Pu V WR T Y ( SLUG-rT/FT3— SEC2 ) 

DRV R 20X20 NODE Y-MOMENTUN RATE ( SLUG-FT/FT 3- , =Er2 ) 

URVDY R 1 DERIVATIVE OF R V wRT Y ( SLUG-FT/FT 3-^EC/FT) 


XXX X X X X X XX XXX XX X X X X X X xxxxxxxx 





X 

X 

X 


012345678^0 12 34S' 
111111 


NOMENCLATURE 


PAGE 2 


VAN 1YP ulM DESCRIPTION (UNITS) 


URVUUX R 1 DERIVATIVE OF PV U WPT X l SLUG~FT/FT3-.SEC2 ) 

ORVVUY R 1 DERIVATIVE OF Pv V wPT Y l SLUG-FT/FT3-SEC2 ) 

OT P 1 PROGRAM TIME STEP (SFC) 

UTPR R 1 PRINT TIME INTFRVAL (SFC) 

Oil R 1 PROGRAM TIME STEP DISPLAYED (SFC) 

DUUX R 1 DERIVATIVE OF U WRT X ( FT /FEC2 ) 

DVUY R 1 DERIVATIVE OF V WRT Y (FT/SEC2) 

02HDX2 R I ? ND DFRIVATIVF OF h wRT X (R/FT?) 

U2HDY2 R 1 2 NO DERIVATIVF Of H wR T Y ( P/FT2 ) 

02UDXY R 1 9 . ND DERIVATIVF Op U WR 1 X,Y (FT/SFC/FT2) 

D2UDX2 R I 2 NU DERIVATIVE OF U WRT X (FT/SFC/FT2) 

□2UDY? Pi 2 ND DERIVATIVE OF U WRT Y (FT/SEC/FT2) 

U2VDXY R 1 ? NU OErJ.VATIVF OF V WRT X»Y (FT /SEC/FT2 ) 

G2VDX2 R 1 2 NU DERIVATIVF Op V wRT X ( FT /SFC/FT2 ) 

w 02VDY2 R 1 ? NU DERIVATIVF OF V WRT Y ( FT /SFC/FT2 ) 

E R 3U SPECIFIC INTFKNaL FNFRGY TABLE ( RTU/LBM ) 

EIlV R 1 SPECIFIC INTERNAL ENERGY IN ( BTL/LRN ) 

EMJECK .1 ALPHAMERIC tfrminatop FOR CARD INPn 

ERHO £ \ INTtRNAL ENERGY (FT-LBF/FT3) 

EO R 1 TNIflAL TANK SPEC. IMTFRNAL FNFKGY ( FT-LPF/SLUG) 

gx r l acceleration v— Component <ft/sec2) 

GY R 1 ACCELERATION Y-ComP^NFNT (FT/SEC2) 

H R 2UX20 NODE TEMjPERATUFE (P) 

HBAR R 1 AVERAGE TEMpFRA TURF (R ) 

hPARX R 1 AVERAGF TEMPFKATUrT FOR OUTPUT (R) 

HH n 1 DOUBLE PRECISION NOUP TEMPFRATI RF (R) 

HMaX R 1 MAXIMUM NODE TEMPERATURE C R } 

HM1N R 1 MINIMUM NODE TFMPERATURE t R I 

HO R 1 INITIAL TANK TFMPERATURE (R) 

i II NODE INDEX X-DIRECTION I TO 20 

I 11 LAST INDEX USED In T-E TABLES 

IFILER I 1 TAPE FILE NUMbFR - KFAP1NG 

XFILEW I 1 TAPt FILE NUMBER - WRITING 

1RECR I 1 TAPt RECORD NUMBER - READING 


x x 


PROGRAM ELENENT 

min' 

ui 23_45b7aR0 1 2345 

X X 
X X 
X 
X 
X. 

X 

X 

X x 
X X 
X V 
X, x_ 

X X 
X X 
X x 
X X 

XX 

X 

X 

X 

X 

X 

X 

X XX XX 

XX 
X 
X 

X 

X 

X 

XXX XX XX 

xx 

X X 

X 

U1234S670QUI2345 
11 1111 



'MOME.IMCLA T U K E 
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VAh 1YP UTM DESCRIPTION (U^ITS) 


IRECw T 1 TAPE KECORO NUMBER - WRITING 

i TAPE I 6 TAPE CONTROL U A T A ARRAY 

ITERM I I TERMINATION FLAG FOR NEGATIVE DENSITY 

1UNITR I 1 TAPE UNIT NUMBFR - RFADING 

iUNITW I 1 TAPE UNIT NUMBER - WPITING 

11 II STARTING i-INDFX AT ROW J 

12 II STOPPING I-INOFX AT ROW J 

IS II ALTERNATE unit NUMBER Fur CARD INPUT TPRU CH2TAP 

J II NODE INDEX Y-OTRECTION l TO 20 

JF1LER I 1 TEMPORARY TAPE INPUT FlLF NUMBFR 

JRECR I 1 TEMPORARY TaPE INPUT RECORD NUMBER 

JUMTR I 1 TEMPORARY TAPE -NPuT UNIT NUMBFR 

J1 II STARTING U-INDFX AT COLUMN I 

J2 II STOPPING J-INDFX AT COt UMN I 

K R 1 THERMAL CONDUCTIVITY (FT-LPp/FT-SEC-R ) 

K II INVeRTFU ROW INDEX 

L R 1 NODE DIMENSION l FT ) 

LABEL R 10 OUTPUT LABEL DAT* ARRAY 

LIMITS I 50 INDEX LIMITS DFFlNING TANK CONFIGURATION 
LINES I 1 NUMBFR OF PRINT LINES FOR CARD INPUT DISPLAY 

LSI AT I I STA1US OF TaPE OUTPUT 

L3 R 1 VOLUME OF NODE (FT3) 

M II NUMBFR OF ENTRIES IN T/E TAbLE^ 

mi ii no - loop Index 

N 11 NUMbFR OF ENTRTtS IN T-E TABLES 

NF I 1 maximum index of tank configuration 

NG I 2U ARRAY CONTAINING STARTING INPICIFS 

NODES R 1 TOTAL NUMbER OF NODES IN TANK CONFIGURATION 

ns i 2u array containing endtng indicies 

NO II MINIMUM INDEX CF TANK CONF TgURA T TON 

P R 2UX20 NODE PRESSURF RELATIVE TO PBAH (LBF/FT?) 

N8AR D 1 REFERENCE PRFSRURE - UPDATFU A < AVFHAGF t LP/FT? ) 

PBARI P 1 AVERAGE FLUID PRESSURE OUTPUT (LPF/IN2) 

PCOL R 1 POTENTIAL COLAPSF PRFSSURt (LdF/IN?) 


program elenent 
amii’ 
ol23456?890 12345 

x 


X 

X 

X 

X 

X 

X 

xxxx 

X 

xxxx 

X 

X 

XXX XX 

X X 


X 



■ 

X XX 

X 

X xx 

X 

X X 

X 

X xx X 

X 

X 

Xx’ 

XVXXX 

X x 

X 

X 

X 


X 


X 


V 


XXX 

X x 

xxxvx 

X x 

X 

X x 

xxxxx 

X X 

XYX 

x x 

X X 

Xx 


X X XX 
X 
V 


0123456789012345 



nomenclature 
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VAR 

t YP UTM 

DESCRIPTION (UNITS) 

PROGRAM 

ELEMENT 







ii mu 





U1 23lt5til?0Ppl 2345 

PFLAG 

R 

1 

PRESSuRt FLAG FOR HEATFR CYC! E OPERATION 

X 



PHI 

R 

i 

PRESSURE UPPER LIMIT FOR HFaTER CYCLF (LPF/FT2 ) 

X 



PLO 

R 

1 

PRESSURE LOWFK lIMII FOR hFAT£R CYCLF (LPF/FT2 ) 

X 



PP 

P 

1 

nOlJOLE PRECISION PPE^SUkE (LPF/FT2) 



X 

PREF 

R 

i 

TANK PRESSURF FOR HEATFR CYCLE (LbF/FTP) 

X 



PROP 

R 

lu 

fluid property data array 

X X 

X 

X Xx 

PT 

R 

I 

SUM OF PRESSURE VARIATIONS ( LdF/FI?) 



X 

PO 

R 

1 

INITIAL TANK PRESSURE ( LPF/FT 2 ) 

X 



PI 

R 

1 

NODE PRESSURE FOR CALL TO OPTO ( LPF/IN2 ) 

X 



K 

R 

2UV20 

NODE DENSITY RElATTVF TO REAR (SLUG/FT3) 

XXX 

X 

XX 

RA1FS 

R 

lbflO 

array contain© fluid state data 

XX xxx 

X. 

xx 

RBAR 

D 

1 

REFERENCE OENSTTY - UPDATED As AVERAGE ( SLUG/FT 3 ) X X 

X 

X XX 

KBAP1 



AVERAGE FLUID DENSITY OUTPUT <LdM/FT3) 



XX 

RCONST 

R 

1 

GAS CONSTANT, P (UNUSED) ( FT-LPF/SL UG-P ) 

X 


V 

A 

RF 

P 

20X20 

NODE INI. ENERGY RELATIVE TO REdAR (FT-LPF/FT3) 

x xxx 

X 

Xx 

KFbAR 

0 

1 

REFERENCE InTEPinAl ENERGY CFT-LdF/FT3) 

X 

X 

X 

KEbARl 

R 

1 

average internal energy 



X 

REPG 

R 

1 

CENTFR LOCATION' IN TANK FOR ACCELERATION HEAD 

X 



RFJ 

R 

1 

SUM OF INTEKNAI ENERGY VARIATIONS (FT-LBF/FT^) 



X 

REX 

R 

21X21 

FACt-CENTEREn VALUES OF RE IN X-DIRECTTQN 

XXXX 



KEY 

R 

21X21 

FACL-CENTERED values of re in y-direction 

XX xx 



RHO 

R 

1 

node density c=;lUg/ft3) 

X 



K1 

R 

1 

real INDEX 1 

X 



RJ 

P 

1 

REAL INOEX J 

X 



RR 

D 

1 

roudlc PRECISION DENSITY (SLUG/FT3) 

X 



RT 

R 

1 

SUM OF DENSITY VARIATION^ <SLUPS/FT3> 



X 

RTPRES 

R 

J 

RETURNED PRESSURE FROM FUNCTION RTPRFS (LbF/FT ?) 

X 


X X 

RU 

R 

20X20 

NODE MOMENTUM X-DIPECTION (Si UG-FT/FT3-SEC) 

XX XX 

X 

X 

RUX 

R 

21X21 

FACE-CENTERED VALUES OF RU IN X-D1RECTI0N 

xxxx 



RUT 

R 

21X21 

FACE-CFiMTERED values OF RU in y-direction 

xxxx 



RV 

n 

2UX20 

NODE MOMENTUM Y-DIRECTION ( SLUG-FT/FT3-SFC ) 

XX xx 

y 

X* 

A 

KVX 

R 

21X21 

FACE-CFnTERED VALUES OF RV IN X-D1RECT TON 

xxxx 



RVY 

R 

21X21 

FACt-CENTERED VALUES OF RV IN Y-DIPECTTON 

xxxx 



RX 

R 

21X21 

face-cfntered values of k in x-directton 

xxx 



RY 

R 

21X21 

FACE-CENTERED VaLUFS OF R IN Y-DIRECTION 

xxx 




U1 2^4 ^678^01 2345 
llltll 



NOMENCLATURE 
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VAR IYP UTM DESCRIPTION (UNITS) 


HO R X INITIAL FLUID DENSITY (SLUG/FT3) 

scale p i scale factor 

state P 2000 ARRAY CONTAINING FLUTU STATE DATA 
STOPFG R 1 END OF RUN FLAG 

T R 1 PROGRAM TIME (CEO 

T R 30 TEMPERATURE TABLE (H) UN ASCENDING ORDER) 

IFMP R 1 RETURNED NODF TE^pera fURF FROM FUNCTION TEMP 

TIME R 1 PROBLEM TIME (VlNuTES) 

TPK R 1 PROGRAM TIME FOR NEXT PRINTED OUTPUT C^EC) 

T SCALE R 1 SCALFD PROBLEM TlMF (MINUTES) 

TSTOP R 1 PROGRAM END TIM£ (SEC) 

TO R 1 PROGRAM S I ART TI M E (SEC) 

U R 20X20 NODE VELOCITY V-DlRECTION (FT/SEC) 

UWALL R 1 WALL VELOCITY A T EXIT PORT (FT/SFC) 

UX R 21X21 FACE-CENTERED VALUES OF U IN X-DIRECTION 

V R 20X20 NODE VFlOCITY Y— DIRECTION (FT/SEC) 

vsc r l fluid absolute viscosity (ldf-sec/ft?) 

VY R 21X21 FACE-CENTERED VALUES OF V IN Y-DlREC7;oP 

WDOT R 1 MASS FLOWRATE FROM FLUID VOLUMF (SlUO/FEC) 

wT R 1 SUMMFD vvEIGHT OF FLUTU (LBM] (OUTPUT) 

WTM R 1 FLUID MOLECULAR WEIGHT (UNllSFQ) (LP/LB-MOLE) 

Y R 20X20 TRANSFORMED OUTPUT DATA ARRAY 

£. R 1 FLUID COMPRESSIBILITY FACTOR (UNUSED) 


PROGRAM ELEMENT 
1111 if 
O’l 2345678*301 2345 

X 

X 

XX X X XX 

X 

X 

XX 

(R) X XX 

XX 

A 

X 

X 

X 

XXXXX X 

X 

xxxx 

xxxxx X 

X 

xxxx 

X 

y 

x X 

X 

X xx 

012345U7890 12345 
11 1111 


APPENDIX II 

COMPUTER PROGRAM LISTING 

This Appendix contains a listing of the GNAT Computer Program followed 
by results from a sample computer run. The input control cards show t.ie 
specified initial state of the fluid and the fluid properties used. Computer 
output on magnetic tape was specified by the ITAPE input data. The standard 
computer printout is shown at initialization and at the first print time 
specified. 
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BZ HUN ''?3|**EFSiS‘l8(3|9BpsF2DA*Ct2>l 


HEJNHJLLER ft 


SH ttSG F|tF RES* TAPE l FHS32 | FST«N 


27 JAN ?1 


8 X r 3«*7| 


2 7 JAN T ‘ 


This page shows typical control cards which are necessary to execute the GNAT 
program on the Univac 1108 system at NASA/MSC. 


IS* 33125 
19 * 33 * 2 * 



5. XUT CUR 

It TH* X 
2 t ill I 


27 JAN 7 1 


1 1 1 33 1 25 • >2 
11(33:25.858 

11 ( 33 : 25.840 


EhO OF FILE *>. UNIT X 
3t TRt X 
1. TOC 


ih;33;3i.bi2 
Mi3X; JR.B1S 


FlEmEnT table 


BlKi 


StHBOLIC 

23 

DEC 

70 

23 Io5 io? 

0 

Ql N3E47q 

n 

(66 

BLK2 


STHBOLIC 

23 

DEC 

70 

23:osu- 

0 

01HH331H 

11 

l ®3 

CROSS 


STHBOLIC 

23 

DEC 

70 

23 J Q5 ; | Y 

0 

01HS03J 4 

11 

IS 

CR2TAP 


STHBOLIC 

23 

DEC 

70 

23;os;2l 

0 

01H5IS54 

11 

100 

CSSTTH 


STHBOLIC 

23 

DEC 

70 

23:05:25 

0 

01H5H3H4 

(1 

IQ? 

EQTANK 


symbolic 

23 

DEC 

70 

23;os:3» 

0 

0 1 H 475QH 

(1 

317 

HdTPP 


symbolic 

23 

DEC 

70 

23:05(30 

0 

01500232 

11 

C3 

HPTO 


SYMBOLIC 

23 

DEC 

70 

23;os;5<' 

0 

0(502111 

(1 

275 

HPPTT 


SYMBOLIC 

23 

DEC 

70 

23:o&:o7 

0 

01512054 

11 

*2 

HpTCV 


$\HB0lIC 

23 

DEC 

70 

23;oa;io 

0 

0)S| 1252 

(1 

31 

HPTcP 


sy. <s“lic 

23 

DEC 

70 

23;o *: 12 

0 

0(5(53(1 

11 

3? 

HpTH 


symbolic 

23 

DEC 

70 

Z3;oa:i3 

0 

0 1 S 14354 

11 

1*0 

HPTPT 


SYMBOLIC 

23 

DEC 

70 

23 j 04| 20 

0 

01523304 

11 

• 5 

KfTV 


SYMBOLIC 

CJ 

DEC 

70 

23:06(23 

0 

0(525551 

(1 

266 

KpTTC 


SYMBOLIC 

23 

DEC 

70 

23:qa;32 

0 

01531770 

11 

22? 

OOT 


symbolic 

23 

DEC 

70 

23 : 04 : 1 1 

0 

0151417 4 

(1 

lie 

OpTpT 


STHBOLIC 

23 

DEC 

70 

Z3(04(*ia 

0 

01550522 

11 

B3 

PTHi, AT 


STHBOLIC 

23 

DEC 

70 

23;o4;si 

0 

0 1 552731 

(1 

l ®D 

ST^P 


symbolic 

23 

DEC 

70 

23104150 

0 

0(557441 

11 

11 

THERH 


StHBOLIC 

23 

DEC 

7C 

23 IQ4 :SY 

0 

0154(031 

11 

371 

tpocb i 


symbolic 

23 

dec 

70 

Z3;07j(2 

0 

01573(14 

11 

»27 

TpOCB2 


SYMBOLIC 

23 

DEC 

70 

23 J 07 J *f S 

0 

0(421130 

11 

2Q4 

H 1 IfJ 


symbolic 

2H 

DEC 

70 

0Y;l7;|0 

0 

0(432(31 

11 

301 

B I NSER 


STHBOLIC 

21 

DEC 

70 

oy;i7;3I 

Q 

□1412371 

1 1 

101 

HF*9 


SYMBOLIC 

2*1 

DEC 

70 

o»:i7:32 

0 

0(415202 

ll 

106 

CROSS 

CODE 

relocatable 

*0 

DEC 

70 

u:3s:i7 

1 

0(4501 ( 4 

21 

1 








0 

01450(14 

n 

*i 

HPDTT 

CODE 

relocatable 

20 

DEC 

70 

u;3si*i2 

l 

0(450234 

21 

1 








0 

01450244 

11 

t 

step 

CODE 

relocatable 

20 

DEC 

70 

tl;35tS3 

1 

0(450114 

21 

1 








0 

01450174 

|1 

1 

HA IN 

COOE 

relocatable 

23 

OEC 

70 

zajaiil* 

l 

0(450544 

34 

1 








0 

0(450432 

11 

70 

B1.K1 

CODE 

relocatable 

23 

dec 

70 

23 : OH I J B 

1 

0)452554 

21 

( 








0 

014524^4 

H 

101 

BLR2 

CODE 

RELOCATABLE 

23 

DEC 

70 

23;oh : t? 

1 

0(4551(1 

21 

1 








0 

0(455111 

|1 

11? 

CR2TAP 

CODE 

RELOCATABLE 

23 

DEC 

70 

23 ;oh ; zi 

1 

0 i 64 q A14 

21 

1 








0 

01440474 

11 

(1 

CSSTTh 

CODE 

relocatable 

23 

PEC 

70 

23(01121 

I 

0(44(202 

34 

1 








0 

0(44(214 

H 

61 

ebtank 

CODE 

RELOCATABLE 

23 

DEC 

70 

23 : 01 : 2 ? 

I 

0(443014 

34 

l 








0 

014431 (2 

»Y 

B3 

HdTPP 

CODE 

relocatable 

23 

OEC 

70 

23:oh;-31 

1 

0(445321 

21 

l 








□ 

01445351 

11 

B 

HPTCP 

CODE 

relocatable 

23 

DEC 

70 

23:DH:32 

I 

0(445531 

21 

1 







0 

0(645541 

11 

3 

HPTCV 

COOE 

relocatable 

23 

DEC 

70 

:3:uh ;33 

1 

0I66S634 

21 

1 










0 

0 I **5*C* 

1 H 

3 

MpTO 

COOE 

relocatable 

23 

DEC 

70 

23;or;3« 

l 

0 1 6*S?-:Q 

2H 

1 








□ 

016*5770 

IH 

1HI 

KpTH 

C. IDE 

RELOCATABLE 

23 

DEC 

70 

23:0‘C3 I ? 

1 

nj*7i*s* 

2H 

1 








0 

01*7170* 

1 H 

*2 

HpTPT 

CODE 

RELOCATABLE 

23 

DEC 

70 

23 JQSJHO 

J 

0 1 673H52 

2H 

i 








0 

0U73507 

n 

? 

OOT 

CODE 

relocatable 

23 

DEC 

70 

23jO‘t;R3 

1 

01*73700 

2H 

1 








0 

D t *7373 0 

l H 

hB 

HpTTC 

CODE 

relocatable 

23 

DEC 

70 

23;0HS*|6 

1 

01*75170 

2H 

1 








0 

□1*75220 

11 

7* 

MPTV 

CODE 

RELOCATABLE 

23 

DEC 

70 

23;oh;h» 

l 

01*77720 

2H 

1 








0 

01*77750 

IH 

100 

OPTPT 

Code 

RELOCATABLE 

23 

DEC 

70 

23;oh;5o 

l 

□ 1 7Q25N0 

2H 

1 








0 

0 1 70Z570 

IR 

ft 

PTHEAY 

COJ„ 

relocatable 

23 

DEC 

70 

23 iOH 153 

1 

01702750 

2H 

1 








0 

01703000 

n 

*0 

therm 

COOE 

relocatable 

21 

DEC 

70 

23 I OH I Si 

l 

0 l 70HS 10 

72 

1 








0 

0 1 70H*20 

IR 

55 

YpOCBt 

CODE 

relocatable 

23 

OEC 

70 

23:05 JDO 

1 

0)70*222 

2R 

| 








□ 

0170*252 

IR 

283 

TPOCBU 

CODE 

relocatable 

23 

DEC 

70 

23S06J02 

I 

0 ] 7 1 *OHH 

2H 

1 








0 

0 1 7 | *Q7H 

IH 

35 

BINSER 

CODE 

RELOCATABLE 

2H 

OE.' 

70 

0YI|7;3l 

t 

D17170H* 

2H 

i 








0 

0171707* 

1 R 

13 

HF#8 

COOE 

relocatable 

2H 

DEC 

70 

*19 1)7132 

1 

D 1 71 73*H 

2R 

l 








0 

0I717H1H 

JR 

f 


entry point t* 6 le 


B|Ns£r 

l B I NSEr/CODE J 

1 000137 

cross 

KROSS/CoDEI 

1 

00003* 

Cr2t*P 

( CkZTAP/ CODE) 

l 

00013* 

C5STTM 

(CSSTTM/COOE) 

1 000751 

EQTANr 

(E0TANK/COPE) 

1 

0013*6 

H0TPP 

( HOTPP/CODE 1 

1 

000071 

HF*fl 

(Hf*0/CODt) 

1 0080*5 

MPOTT 

jH e OTT/CODE| 

1 

000071 

HFTcP 

(MPTCP/CODE) 

l 

00001* 

HPTcV 

(MPTqV/COOE 1 

1 00001* 

hptd 

IhPTO/COOEl 

l 

oooshn 

HPTh 

ihpTh/ CODE l 

1 

000557 

mptft 

1 HBTp T/COOE 1 

1 00i0*H 

MPTTC 

IHPTTC/COOEI 

1 

OOOSqH 

hpTv 

IHPTV/CODEI 

1 

ODOR* | 

OOTIC 

(00T/COD£t 

1 U0060H 

OOTIT 

t OPT /code ) 

1 

000*13 

OOTP 

* OOT/ CODE 1 

1 

000575 

OPDCP 

| IHErM/CODE ) 

1 000703 

"°0CV 

(THERM/COOE} 

1 

0Q0*H 3 

OPOH 

ITHERM/CoOEI 

1 

□00*03 

OPDJC 

J TMERM/C0DE1 

1 0005H3 

JPD1T 

(THERH/cOOE) 

1 

000503 

OPDT 

( THERM/CODEt 

1 

0D07H3 

OPT CP 

1! THERH/CODE) 

1 000723 

OPTCY 

(THERH/cOOE) 

1 

000**3 

OPTo 

(TmERM/COOE) 

1 

000H*3 

OFTh 

therm/cooe) 

1 000*23 

OPTIC 

1THERM/C00E1 

1 

0005*3 

OPT jT 

1THERH/C0DE> 

l 

000523 

OPTpT 

(OPTPT/COOE) 

I 000052 

OpTTC 

t THERH/COOE 1 

1 

00 1003 

OPTv 

ITHERM/CoSE) 

l 

□007*3 

PTHE*T 

IPTHEAT/CODEI 

1 000*17 

STEP 

(STEP/CODE’ 

1 

000023 









SLOCK 

TABLE 






blkpth 

| BLK] /CODE ) 

3H BANK 

2 oeperden r 



SP«EAT 

1BLK2/C0DE1 

3H BANK 2 

DEPENDENT 

TPCB 

XPOCBI/CODEj 

3H BANK 

2 DEPENDED 



TPCB 

ITPOCBZ/CODE, 

3R BANK 2 

depe 

NDENT 


cobol library Table empty 
pro.e»ure name table empty 


end CUR LCC I102-Q03B LS 



9 ROR HA|N. HA IN 

UNJvAC I 10® FOR7RAN V L£yEL 220* 00 |B F501BH 
This compilation aas done on 27 jan 71 at nj33;3i 


onatooo’ 


27 JAN 7 1 


11133T 31*1*0 


MAIN PROGRAM 

STORAGE USED I BLOCK 1 NAME, LENGTH) 


L J 1 *Co£>E 001217 

OBOD «DATA 000137 

0002 <BLaNK 000000 

0003 DATA 0071*1 

0PC1 D I 000011 

OOOS 02 DOOOIO 

000* C5TS 00000* 

0007 trnsht ooooos 

0020 OSTATE 003)00 


external references iglock, nahei 


0011 

crztap 

0012 

ROTAPE 

0013 

taper* 

0011 

TAPEPS 

0015 

TaPEPR 

Oil* 

beta 

0817 

OPTd 

0020 

rtpres 

□ •21 

heana 

0822 

0 IFF 

0823 

BCOUT 

002*4 


002S 

0IFF2 

002* 

heater 

0027 

temp 

0030 

BULK 

0031 

output 

0032 

TAPEQF 

0033 

NRNLS 

0031 

N*DUS 

0035 

N I 0 I S 

003* 

NIOZ* 

0837 

NSTQPs 


STORAGE ASSIGNMENT FOR VARIABLES IBL0CK, TyPE, RELATIVE LOCATION, NAME) 


0001 

000001 

IDOL 

0001 00003 1 

IOSL 

3001 000123 

USL 

0001 


Q0011 3 

1 7 OL 

OOOi 0000*1 

I *0L 

OoO I 

000117 

200L 

0001 200071 

20" 

0001 000117 

21SG 

0001 


000212 

2l*G 

OOOi 000332 

27?G 

000 , 

0003S ) 

277s 

ODD) QOOS20 

3', G 

0001 0C0531 

3*1 G 

0001 


000*3 * 

3***G 

OOOI 000*52 

2/lG 

oooi 

QQQ7ZQ 

1C1G 

0001 0007)1 

ills 

0001 001013 

127fi 

0001 


03 l 0*5 

1316 

OOOI 0011*1 

180l 

0001 

001201 

SOOL 

0001 00121* 

SSOL 

0000 000070 

7 lOF 

0003 

R 

000000 

A 

1)003 000)32 

AVE 

oou 

R 000000 

beta 

0000 K 000021 

c 

000* R COOOOO 

CST1L 

000* 

ft 

000003 

CSTIlZ 

uOU* R 00000 1 

cstzl 

000* 

R P00001 

CST21.2 

000* R 000032 

CST1L 

000* R ooooos 

CST1L2 

0003 

R 

DODQOQ 

ctl 

0001 a QOODOO 

BPOX 



OOOH 

H 

oooooi 

OPDY 

0000 

R 

000005 

OCBC 

oooo 

R 

00004& 

DnH 

oooo 

R 

000027 

D8HEAT 

0007 

ft 

000001 

081 

□ OQ7 

R 

000002 

0q2 

0007 

R 

O00QQ3 

093 

□ 007 

R 

000004 

084 

0010 

R 

□oaooo 

ds 

0010 

R 

000620 

ORE 

oous 

R 

000002 

OREOpX 

COOS 

R 

□00003 

PREVoT 

0010 

R 

00 1440 

DRU 

0005 

R 

OOOOOO 

DRUB* 

□ 004 

R 

000010 

DRUUOX 

0004 

R 

OOOU13 

oruvdy 

0010 

R 

0022*0 

DRV 

0005 

R 

□00001 

ORVOT 

□ 004 

R 

ODOO 1 2 

drvdd* 

0004 

R 

□□00 1 1 

ORVYDY 

0003 

R 

OOOODl 

dt 

0003 

R 

000003 

dtpr 

0003 

R 

000037 

OTI 

□ 005 

R 

Q0OU04 

dudx 

0005 

R 

000005 

OVOY 

COOS 

R 

ooooo* 

D2HDX2 

0005 

R 

000007 

D 2 HDY 2 

0004 

R 

coooa* 

P2U07Y 

0004 

R 

000002 

02uDX2 

3004 

H 

000004 

02U0T2 

scat 

R 

000007 

02V0XY 

0004 

R 

000003 

0 2 V0X 2 

0004 

R 

000005 

D2VdT2 

OOOO 

R 

□00012 

ENDEC* 

OOOO 

R 

000046 

ERHo 

0000 

R 

000032 

EQ 

DQOO 

R 

000030 

sx 

oooo 

R 

oooo-M 

GY 

0003 

R 

0Q16Q4 

H 

□ 000 

0 

000001 

HH 

0003 

ft 

030034 

HO 

0000 

l 

000015 

I 

0003 


oaoo|6 

ifiler 

□ 003 

| 

OOOO 2 1 

IF! Le* 

000$ 


000017 

|f?Ec J 

t ON ] tr 

i> 2 1 

0003 

t 

000022 

1 REC* 

0003 

i 

oono ; 

ITaPE 

OOOO 

1 

000042 

1 TERH 

0003 

I 

00001$ 

jOniTR 

0003 

I 

I 

000020 

0000 

I 

□oooll 

1 i 

oooo 

I 

□OUu< j 

IS 

OOOO 

I 

□00020 

J 

□ OOU 

I 

000016 

Jl 

OOOO 

000017 

0007 

R 

oocnoo 

K 

0003 

R 

00003 | 

E 

0003 

R 

000036 

LABEL 

0003 

000050 

UNITS 

OOOO 

1 

000014 

LINES 

0000 

R 

oooooo 

E3 

Q003 

1 

R 

000121 

Nf 

0003 

I 

000050 

NO 

0003 

R 

000122 

NODES 

0003 

1 

000074 

NS 

0003 

i 

000120 

HO 

0017 

□OOOOQ 

OPTO 

0003 

R 

000144 

P 

□ 003 

D 

000132 

pbaR 

0003 

s 

00003$ 

PFL*s 

□ DOO 

R 

800043 

PFLAGl 

0000 

R 

000022 

PHI 

OOOO 

R 

□00023 

PLQ 

□ 000 

R 

000Q44 

pRef 

□ 003 

R 

000024 

PROP 

0000 

R 

000021 

PO 

0000 

R 

□00033 

Pi 

0003 

R 

000764 

R 

□ 003 


004q64 

rates 

0003 

0 

000134 

Rbar 

0003 

R 

800032 

RCONST 

0003 

R 

004704 

RE 

0003 

0 

000126 

R£BA r 

oooo 

R 

00003$ 

REFG 

GOOD 

R 

000040 

rho 

0000 

R 

00003d 


ooca 

R 

000037 

RJ 

ODOO 

D 

000003 

RR 

0020 

R 

OOOOOO 

RTpRES 

0003 

R 

005524 

RU 

0003 

R 

006344 

RV 

oooo 

R 

000034 

RD 

0003 

R 

000014 

Scale 

0003 


000144 

state 

0003 

R 

000023 

ST0PF6 

0003 

ft 

800036 

T 

0027 

R 

OOOOOO 

TEMP 

good 

R 

000041 

t»*r 

0003 

R 

000040 

TSCALE 

□ 003 

ft 

0000D2 

TSTOP 

0003 

R 

ODOOOO 

TO 

0003 

R 

002421 

U 

0003 

R 

□03244 

V 

OOOO 

R 

000025 

VSc 

OOOO 

R 

000026 

• dot 

0003 

R 

000024 

*Ti1 

0003 


000030 

z 














00100 

1 * 

c 


GNATODIO 

00100 

2* 

C 

P J HEINHILLER r R * HH H21 TEL s*l 3133 a 2701 

fcNATOOl 1 

ODiOO 

»• 

c 


GNaTOO 1 2 

00100 

*• 

C 


GNAT00I3 

oaiai 

5* 


REAL Li Ki NODES. LAbELi Li 

SNaT0014 

00103 

66 


COH'ON /DATA/ A ( 3700 1 

GNaTOO 1 5 

cnio4 

7* 


COMMON /Dl/ DPDXiPPDYiD2!J0X2.d2VDX2iD2U0Y2i02V0Y2i 

GNATOO 1 6 

004 04 

8* 


1 D2U0XY .D2VIJXT .DR080X .DRVVdTiDRVUdXiORUVOY 

GNAT00I7 

□ 0105 

*• 


COMHUN ff)2/ DRUOXtORVOYiDREUDX.DREVOY.DUOX.DYCTi0 2 HOX2,02HDT2 

GNaTOPIS 

00106 

1B» 


CONMON /CSTS/ CSTlLiCST2LiC5T4LiCSTlL2,C5T2L 2 rCST4L2 

GNaTDOI* 

□ 0107 

1 1* 


COMMON /TRNSMT/ K .DQ 1 i MRE i0 93 , D34 

GNAT0Q20 

00110 

12* 


COMMON /OSTATE/ CRI20.20). DRe(Z0i 20), 0RU120.Z0), 0RV(20 t 20) 

GNAT0021 

00111 

13* 


equivalence 

SNA . 0022 

001 11 

14* 


I Ultl.CTU 

GNAT0023 

00111 

1ST 


2 t ( A l 2 1 1 , PROP 1 

GNAT 0024 

001 11 

16* 


3, (A (31 1 , LABEL) 

SNAT00Z3 

001 1 t 

11* 


3, (At 4] ) .LIMITS) 

GNAT0QZ6 

001 11 

>*• 


4, ( A ( * l 1 i AVE1 

GKAT0027 

00 I l 1 

1*6 


5, UIIOI 1 iSTATEl 

GNAT0C2S 

0011 1 

20* 


6, (A (21011 .RATES) 

GNAT002? 

00112 

21* 


DIMENSION AVEI101) RATES ( 1 6 r .O ) i STATE 1 2000 1 

GNAT0030 

0011$ 

22* 


DIMENSION P(Z0,20) ,R( 20,20 I.H 120,20 1 ,Ul 20,20 l.V 120.20) 

5NATD031 

00114 

23* 


DIMENSION RE! 20 . 20 ) , H0(20i RV(2g,2D) 

GNAT0032 

Dans 

246 


EOUU'ALENCE 

GNAT0033 

00115 

25* 


1 ( ST ATE 1 11, VI 

GNaT0034 

00 115 

28* 


2, (STATcnOl | ,Ri 

SNATQ035 

00115 

2?6 


3, (STATEIBOll iH) 

GNATOO 56 

00115 

28* 


4, ( ST ATE ( 1 20 1 1 .U 1 

GNA T 003 7 

□ 0115 

2** 


5, (STATEI 16011 ,V) 

SNAT003B 

00115 

30* 


7, (RATES (401). RE 1 

GNAT003? 

00115 

31* 


B, ( RATES ( P 01 ) *RU 1 

CNATC04D 

&0 1 1 s 

3 2* 


?, (RATES( 1201 1 .RV I 

GNAT '041 



□ oil* 

3»* 


DOUBLE PRECISION PSAR, RBAR i REBA« 


CNaT0Q12 

00117 

31* 


double PRECISION hh, rr 


SN*T0013 

00120 

35: 


EQUIVALENCE UVEtll.PBARl. 1 AVE(3 1 iRBaR I 1 (aVEISI , REBAR) 


SNaTOOII 

00121 

36* 


dimension labeliioi 


SNaT0015 

□ 0122 

37* 


equivalence 


6NAT0D16 

00122 

36* 


1 {LABEL t 1 ) * T 1 


CNATC017 

00122 

3T* 


2 a (LABEL (21 > OT 1 ) 


SNaTOOIS 

00122 

10* 


3, (LABEL131 iTSCALEI 


GNaTOOH 

00123 

11* 


DIMENSION LIMITS150I 


S N*T0050 

00121 

12* 


EBU [ VALENCE lUHlTSt ll.NGl . (LIMITS 121 1 ,NSJ 


5 NaT0051 

00121 

13* 


l 1 (LIHITSlIl 1 iNOI . (LIMITS (12) ,NFl . ( L 1 M I TS 1 13 J t NODE S 1 


6 NaTD0S2 

00125 

ii* 


DIMENSION PROP ( 1 0 ) i CTL ( 20 1 i NG(2Q), NS(20l> 0BBC(1| 


GNAT0053 

00126 

I5t 


EQUIVALENCE (pROPtil.iTHli (PROP ( s 1 • 7. 1 i (PR0P(61,Ll* 


SNAT0Q51 

00126 

16* 


l t«R0PI7>i RC0NS7 1 • 1 PROP 1 1 1 (HO 1 


SNAT0055 

00127 

17* 


EQUIVALENCE IcTL 1 1 1 tin J i (cfL<2> ,0T1 , (CTL131 .TSTOPl , 


sn^toos* 

00127 

1#* 


1 ( CTL < 1 1 tDTPR I , (CTfc (20) iSToP ps > 


SNATD0S7 

00127 

If* 


2 , ( CTL ( 1 3 1 f SCALE 1 


GNAT0QS8 

00130 

58* 


DIMENSION I T APE (61 


GNaTDOSI 

00131 

51* 


EQUIVALENCE (CTLllll.iTAPEUII 


GNaTOOAO 

00131 

521 


1» (ITAPEU) iJUNJTR) j 1 JTaPE(2) , IFJLER1 i IJTaPE(3),IReCRJ 


GNaTODAJ 

001 M 

si* 


2 t ( | TAPE 111 1 1 UNI TW1 > (ITAPE(S)'IFILEW) r 1 1 TAPE ( 6) , | «EC«1 


GNa T0Q62 

00132 

51* 


EQUIVALENCE (PROPdO) iPFlAS) 


GNAT0D63 

00132 

55* 

c 



GNAT0041 

00131 

54* 


data PROP / 32.1 *2d0 . .08, % . 2E*3 » .2, .1, 1*0. / 


SNaTODAS 

00135 

57* 


DATA CTL / O.i .001. 60,, .01. O.i O.t 0.. 0, , 

0 . > 0. , 

GNAT0066 

00135 

581 


1 7*0. / 


GNAT0067 

00137 

5?« 


DATA NG / Bi 6 i 1> 3. 3, 2i 2. I. ti |. |> ]. Ii 2. 2i 3. 

3, 1, 

GNATC06* 

00137 

60* 


1 6, 8 / 


CMAT006* 

00111 

6 1 * 


DATA NS /l3iIS.l7,ie,ja,)V l |V,2o.20.20t2a.20>20<lf . 1 f . 1 S . 1 3 » 1 7 , 

GNATQ070 

00111 

62* 


1 15.13 / 


GNAT007J 

00113 

63* 


DATA ENDECK. is. LINES / 6H SEND , 10 ,S? / 


SNAT0072 

00117 

61* 


DATA PFLAG /O./ 


GNAT007 3 

00117 

*«• 

c 



GNAT0071 

00117 

66* 

c 

STATEMENT EUNCTlSNa 


GNAT007S 

00151 

671 


fh ( I , j) • R( i ■ j) * Rear 


GNAT0076 

ODZSI 

(*• 

c 



GNAT0D77 

00(52 

4 f 1 


NAMELIST /INPUT/ PROP, CTL, NS, NS ; iTAPE 


GNATQ07® 

00152 

70* 

c 



GNAT0071 

00152 

71* 

c 



GNaTODID 

0O|52 

72* 

c 

To input state from tape enter iunitr «, itile *>* i^c* * 

itapeh 

GN a TD081 

00(52 

73* 

c 



GNAT0082 

00 1 S3 

71* 


IUNITR ■ 0 


GNATD0B3 

00151 

751 


100 CONTINUE 


GKaTOOBI 

00155 

7»* 


CALL CR2TAP ( ENDEC K » 15, LINES 1 


GNaTOOBS 

00 156 

77* 


READ (15.INPUTI 


gnaTOObf 

0016) 

7»* 


IF ( IUNITR.LE.O) CO TO JQ5 


GNAT0007 

00163 

7?« 


CALL HDTAPE 


GNATOOBS 

00161 

ao* 


CALL CR2TAP ( ENOECK, IS, LINES | 


GHATOOif 

00165 

SI* 


READ (15, INPUT) 


GNaTOOI 0 

00165 

821 

c 



GNATDOf 1 

00170 

S3 1 


105 CONTINUE 


GNaTDQIZ 

00170 

81* 

c 

POSITION OUTPUT TAPE IF REQUIRED 


SNA T0093 

00171 

85* 


IF ! IUNITR. LE. 01 CO TO 180 


GHATD0T1 

00173 

86* 


kALL TAPER* (O.lUNim 


GNATOOfS 

00171 

S7i 


CALL TAPEPS (0, IUNITR, |F:lE*-| , IRECA-1 1 


GNaT00*6 

00175 

as* 


CALL TAPEPR lO.IUNITp] 


GNAT0017 

00176 

Bf * 


1 8 p CONTINUE 


GNATOOfS 

00177 

»0* 


OTI ■ OT 


GNaTDOT? 




■/a 


0017? 

41* 

C 


GNaTDIOO 

OD2DO 

92* 


NO * 20 

GNaTOIOI 

0020, 

9J* 


NF ■ 0 

GNAT0102 

00202 

94t 


NODES * Ot 

GNaTOIOI 

00203 

95* 


00 |25 1*1,20 

GNAT0104 

00206 

9B» 


01 * NGU1 

GNaTO 105 

00207 

97* 


02 * NSl J 1 

GNAT01Q6 

002|0 

9Bt 


IF t 01 • LE. 0 I GO TO 125 

GNaTO 107 

002*2 

99* 


NO * NINO '.NO, Jl i 

GNaTOIOB 

002)1 

109* 


NF * MAXO (NF.021 

GNaTO 1 09 

002 1H 

101* 


00 120 0*01 >J2 

GNAT01 |0 

00217 

102* 


120 NODES - NODES • 1, 

GNATO lit 

00221 

103* 


(25 CONTINUE 

GNaTOI 12 

00221 

104* 

C 

SET CONSTANTS FOR Ol-SSICN IN DIFFERENCE routines 

GNaTOI 13 

00223 

105* 


L3 * t**3 

gnatohh 

00224 

J04* 


CSTIL * 1* / L 

GNATOI 15 

00225 

107* 


CST2U ■ .5 7 L 

GNATOI 16 

03224 

toe* 


CST4L * ,25 / L 

GNATOI 17 

00227 

1D9» 


CSTIL2 * 1, / L**2 

gnatoi is 

00230 

110* 


CST2L2 * .5 / L**2 

GNATOI 19 

00231 

114* 


CST4L2 ■ .25 / L**2 

GNATOI 20 

00231 

112* 

c 

CONVERT units 

GNAT012I 

00232 

1 13* 


PO * PROP ( 8 J • 144. 

GNaTOI 22 

00233 

1 1 4e 


PHI * SOO. • 144. 

GNaTOI 23 

00234 

US* 


PUD « 840, • 144. 

GNATOI 24 

00235 

11** 


C ■ PROP 12 1 « 778.194 • 32.2 

GNaTO 1 25 

00234 

1171 


K « PROP 1 3 1 • 778,154 / 3600, • SCALE 

GNaTO 1 26 

00237 

1 1*'> 


VSC * PR0P141 * 2,3886*3 • SCALE 

GNaT0|27 

00240 

1191 


■DOT » CTL t S > ; 32.2 t 3400, « .075 * SCALE 

GNATO 128 

00241 

120* 


04HEAT 1 CtLUI • 778.156 73400* / 1 7 * 7 1.* SCALE 

GNATOI 29 

00241 

nil 

c 

units vsc . L8F B SEC/fT-2 Input as poise 

GNAT 31 33 

00242 

122* 


HCONST * 32*2 • 1545, / »TH 

GNaTD 1 3 1 

00243 

!23« 


OX ■ 32,2*CTL( 1 1 1 • SCALE **3 

GNaTO 1 32 

J02H4 

124* 


ST * 32.2‘CTLIU 1 • SCALE •*2 

GNaTO 133 

□ 0245 

125* 


DO 110 1-1.4 

GNATO 134 

00250 

124* 


DQSC(i) - CTL (1*41 • T7B.156 / 0400. • SCALE 

GNaTO 1 35 

00251 

127* 


t 10 CONTINUE 

GNaTO 1 36 

00253 

1 20* 


DQ1 - OQBCll) 

GNaTO 1 37 

00254 

129* 


002 *-Da3C(2l 

GNATOI 38 

00255 

130* 


003 - DQBC (3 1 

GNATO 139 

00254 

131* 


004 — 0QBC14) 

GNaT0|40 

00254 

132* 

c 

set initial conditions 

GNaTOIHI 

00257 

1331 


IF 1 IUNIT3.GT.0) GO TO 17D 

GNATOI 4 2 

00241 

1349 


T * TO 

GNaTO 143 

00242 

1 35* 


PBAR « PO 

GNATD | 44 

00263 

134* 


EO - BETA t TO 1 • 778. ISA * 32,2 

GNATO | 45 

00243 

117* 

ct 

RO - PO 7 U* NC0N5T* HOI 

GNaTO I 44 

00244 

138* 


PJ - PO 7 144. 

GNaTO | 47 

00265 

139* 


RO - OPTD' Pl.HO 1 732.2 

GNATO 1 4 B 

00246 

140* 


RBAR ■ RO 

GNATOI H 9 

00247 

14J9 


REBAR - RO • £0 

GNaTOISO 

00270 

142* 


REFG - FLOAT 1 NF- -U 1 77. *1. 

GN A T015I 

□ D271 

143* 


DO ISO 1-NC.fE 

GNATOI 52 

00274 

144* 


J! * NG(I) 

GNATO 1 53 

00275 

l 45 1 


J2 ■ N5[ n 

GNaTO 1 54 

00274 

1449 


DU ISO J-JltJ2 

GNaTO 1 55 

0030 1 

147* 


R I - 1 

GNATO J 56 

00302 

148* 


RJ - J 

GNaTO 1 57 



003Q3 14** PU.JI * - RO • l G**L* (REfG-RI J + G T *L* (REFG-RJl I 

0G303 I50» C H 1 1 i J ] * IPtl.J}* PBARI / ( I* RC oN sT • HQ I * RBAr 

0030“! |5|« PJ ■ | P(|,J)»PBAR I / |44. 

00305 J 52* RII.JI ■ QPTDt PiiHO J / 32,2 - RO 

00306 153* RUll.JI ■ 0. 

00307 154* R V 1 1 ,0) ■ 0. 

00310 ISS* Ft E ( | , J ) - E0 * R<||J) 

003)1 156* HI t ,J) * HO 

003)1 157* C 

00312 159* RHO ■ R ( { , J ) * Rg*R 

00313 ISS* Plltjl ” HTPRESt RHO, H<HJ* I * P8AR 

00313 | AO * C 

□ 03|4 IMr U( I , J | m 0* 

003 |S 1*2* ¥ t 1 * Jl ■ 0 • 

00316 *i,.» 1 So CONTINUE 

00316 164* C 

00321 165* | 70 CONTINUE 

00322 1661 TPR ■ T 

00323 |6?» ITERH » 0 

00324 1666 50 TO 49Q 

00324 l*?» C SKIP TO 490 * OUTPUT INITIAL CONDITIONS 

00324 1706 C 

00325 17|« 200 CONTINUE 

00325 172* C 

00324 173* PFLAGl • PFLAG 

00327 379* PREF ■ PBAR 

00330 175* IF [ PREF .ST. PHI j PFLAG • Q. 

00332 176* IF ( PREF ,LT, PLO | PFLAG ■ 1, 

00334 *D I AONOST tc* THE T£*T FOR eOUALITT BETWEEN NON-|NTeSE«S N*T NOT BE HEaNI 
00334 |77* IF I PlkAS .He. PFL AGJ 9 RRlTE U.VIOJ T. PFLAG 

Co 003SI 1786 *10 FORHAT UKO, E12.6, 52, 7HPF LAG *, F3.0 f ) 

<JT 00341 1 7 * * C COHPUTe 1ST DIFFERENCES for VELOCITY CALCULATIONS 

00342 ISO* CELL HEANA 

00343 1 a I • 00 270 I *N0 iNF 

00346 182* Jl * N£( 1 I 

00347 1 S3 * J2 * NS 1 1 1 

903S0 IBS* 00 270 J*Jl,J2 

00353 IBS* CALL OIFF ( t >J 1 

00353 196* C COHPUTE RATE OF CHANGE 0£ VELOCITIES 

00354 1 97* ORU(I.J) * - QRUU0X - QRUVoT . DPuX * GX* FRll.J) * 

00354 |9C* 1 VSC ■ (D2UDX2 + D2UDY2 ♦ 1 0ZUD*2*i)2V0XT t /3> I 

□0355 IS** DRV I I . Jl ■ • DRVUOX . DRVVDT . DPDY * CT* FRll.J) * 

00355 |*0* 1 VSC • [D2VDX2*D2¥DTZ * t D2VDY2*D2U0ST ) /3, | 

0035* 1*1* 270 CONTINUE 

0q36l 1*3* IF [ RDOT.GT. 0.) CALL BCOUTU.SOOT) 

03361 1*3* C UPDATE VELOCITIES 

QD363 1*4* DO 300 I»N0iNF 

00366 1951 J) * NCtll 

00367 1949 JZ * NS U ) 

00370 1*7* DO 300 J-.Jl.J2 

00373 1*9) RUlt.J] « RUI|,J| * oT * 0RU|],J1 

00374 |*f* RVlIjJi ■ RVII.J) » OT * DRV t I , . I 

□0375 200* Ull.Jl * RU { 1 , J) / FR(t,Jl 

0Q376 201* VII, JJ * RV 1 1 1 J 1 / FRll.J) 

00377 202* 300 CONTINUE 

-HP377 203* c COHPUTE DIFFERENCES FOR HEAT RATE EGUAT ION 

00402 204* CALL MEANS 

0C4Q3 205* DO 350 t-Np,NF 


GNftTOlSS 
GltATOIS* 
GNaTD|6Q 
GNaTD I 6 1 
SNATDI6Z 
GNATOi 63 
GNAT0164 
GNaTO) 65 
GNaTO l 66 
aN ft T0>67 

5NAT0168 

GNATO | 6 * 
GNAT0170 
GNaTO k 7 1 
GNATO 1 72 
GNATO 1 73 
GNATO 1 74 
GNaTO 1 75 
GNaTO 1 76 
GNATOI 77 
GNATO 179 
GNAT017? 
GNAT0190 
GNaTOSSI 
GNaTO 1 92 
GNATOI 93 
GNATO 194 
GNaTOIIS 

» 

&N ATO 1 96 
GNaTDI 97 
GNaTOISI 
GNaTOIB? 
GNaTOITO 
GNATOI* I 
GNAT 01*2 
GNAT0193 
GNaTO 1*4 
GNaT01»5 
GNaTOI** 
GNATO 1*7 
GNATO) f 9 
GNATOI** 
G n AT0200 
GNAT0201 
GNaT0202 
GNAT0203 
GNAT0204 
GNAT0205 
GNAT0206 
GNaT 0207 
GNAT02D9 

GNaTOZO* 
GNaT 02 1 0 
GNaTOZI 1 
GNATQ212 
GNaT02|3 
GN aTOZ 1 4 



00406 

20»* 



J1 * NG | 1 ) 

Q0407 

'•5Jt 



J2 • NS{ 1 ) 

DoS 10 

209* 



DO 350 J»J| .J2 

DOS 1 3 

209* 



CALL 0 l FF2 1 1 ,J1 

00413 

"MU 

c 

COMPUTE OENSITT AND HEAT RATES 

00414 

2 I l * 



ORIliJl * - DRUDX * ORVD* - RbAP • IDODX*DVOYI 

004|S 

2 1 4 * 



DRE(|iJ) * - DREUD X - 0REVD7 - IPBAR * R EBArMIdUOX 

OOllS 

2|3* 


i K • JD2HD*2*02H0T2l 

304)4 

Z 1 4 C 


35q 

continue 

00421 

ZlS* 



IF I ftOOT.GT* 0.1 CALL SCOUT 1 2 • WdDT 1 

00423 

21*. 



OOH » DOHEAT • FPL AG 

00424 

2| 7 * 



IF ( OOH . GT *0 • 1 CA l l HEATER ( 0SH.L3 1 

00424 

Z IS* 

c 

update HEAT AND MASS 

00424 

2l?t 



DO 450 J-AO.NE 

00431 

ZZO* 



Jl * NG t l 1 

00432 

22 1 ft 



JZ ■ NS [ I) 

00433 

222* 



DO 450 J.J1,J2 

00434 

Z23* 



nil .Jt u ri 1 1 ji ♦ or* oR( t . ji 

00437 

224* 



RE 1 1 ■ J 1 ■ RElI.J) + DT * DREtt.J) 

90440 

225* 



RS * R I I ,J) * RBAR 

00441 

22*. 



UII.J) * RUU.JJ / RH 

0D442 

227* 



VUeJ) « RVII.JI / RR 

00443 

228* 



ERMO » REll.Jl * REBAR 

00444 

224 ft 



RHO » RR 

0Q44S 

230* 



HCl.J) * TEMP) ERMO, SHO ) 

0044* 

23)* 



HH * Mil ,J) 

0044* 

232* 

c 


CALL PRESS 1 PU.JI , R||,J| , HH 1 

00447 

203* 



PtJ.Jl ■ RTPRE5IRH0.HF J ,J1 I - PbaR 

00447 

234ft 

c 



0Q4S0 

235ft 



IF I RHO , 0 .) ITERM * l 

00452 

23** 


450 

CONTINUE 

00455 

237» 



CALL 4ULK 

0045* 

238* 



IF 1 1 VERM .EG. 1) GO TO 480 

004*0 

239 ft 



IF 1 *LT* TpR-.OQQOll GO TO 500 

004*2 

240* 


480 

CONTjNUE 

004*3 

24 l ft 



TSCALE ft T * SCALE / *ii. 

004*4 

242* 



CALL OUTPUT I T. DT, TSCAtE. 0,, HO 1 

004*5 

243* 



[F [ ITERM ,E9, | ) 60 TO 550 

004*7 

744* 



TPR ■* TPR ♦ OTPR 

00470 

245* 


500 

CONTINUE 

03471 

244* 



T » T ♦ OT 

0057? 

247ft 



IF 1 T.LE. TSTDP*DT*l,E-5 l Go TO 20Q 

0D474 

24fc ■ 


5S0 

CONTINUE 

00475 

24\ ■ 



IF I 1UNITN.GT.0I CALL TAPEOF IQ.lUNITRl 

00477 

25Lft 



IF <IUNtTW.GT.01 CALL TAPER* IQt(UNITM) 

00501 

251* 



IF t STOPFG <LT. 1* 1 GO TO 100 

00503 

252ft 



STOP 

00504 

258 * 



END 


END OF UNtY*C 1108 FORTRAN V COMPILATION. 1 'DIAGNOSTIC* HE5SAGEIS1 

MAIN STH&QLlC 24 DEC 70 Q9l|7! 

MAIN C<jO£ RELOCATABLE 23 DEC 7p 23I04I 


GNAT0215 
GN A T 0216 
GNaTOZ 1 7 
GNAT021B 
GNaT0219 
GNjiTOZZO 
GNIVT0221 

GNaTQ222 

GN, v T0ZZ3 

GNftT0224 

GNAT022S 

GNAT0226 

GNAT0227 

GNAT0228 
GNATD229 
5NFT0230 
GNAT0Z3I 
GNaT0232 
GNAT0233 
GNAT 0234 
GNAT0235 
GN*T0236 
GNATD237 
GNaT 0238 
GN*T023? 
fiN*T0Z40 
GNAT024 1 
SNaT0242 
GNAT0Z43 
SNAT0244 
GNaTDZRS 
GNaT024* 
GN*TD247 
GNAT0248 
GNAT02“9 
GNAT02SD 
GNAT02S 1 
GNAT0ZS2 
GNAT0253 
GNATOZSH 
ENAT0Z55 
GNaTOZS* 
GNaT 02S7 
SNA tGZSB 
GNaTOZS* 
GHaTOZAO 
GMATOZM 
G»aT02*2 


18 

0 

01*32134 

14 

304 

•DELETED] 

1* 

1 

01*5054* 

3* 

l 

IDtLETEBl 


0 

01*50*32 

14 

70 




B FOR HE AN A , mEaNA 

Off I U AC HOj FORTRAN V LEVEL 2204 00IB F501BH 
This compj .aTion mas done on zt jan 7i at tmoise 


SNAT02A 


27 jan 7 i 


)1i33:3B.13S 


SUBROUTINE NEANA 

ENTRY POINT 0Q03S*) 

STORAGE USED 

BLOCK, 

NAHE, LENGTH) 

00(1 

•CODE 

000373 

0000 

• DATA 

000017 

0002 

•BLANK 

000000 

0003 

data 

007 I 41 

0001 

nooaal 

010172 


EXTERNAL REFERENCES (SLOCK, NAME) 
0005 NE&R3S 


STORAGE 

assignment for 

variables 

(BLOCK 

. ttpe. 

RELATIVE 

LOCATION 

. NAME) 









0001 


ooooos 

1 I 3G 

000! 


0000*1 

12HG 

0001 


DOOMS 

20b 

□ 001 


□DO 

30L 

0001 


D0020* 

sol 

000 T 


000257 

70L 

□ 001 


0D03Q0 

bol 

0001 


□00320 

*0L 

0002 

R 

000000 

A 

DODO 

l 

000003 

1 

0000 

1 

oooool 

II 

0000 

1 

000002 

12 

0000 

I 

000000 

J 

0003 


000050 

limits 

0003 

I 

000 1 2 1 

NF 

0003 

I 

oooqso 

NG 

0003 

1 

00007*) 

NS 

0001 


000120 

NO 

0003 


0010*1 

RATES 

ODOR 

R 

0047 1 0 

RE* 

ODOR 

R 

0074(J| 

Ret 

0002 

R 

005521 

RU 

DOOR 

R 

000000 

RUX 

0001 

R 

00047 | 

RUT 

0003 

R 

004311 

RV 

0001 

R 

00)542 

RVX 

0005 

R 

0Q21S3 

RVT 

DOOR 

R 

□05124 

RX 

DOOR 

R 

00*0*7 

RT 

0003 


COO 1 IN 

STA 

0009 

H 

002121 

U 

0001 

R 

003311 

UX 

0003 

R 

0032*1*) 

V 

CSQ1 

R 

001235 

VT 






00101 

1 • 


SUBROUTINE HEAnA 

GNaTOZAI 

00)01 

2* 

C 


GNAT024S 

00)01 

3* 

c 

Interpolate values of velocity and hqhentuh at nooe faces 

GNaT02** 

ODIJt 

1* 

c 


GNAT0247 

OOtCi 

6» 


CQHHON /DATA/ A137001 

C-NAT024E 

00101 

,♦ 


COMMON /NODRAL/ RUX 1 2 1 ,2 i , RUT (2 1 i 2 1 ) ,RVX 1 2 1 .2 1 1 , RVT (2 1 . 2 1 ) 

6NAT024? 

00(01 

Jt 


1, UX12I *2| 1 ,yY (2! ,2U ,RX (21 ,21 1 ,RT (2| ,21 ) ,REX 121 ,2) ) ,REY ( 2 1 ,2) } 

GNAT0270 

G0I01 

a* 

c 


GHATQ27I 

00105 

1* 


DIMENSION RATES 1 1400) .STATE ( 2003 1 .LI HITS ISO) .NS I 20) «nS{20) 

5NAT0272 

0OIC4 

10* 


01 HENS I ON U(20,Z0 1 , V < *0 . 20 ) *RU ( 20 * 20 1 ,R V (20 . 20 ) 

SNA : 027* 

□ 0107 

11* 


EQUIVALENCE (Aill ) .LIPITS! , I A ! 1 0 1 ) , STATE 1 , (A (210! I .RATES ) 

GNAI 0271 

001 10 

■21 


EQUIVALENCE (STATE ( ) 20 1 ) ,U 1 • I5TATEC140U.V1 

GNaT 0275 

00110 

I3« 


!, (RATES (SOI 1 ,RU) , (RATE5 ( 1 201 > iRVl 

GNAT 027 4 

0011 1 

lit 


EQUIVALENCE ( L IH I TS I l * ,NG I , (LIMITS (21) »N5 | 

GNAT0277 

□ 0111 

15» 


1, l L 1 H ITS ( 1 i » ;0 ) » (LIMITS (12 )*NF) 

GNAT0Z7B 

00111 

14* 

c 


GN AT027V 

00112 

1 ' * 


00 100 »*N0.<1F 

GNAT02B0 

□OHS 

]8« 


11 * NG ( J ) * 1 

GNAT028I 

00114 

IV* 


12 ■ NS ( J ) 

GNATOZBZ 

-nan? 

20c 


UX(U-1 ,J) * 0, 

GNATD283 

00120 

21* 


UX ( (2*1 , J) * 0. 

GNATD2S1 

0012! 

22* 


VTiJ.lMI? * 0. 

GNaT02BS 



00)22 

23 r 


VY(J.I2*1 1 

■ 

0. 



00123 

24* 


DO 100 1*11 

* 

12 



00126 

25* 


OXt 1 > J 1 * 

1 

Uli-1 . j) * Ut 1 ,J) 1/ 2. 



00127 

264 


VTIJ.I I « 

1 

V(J.I-l) * V(J,|I 1/ 2, 



00130 

27* 


1F( IUXII,J1 

,(,T,0.) .AND. < I .GT. : U 1 GO TO JO 


D0132 

2B* 


I F t lUXII.J) 

.UT.G.I. A ND.H.LT.12) 1 GO TO 30 


00131 

29* 


RUXIJ.JI • 

1 

HUt 1-1 t J) * RU( I , J) 1/2. 



00135 

30* 


RVMliJl * 

( 

RV 1 1-1 >Jl * RV 1 I .Ji i/2* 



00136 

31* 


GO TO 50 





00137 

32* 

20 

CONTINUE 





aomo 

33* 


RUX ( 1 i J) * 

t 

-*RU 1 | *2 . J 1 *6»*RU( 1-1 . Ji 

*3 »*RU 1 1 . J | 

1/8. 

00191 

34* 


RVXII.ji * 

( 

*RVU-2iJi »6.*RVl 1-1 .Ji 

*3 • *RV 1 1 t JI 

J/B. 

QDI 42 

35* 


GO TO 50 





00143 

36* 

30 

CONTINUE 





0DI49 

37* 


RUXIliJ) * 

( 

3.*RUt I-i 1 J> *6,*RUU.J) 

-RUt j*l ,J1 

1/8. 

00195 

36* 


RVXll'J* * 

( 

3**RV<l-l.Ji **s*RVll<Ji 

— R V 1 | * 1 1 J 1 

J/B, 

CO 1 46 

39* 

SO 

continue 





00147 

40* 


IF t lVT(J, 

1 | ,GT.0,1.AN0, ( 1 ,GT,I1 1 1 

GO TO 70 


00151 

41* 


IF 1 (VTIJ. 

n.U*0H.AN0.U*LT?!21 1 

GO TO BO 


00153 

42* 


Rurijii) ■ 


RUfJ.I'll * 9U» J. I 1 i/2* 



00154 

434 


RVTtjil) - 


RVl J. l-l ) * RV t J. I ) 1/2. 



00155 

44* 


GO TO 90 





00156 

45* 

70 

CONTINUE 





00157 

46* 


RUT 1 J * 1 1 * 


«RU I J . 1 *2 } *6.*RutJ.]-H 

*3 • *RU 1 J • 1 1 

1/8. 

00)60 

474 


RVT 1 Ji | t " 


»RV*Jil-2i +6 .*Rv 1 Ji l-l ) 

*3 • *RV 1 J 1 1 1 

)/8. 

00161 

46* 


GO TO 90 





00162 

404 

80 

CONTINUE 





00163 

504 


RUY|J,|) * 


3 • *RU 1 J . 1-1 1 *6,.rUIJ.M 

• RU(J> 1*1 1 

1/8. 

00164 

51* 


RV Y ( J » 1 1 * 


3.*RVtJ.I-li 46..RVtJili 

T RV 1 J 1 1 * 1 1 

1/8. 

00165 

524 

90 

continue 





00166 

53. 

|00 

continue 





00171 

54* 


return 





00172 

55* 


END 






END or UNIVAc hob FORTRAN V COMPILATION 


0 'DIAGNOSTIC* HFSS AqE ( 5 ) 


GNAT0286 
GNAT0287 
GNAT02B8 
GN A T0289 
GN A T029D 
GN a TC 29 l 
GN A T0292 
GNAT0293 
GN A T02»1 
GN A T02« 
GN A T0246 
GNATD297 
GNXT029B 
GN A T0299 

gn a tqjqq 

GNAT030 | 
GMT03D2 
GN A TD303 
GN a T 0304 
GNATD30S 
GNAT030* 
GNAT03D7 
GNHT030® 
GN A T03C9 
GNAT03I0 
GNAT03H 
GNUEUI2 
GN A T03|3 

GNATOil 4 

GNAT03IS 
GNAT0316 
GN A T03| 7 
GN A T03|B 



B FOR MEANS , MEANS &NAT031* 27 JaN 7) l *i 1 33 ; NO,*! ?N 

univac mos fortran v level 2204 oojb fsoish 

This COmPILATJON fiAS done ON 27 JAN 7] AT IR133JR0 


SUaROUTlNE KEANS 

ENTRT POINT 000332 

STORAGE USED 

[BLOCK, 

NAHE, LENGTH} 

0 00 1 

•CODE 

0O03N7 

aooo 

• DATA 

OODON3 

0002 

•BLANK 

000000 

0003 

DATA 

007 l AN 

DOOM 

nodral 

010N72 


EXTERNAL REFERENCES [BLOCK, NAHD 
ODDS NERR3J 


storage 

assignment for 

variables 

(BLOCK 

. TTPE, 

RELATIVE 

LOCATION 

t NAHtl 









0001 


000005 

II3g 

0001 


oooom 

1205 

0001 


0001 25 

20U 

ooot 


0001R* 

3Ql 

ooot 


000144 

SOL 

L-OI 


000237 

70L 

□ 001 


□002*0 

«0L 

000 I 


000300 

»0L 

0003 

fi 

000000 

A 

0000 

I 

000003 

I 

0000 

t 

000001 

It 

0000 

I 

000002 

12 

0000 

I 

000000 

3 

0003 


000050 

LIMITS 

0003 

1 

000121 

NF 

0003 

I 

000050 

ng 

0003 

l 

00007R 

NS 

0003 

I 

oaouo 

NO 

0003 

R 

00074R 

R 

0003 


00 04N 

RaTi 

0003 

R 

00R70R 

RE 

OOOH 

fi 

00471a 

REX 

DOOR 

R 

007401 

RET 

oaov 

R 

000300 

RUX 

DOOR 

R 

000*7 l 

RUT 

ODOR 

R 

001542 

RyX 

OOOH 

R 

□ 02 V S3 

RVY 

ooov 

ft 

005124 

RX 

ODOR 

R 

00*017 

RY 

0003 


OOOINV 

sta 

0003 

R 

002H2R 

U 

odor 

R 

Q032RN 

UX 

□ 003 

p 

0032NR 

V 

ODOR 

R 

00RZ35 

VT 






oa id 

u 


SUBROUTINE MEANS 

GNATD320 

00 10 1 

2» 

C 


GNAT0321 

00101 

3* 

C 

INTERPOLATE VALUES OF VEbOCITT , DENslTT, AND ENERGT AT N<>0E FACES 

GNAT0322 

00101 

R * 

C 


GNAT0323 

□ 0103 

S» 


COMMON /DATA/ A 1 3700 ) 

GNAT032R 

00 IOR 

*• 


COMMON /NODWAu/ RUX(2(,2i},RUT(2U2n,RVX(21r2|},RVT(2t,21) 

GNAT03Z5 

DOIOR 

7 • 


1 1 UXI21 ,21 l.vr 121 1 2 ) 1 (fix 1 21 .211 ,RV*2l ,21 1 .HEX (2 1 ■ 2 1 J • RET < 2 1 , 2 1 1 

GNAT0324 

OOIOR 

St 

L 


GNAT0327 

00105 

V* 


01 HENS I ON RATES (1*00) .STATE <20001 l L I M ITS (5n ) ,NG ( 20 ) , NS [ 20 ) 

GS A T032B 

0010* 

10* 


DIMENSION U 120.201 , V 120, 20), R 120,20) ,RE( 20, 201 

SNAT03/V 

00107 

11* 


ESUIVALENCE (AIR)), LIMITS), [ A (1 0 1) , STATE ) , 1 A (2 10 1 1 ,RA TE5 1 

GNATD33Q 

ooi to 

121 


EfiUlvALENCE (ST ATE 1 | 2D 1 ) ,U ) , [STATE 1 1*01 1 .VI 

GNAT033I 

001 10 

13* 


l, (STATEIROl ) ,RI , 1 RATES (RO II ,RE ) 

GNAT0332 

00111 

1R* 


eouivalence (biniTsiiiiNGi , (limitsizjmns; 

GNAT 0333 

0DU1 

IS* 


I. CLIMITSIRU .NO) , ILIHJTS(R2) ,NF! 

GNA T033R 

00111 

1 *• 

C 


GNaT 03 3 5 

00112 

1 7 v 


DO 100 J-NO.NF 

GNATD334 

00115 

IS* 


11 * NGtJI * 1 

GN*T0337 

□ 01 1* 

1»» 


12 - NSIJ) 

GNATD33* 

'TTUj 17 

20* 


DO 100 1-11)12 

GNAT033V 

00)22 

21* 


UX(I,U) ■ { * U|I,JJ )/2, 

GNAT03R0 

00123 

22 • 


VYIJ.il - ( VIJ.I-1) • HJ.l) I/P. 

GNATD3R1 



00121 

23* 


IF 1 1UXI1 

-J>. 

GT.O.I.AND.l J.GT.11I 1 GO TO 20 


00|2* 

29* 


I F I 1UX<) 

. J I • 

LT.Otl.ANO.ll.LT. )2) 1 GO TO 30 


0Q13U 

25* 


R XI1.J1 

* ( 

R ( I"! U) * R t l .J) 1/2. 



0013] 

2** 


REX 1 | i J) 

* t 

REU-t.J* ♦ HE)1,J) 1/2. 



00)32 

27<r 


GO TO 50 





00133 

2B. 

20 

CONTINUE 





00139 

291 


R X) I iJl 

« ( 

•R (1*2, Ji **.*R (I.|,Jl 

*3.*R 11. J) 

I/S. 

00135 

30* 


REX ( | i J 1 

- ( 

-RE)l-2,J) **<*RE(l-ltJl 

*3 ..RE1 1 ■ JI 

I/S. 

0013* 

3|» 


Go TO So 





00117 

32< 

30 

CONTINUE 





00190 

33. 


R XI), JI 

« ( 

3..R tl.I.Jl *A.*R II, JI 

“R (1*1, J) 

I/S. 

oont 

39* 


REX 1 | • J) 

- 1 

3,»REll*i,JI ♦* ■ *RE I 1 ■ J I 

-RE) 1*1 i J) 

1/8. 

00112 

35c 

50 

CONTINUE 





0Q113 

3** 


IF | (VT|J,|J 

,GT,Q<), AHD.lt, ST, I|l 1 

CO TO 70 


00115 

37* 


IF ( (V y (J,tl 

,LT,0. | .AND, l 1 ,LT. 121 I 

GO 70 SO 


00197 

3B. 


H Ttj,J> 

* S 

R <J»l*ll * R )J.I» 1/2. 



00150 

39* 


RETS Ji|l 


RE) Ji 1*1 1 4 RE 1 J . 1 1 1/2. 



00151 

90* 


GO TO 90 





00)52 

Ilf 

70 

CONTINUE 





00)53 

92* 


H TC J,) 1 

* 1 

-R (J.I-Z1 **,*R 

*3.*R (Jill 

1/8. 

00)59 

93* 


RET) Ji 1 1 

• t 

*RE ) J ■ 1*2 I »*.»RE)Ji 1*1 1 

*3 .*RE t J , 1 1 

I /Si 

00)55 

99* 


GO TO 90 





OOIS* 

95* 

BO 

CONTINUE 





00157 

9*« 


R Tlj t |l 

* ( 

3.*R IJ.l-l) t*..R ( J , 1 1 

* R (J.l+1) 

1/8. 

DO 1*0 

97* 


RETCj.)) 

- ( 

3.*REIJ.I*D **..RElJ.ll 

- RE 1 J » 1 *i f 

1/8. 

001*1 

IS* 

»a 

CONT )NUE 





oauz 

99* 

too 

CONTINUE 





001*5 

SQo 


RETURN 





001** 

SI* 


END 






END OF UN J VAC I1DB FORTRAN ¥ COMPILATION. 0 .DIAGNOSTIC* MESSAGEISJ 

X) 

o 


gnaT039Z 

GNAT0393 

GNAT0391 

GHAT039S 

GNATD39* 

GNAT0397 

GNaT039B 
GNATQ399 
SNAT03SD 
GNAT035I 
GNAT0352 
GNaTD353 
GNAT03S9 
GNAT035S 
gnaTOOS* 
GNaT03S7 
GNATQ3S8 
GNAT03S9 
GNaTOIaO 
GNaTOJAI 
GNAT03*2 
SNaT03*3 
GNaT0349 
GNfTOTAS 
GNA ’ •>* 

GNA ' 7 
GNAl- ..A 
GKaT03*9 
GNAT0370 



B FOB DjfFi OJFF 6NATD371 27 jAN 7| I4l33i 

WN1VAC HOB FORTRAN V LEVEL 2204 QO10 F301BH 
This compilation «as done on 27 Jan 7i at ]4)33;42 


SUBROUTINE PIFF 
STORAGE used (SLOCK, 

0001 ‘CODE 
□□00 «OAta 

0002 *8LaNK 

0003 DATA 

0004 U1 

0005 NOC*AL 
000b CSTS 


External references 

0007 NERR3* 


ENTRY POINT 000413 

NAME, LENGTH) 

00047M 

000032 

000000 

007)44 

OODG14 

0)0472 

000004 


(BLOCK, NAME) 


storage assignment for variables iblocki type, relative location, hahei 


0001 


000152 

7) ol 

ooo* 


000243 

2IIL 

0001 


000243 

zhl 

NODI 


OOQ3o 3 

22ql 

000 ! 


000374 

ZZ1 L 

000) 


□004 ) 4 

Z22L 

0001 


000434 

23QL 

0001 


□00473 

23 JL 

0001 


0005)5 

Z4Ql 

0001 


000554 

241 L 

0001 


000575 

3Q0l 

0003 

R 

000000 

A 

000b 

R 

000003 

CSTl* 

0004 

R 

000000 

CSTlL 

000b 

R 

000003 

CST1L2 

0004 

R 

000001 

CST2L 

0004 

R 

000004 

CST2 L 2 

000b 

R 

□00002 

CSTML 

□ 004 

R 

000005 

CST4u2 

□ 004 

R 

300000 

DPOX 

□ 004 

R 

000001 

OPPY 

0004 

R 

000010 

DRUUoF 

0004 

R 

000013 

PRUVOY 

DOOM 

R 

000012 

DRVUo* 

0004 

R 

□0001 1 

ORVydT 

0004 

R 

ooqoqa 

DZUOXY 

□ 004 

R 

000002 

D2UDX2 

0004 

R 

□00004 

02U0YP 

0004 

R 

oooao7 

D2V0XT 

0004 

R 

000003 

D2VoA2 

0004 

R 

000005 

P2V0Y2 

0003 


001404 

H 

0000 

I 

000001 

U 

0000 

I 

000002 

12 

□ODD 

I 

000003 

Jt 

ODOO 

I 

DOOOOM 

JZ 

ogoo 

R 

000000 

L 

0003 


QDOQSO 

LIMITS 

0003 

1 

000050 

NG 

0003 

1 

000074 

NS 

0003 

R 

000144 

P 

0003 


004044 

HATES 

0003 


0047Q4 

RE 

0005 

R 

004710 

REX 

0005 

R 

□07401 

Ret 

0003 


005524 

RU 

0005 

R 

000000 

RUX 

0*105 

R 

0QQb7l 

RUT 

0003 


00*344 

RV 

0005 

R 

0015*2 

RV* 

0005 

R 

002453 

RVY 

0005 

R 

00512b 

RX 

OCE’J 

R 

004017 

RY 

0003 


000)44 

STATE 

0003 

R 

002424 

U 

0005 

R 

003344 

UX 

0003 

R 

003244 

V 

□ 005 

R 

004235 

VY 










00101 

1 ■ 

□ 0103 

2* 

00104 

34 

00104 

44 

□ OIOS 

5* 

00105 

4* 

00104 

7* 

00107 

8* 

00110 

9 * 

001 11 

10 11 

001 11 

11* 

■ssouz 

1 2* 

00113 

1 3 * 

00U3 

14* 

00114 

15* 


SUBROUTINE DIEF III 71 
COHMON /DATA/ A(3700^ 

COMMON /Ol/ 0 PPXiPPDYiD2U0X2,d2VDX2)D2UPY2,D2VDT2i 
1 p2UDXY ,02VSXY iORJUDX iDRWoT iORVUOX ,DRUVOT 
COMMON /NOD*AL/ RUX ! 21 , 2 1 I ,RUY (2 1 .2 1 ) YX ( 2 J , 21 ) , RvY 1 2 I , 2 I I 
| , UX(2| ,21 I ,VYl2l i2l I iRX(2| i2l > RY(2) «2) ) ,REXI2| ,21 I > RE Y ( 7 ( , 2 1 I 

COMMON /CSTS/ CST IL , CSTZL iCST4L iCST 1 l 2 , CST2L2 1 CST4’. i 
EQUIVALENCE (CSTL2,CST1LZ1 

DIMENSION RATES ( 1 400 I ■ RE I 20 , 20 I , RUI20tZ0), RVI20.20) 
EQUIVALENCE (A(2!DI), RATES), (RATES ! 40 I I , RE I 
], (RATES (001 ) ,RU> , (RATES I 1201 I .RV) 

EQUIVALENCE I A(4U (LIMITS) ■ ( A I I 0 11 i ST ATEI 
DIMENSION LI HITS (SO) ,NG 120) .NS (20 ) tSTATE l 2000 ) 

1, P (20 » 20 ) • HlZO, 20) „U (20,20) iV(20 |20) 

EQUIVALENCE ( L I HI TS ( l I ,NG I , ( L I N I TS l 21 ) ,NS ) 


GNAT0372 
GNAT0373 
GNATC374 
GNA10375 
GNA r 0374 
GNAT0377 
GNAT037B 
GNAT037V 
GNAT03B0 
GNA T03B 1 
GNAT03B2 
DNaT03B3 
GNAT03B4 
GNAT03BS 
GNaYMSB* 


-/.*|S 



0DII4 

l 6 1 


1. (STATE! 11. PI 

00114 

179 


3. (5TATE1 80U.H) 

00114 

!»• 


4. (STATE ( 1201 ) ,U) 

-0114 

19* 


S, (STATE (16011 ,V1 

00115 

20» 


REAL L 

00116 

21* 


L * l./CSTIL 

00116 

22* 

c 


00117 

23* 


ll * I6(JJ 

00120 

24e 


12 - NS|j) 

00121 

25» 


31 * NGI I 1 

00122 

26» 


J2 - NS III 

001ZZ 

27* 

c 


□ 0123 

2B* 


DRUUPX * IRUX||*1 ,J1.UX( 1*1 ,J) - RUX " 1 , J 1 »UX l I . J 1 ) /L 

00124 

2V. 


DRVP-JX * (RVX(|.l,Jl«uX(l*i,Ji * HVX{I,J).tfX{J,JlJ /L 

00125 

30* 


ORli-OT * (ROT! 1 , J*1 l.yrl | , J»1 1 « ROY (]tJl*VT(IiJl! 7L 

00126 

31* 


DR WOT ■ (RYT(1 iJ»l>*VTtl .J*U - R9Y l J , J 1 *VT ( I , JV 1 7 l 

00126 

32* 

c 


00127 

33e 


IF II.EQ.n | GO TO 2 1 0 

00131 

34* 


IF ll.EQ.I2l GO VO 220 

001 J3 

35* 


IF 1 J.EQ.31 1 GO TO 230 

00135 

36. 


IF (J.EQ.J2) GO TO 240 

00135 

’H 

c 

COMPUTE AT CENTRAL POSITIONS 

00137 

36» 


OPDX * CST2L *1 Pt|*l.j) - Pll-I.JI 1 

00140 

39. 


DPOT * CST2L “< 8 1 I 1 J* 1 ) - pll.3-1) 1 

00141 

40* 


02UDX2- CSTL2 •( U(I.l.J) .2,.U(|.J> *0(1*1.31 1 

00142 

4 1 » 


0290X2* CSTL2 •( V«1*1.J) -2..y(l.jl *Vl|*t.J) 1 

00143 

42* 


0200Y2 » CSTLz *( yll.J-U -2.«U(I.JI »U(I,J*1) 1 

00144 

43. 


02V0Y2 ■ CSTL2 =1 VlI.J-ll -2..V(l,J) *9(1. J*l) ) 

03145 

44* 


D20DXY * CST4L2 * 1 U l l » 1 . J* 1 1 *U ( 1 - 1 . J- 1 ) -U ! I ♦ l . J* 1 1 -U U » 1 ♦ 1 ) 

00144 

4»» 


02VDXY * CST4 U Z • 1 9 1 I *1 . J* l 1*9 t I -1 . J-l >-V ! I *1 . J-l 1-V 1 J -1 . J*1 ) 

00147 

46* 


GO TO 300 

00 1 47 

47* 

c 


DO 147 

48 • 

c 

cohpute at left hall 

00150 

49. 


Z 1 0 cOMTlNvE 

00151 

SO* 


OPDX - CST1L *1 P 1 I *1 .3 1 - PII.JJ 1 

00152 

519 


D20Dx 2 " CSTL2 *t U(l*ltj) - 3.*U<liJi * 

00153 

52* 


02VDX2 * CSTL2 •! 911*1, J) - 3 , *911,31 1 

00154 

53* 


IF U.EQ.Jl) GO TO 2(1 

00156 

54* 


IF (J.EQ.J2) GO TO 2)2 

00160 

55* 


DPDY * CSTZL *1 Pl|.J*l) - ptl.J.ll ) 

00)6| 

56* 


D2UDT2 * CSTL2 ■! U t j * j- 1 » -2,*U(J,J) *Ut|,J»|) ) 

00162 

5?* 


D2VDT2 « CSTL2 *1 Vll.J-l) •2.*V(1.J> *V'1,J*11 ) 

00163 

58* 


02U0XT * CST2L2 • ( U 11 * l , *♦ 1 > *U ! 1 , J-l 1«U 1 1 ♦ 1 , J-l 1 -V{ J , J*1 1 ) 

00164 

59* 


DZ VDXT * CST2L2 • 1 V 1 1 * 1 • J*1 1 *9 1 1 ,J-I J-9 1 J *1 1 J-l 1 -9 1 i . J* 1 1 ) 

00 It 5 

608 


GO TO 300 

00166 

61* 


2 1 | CONTINUE 

00167 

62. 


02UDXY * CST1L2 *1 0 ( 1.1 , J. 1 J .0 U , J | -0 ( 1 , J t l )-0 ( 1 ♦ 1 , J 1 ) 

00170 

o3« 


02V0XT ■ CST1L2 *( V I J*l ,J*1 l*V( 1 • Jl-» 1 1 ,Jtl )-9l |M . Jl 1 

00171 

64* 


GO TO 23 | 

SJ172 

65* 


212 CONTINUE 

00173 

66* 


020DXT • CST1L2 •! U| 1*1 , J>*V( I , J-l l-WI 1*1 ,J-l l-0| 1 , J) ) 

00174 

67* 


0Z90XT ■ CST1L2 •! 9 1 l * l . J 1 *9 ! 1 . J* 1 1-8 l I *1 . J-l 1 -9 1 I . J I 1 

00175 

60* 


GO TO 24J 

00175 

69* 

c 


00175 

70. 

c 

compote at right wall 

00176 

„♦ 


220 CONTINUE 

00177 

72* 


OPOX • CST1L * 1 PII.jl - PU-t.JI 1 

00200 

78* 


□2UDX2 « CSTLZ * ! Jlj-i.J) - J.'UlI.Jl t 


SHAT0387 
GNATO3B0 
GN;,T0J89 
gnaTOJbg 
GSATQ 39 l 
GNAT0392 
GNATD393 
GNAT0394 
GNaT039S> 
<5H*T039* 
GNATQ397 
GHAT0398 
GNATD399 
GNaTOIOO 
GNaTOIO) 
GNATD4Q2 
GNAT0403 
GNaT0404 
GNaTOHQS 
GNaTDIO* 
CXATG407 
GNATOHCe 
GNAT0409 
t'lAtOHIO 
GNaTO**; | 
GNAT0412 
GNAT04] i 
GNaT0414 
GNaTONIS 

F (A T04 1 * 

GNATQ4I7 
GHaTOH IB 
SNaTQ 4 l Y 
GN7.T04Z0 
GNaT042I 
GNAT04Z2 
CNAT0423 
GNAT0424 
CNaT042Si 
GNAT04. 6 
SNAT0427 
GNATC4Z8 
GNAT0429 
SNaTOHJO 
GNAT0431 
GNaT 0432 
GNaTD433 
GNaT0434 
CNAT0435 
GNaT0436 
GNaT04 3 7 
GNaT043B 
GNAT093? 
GNATQ44Q 
SNa7&441 
GN*TD442 
GNaTDH43 
GNATQM 49 



VO 

OJ 


00201 

71* 



00202 

75* 



00201 

76* 



0020A 

77* 



00207 

76* 



002 1 0 

79* 



0021 1 

SO* 



00212 

81* 



002 1 3 

S2* 



00211 

S3* 


22 1 

002 IS 

SI- 



002) A 

85* 



00217 

BA* 



0Q22P 

87* 


222 

00221 

88* 



00222 

89* 



00223 

90* 



00223 

91* 

c 


00223 

92* 

C 

CO 

00221 

95* 


230 

0022$ 

91* 



0022A 

95* 



00227 

9** 



00230 

97* 



00231 

98t 



00232 

99* 


231 

0Q233 

100* 



00231 

lot* 



0023$ 

102* 



0023$ 

103* 

c 


0023$ 

101* 

c 

CO! 

O023A 

IDS* 


290 

00237 

lot* 



00210 

107* 



00211 

108* 



00212 

109* 



00213 

110* 



00211 

111* 


211 

00/1$ 

112* 



0021 A 

1 13* 



00216 

1H* 

c 


00217 

IIS* 


300 

QQ2S0 

l 1 A * 



00251 

l|7* 




D2VDX2 - C5TL2 =1 VlI-l.JJ - 3 r*V<)>J) » 
IF IJ.E«,J1J GO TO ZZ! 

IF IJ.EU.JZ) GO TO 222 


DPOT 
OZUOTZ 
□2VDY2 
J2UDXT * C5T2L2 
02V0XT ■ CST2L2 
GO TO 300 
CONTINUE 
D2U0XT * CSTIU2 
D2V0XT * CSTIl.2 
GO TO 23 J 
CONTINUE 
D2UTXT • CSTIL2 
O2V0XT * C5TIU2 
GO TO 211 


CST2L *1 e T I iJ*l 
CSTL2 •( UtliJ-1 
C5TL2 *1 YlliJ-t 


• I 
•I 


U( 1 , J 
¥<I.J* 


I * Pi | , J-l I I 

> «2.*U(1,J) *U 1 1 , J* I I I 

> -Z>*Vll»J) *¥U,J*1J ) 

i )*ui i-i , j-i i-uc i-i ,j*n-ui j i j-i i •, 

l*¥t l-I.J-ll-¥t l-l.J-ll-Vl I. J-l I ! 


uil *J*l l*U| I-l Ijl-Ull-I ,J*I >-Ul I |J) ) 
¥ 1 1 .J-l l*¥l (-1 «jl-V lj"l • J*l >’¥<| « j> J 


U< t , 
¥ l I 


♦Ut l- l «J-l »-Ul t-l .Jl-UlI.J- U 
♦Vlj-1 iJ-lI-VlI-ltJ)-¥ll«J-l J 


ipute AT bottom mu 

CONTIN' £ 

DPOX * CST2L •( PI 1*1 aO) 
D2UDX2* CSTL2 *1 UII-I.JJ 
0ZV0X2- CSTL2 *1 VII-I.JI 


DZUDF ' 
D2V0XT 

opdt 

02U0Y2 
D2YQY2 
GO TO 300 


CSTZL2 • I Ul »*l , 
CST2L2 •! ¥11*1 , 
S5TIL *1 PtI.J*| 
CSTL2 •( UlI,J*l 
CSTL2 • I Y ( I i J* l 


IPUTE AT TOP 
CONTINUE 
DPOX - CST2U 
D2U0X2* CSTU2 
D2¥DX2- C5TLZ 


MALL 


D2U0XY 
02VDXY 
DPOT 
02UDYZ 
D2Y0Y2 

CONTINUE 

RETURN 

END 


( PI 1*1 .Jl 
•I Ull-t.J) 
• I V 1 1 - 1 i J ) 


C5T2L2 •( U(U1, 
CSTZL2 *1 V{I*I, 

cstil •( eij.j) 

CSTU2 •( U 1 1 1 J- 1 
CSTL2 • I VtliJ-1 


- PII-I ,J! I 
-2-»U 1 1 • J I ♦UU-liJl » 

-Z..YU.JJ -VII*!, JJ I 
J-U-Ull-l ,J)-U(I-l,J*lJ-Ull-lU» 
j-l>*¥<|-liJ>-¥<I-iiJ*lJ-¥l!*liJ> 
I • PI l ,01 > 

J • 3, *U 1 1 i J ) J 
] - 3 r *¥ I I i J 1 I 


• PI I-I .JI I 
*2 r-UI I * Jl AUll-l.jl I 
•2.-VII.JI -Ytl-l.JJ l 
JJ*U ( Ul , J-l )-Ul I, l ,Jl-U 1 1- l , J-| J ! 
J|*V|I»l,J*ll-¥(l-I,J}.¥lI*l,J-l| J 
- P 1 1 . J-l I l 
I - 3 • *U I l • J 1 1 
) - 3**YlI,JJ I 


ghatohis 

GHaTQIYA 

GIIAT0N17 

g>m rot is 
gnatout 
ghatoisq 
ghaTdssi 

GHaT01S2 

GNAT0»53 

guaToisi 

5 U a T015S 

GNAT015A 

GNAT0157 

GUaTOISS 

snaTOis? 

ghatoyad 

gmatoiaj 

GHAT0YA2 

G’IaT01A3 

GNaTDIAI 

GNaTOIAS 

GUAT09A* 

G?AT0YA7 

GNaTOIAB 

gnatoiai 

SNAT0170 

GNAT0171 

G«aT0172 

GKATOY73 

GNAT0171 

GNAT0175 

GNAT017A 

GNaT0177 

GNAT0Y78 

GNAT0Y79 

GNaTOIsQ 

GNATOIgl 

G«aT0Y8! 

GHAT01B3 

SNAT0TB1 

GNaTOYSS 

GNaTOTBA 

GNaT01«7 

gnaToiss 


END OF UNI ¥Ac 1108 FORTRAN ¥ COMPILATION 


0 -0 1 ASNOST 1C' MESSAGE ( S I 



8 FOB DIFFZ, DIFF2 

UNIVAc HOB FORTRAN V LEVEL 2206 ODU FsOIBH 

This compilation ras done on 27 Jan 7i at i4j33)4* 


sn*tonbt 


27 JAN 71 


14(33145.982 


SUBROUTINE P|FFZ ENTRY POINT 000236 

storage used [Block, n*ne, length) 


OOOI 

•code 

□00275 

0000 

• data 

000032 

0002 

•blank 

000000 

0003 

DATA 

0071*4 

0004 

02 

000010 

0005 

N0DNAL 

0)0472 

000* 

TRNSMT 

000005 


external references (block, nahei 

ftOQ7 NERR3S 


00Q1 000203 )2o L 0001 00012* 2Ql 

0000 R 000001 CSTL2 0000 N QQ000O CST1L 

000* R OOOOOR 0S4 0001 « 000002 £,flSupX 

□OOR R OOOOOH dUdX 0004 R OOOOOS QVuT 

0000 | 000002 |1 0000 | 000003 j2 

0003 R 000031 L 0003 0000^0 L<«lTS 

0003 0007*4 R 0003 0040*4 RrTeS 

OOOS R 007*01 RET 0003 005524 RU 

0005 r 00|5*2 RVX OQOS r 002453 r V T 

0003 D02424 u OOOG r 003344 ux 


storage aSsignhent f or variables (bloc*. type, rela'ive location, namei 


oooi 10 
009220 
R 00000 1 
R 000003 
R OOOQO* 

J 0000D4 

1 000050 
000(34 
R OOQOOD 
R 005)2* 
003244 


10 L 

2 Q 0 L 

out 

OREVdT 

02 HOX 2 

J 1 

NG 

RbAR 

RUX 

R* 

V 


ooot 

0003 R 
000 * R 

0004 R 
0D04 R 

0000 j 

0003 1 
0003 

0005 R 
0005 R 
0005 R 


000143 IDOL 
000000 A 
000002 DOZ 
000000 ORUOX 
000007 d2hDY 2 
000005 j2 
000074 N 5 
0047Q4 RE 
□00*71 rut 
00*017 RY 
004235 VY 


0001 

0003 
000* R 

0004 R 
0003 R 
000 * ] 
0003 

0005 R 
0003 
0003 


0001*5 11QL 
000(32 A VE 
000003 0 g3 
00000 1 BRVOT 
00**04 H 
OOODOO K 
000024 PROF 
00*710 R E x 
00*344 RV 
000144 STATE 


00101 1 * 

00103 2 * 

00104 3 * 

00105 4 • 

DO 1 05 5* 

00 ) 0 * *■ 

0D107 7 • 

00110 B« 

OollO 1 * 

Oo I 1 0 10 * 

ooin ii* 

001*2 12 * 

00)|2 13 * 

00113 14 . 

430113 15 * 

00113 1 * * 


SUBROUTINE D1FF2 ( I t Jl 
COMMON /DATA/ A { 3 700 I 

COMMON /02/ ORUOX,ORVOTiOREUDE>DREVDT>DUDX,DVDY>D2hDX 2,DZHQT2 
COMMON /NOD*AL/ RUX (21 ,Z| I ,flUY(2l ,21 ) ,RVX 121 ,Z| I ,RYT (Zl ,21 ) 

1 ■ UXI2I ,2) ) ,VY(2l .21 t ,RX(2I ,21 1 ,RY(Z] .21 I ,REX(21 ,21 ) ,RET(2l .21 ) 
COMMON /TRNSMT / K,OQ1,002 iDG3,0(]4 

DIMENSION RATES(l*001i REI20.2Q], RU120.20), RVI20.20) 

ESU] VALENCE (a* 3 !?! : .RATES) , JR ATES 140 1 J i RE) 

1, (RATESCBOl I iRU) , RATES l 1 20)) .RV I 

2, (Ai21 ) .PROP I , (PRCP(*).L> 

EQUIVALENCE U(4' , ,_;n)T 5 I , (A(10) ) iSTATE) , { A(?IJ, AyE 1 
DIMENSION LIMITS ISO I ,NG I 20) >NS(20) (STATE 1 2000) . A V E ( 1 , PROP U 0 1 

1 , R(2D,Z0) ,H(20,20> ,U (20,20) ,V (20 ,20) 

EQUIVALENCE (L - M I TS | 1 | , NG ) , ( L I M I TS | 2 1 I ,NSj , [AVE[3| , Rti AR ) 

2, ISTATE ( -luil.HI 

3, (STATE ( 80l),H) 


GNaT 04?0 
GNAT04 V l 
GNAT0442 
GNAT0493 
GNAT04V4 
GNAI04V5 
GNAT04V* 
GNAT04V7 
GNAT049# 
GNAT049V 
GNATOSOO 
GNATQ50I 
GNAT05D2 
GNaT0S03 
GNa T0504 
GNAT050S 



00113 

17" 



H > ISTATEllZOl 1 .Ul 

00113 

l 8 * 



s, c«taTeU*ou t Vl 

00S1H 

l*" 



REAL L 

001 IS 

20" 



CSTIL * 1. / L 

0011* 

21" 



CSTL2 * CSTIL • CSTIl 

DOHA 

22" 

C 



00117 

23" 



11 ■ NGIJ1 

00120 

2H" 



12 ■ NS(J| 

00121 

25* 



Jl ■ NGtl 1 

00122 

2A" 



J2 • NS 1 1 ) 

00122 

27" 

C 



00123 

26* 



DUOX . (UX(1»1,U) - UX[I,Jj) /L 

0012'* 

2*" 



OVD* » IVY(I,J"11 - VTtl.JM /{. 

00125 

30* 



DROP* ■ t R X i I » 1 iJ)"UX(l*l tJ> - RXf 1 ,jl "UXII .J) 1 /L 

00 1 2 A 

31" 



DRVDr * IKTCl,J*U"Vr(l iJ"|l - RTtl«JI"VTIl»J)l /L 

00127 

32" 



DREUDX ■ (REX{]"1.J)"UX{]*I,J) • REX ( I , J )"UX ( I , J) I /L 

□ 01 30 

33" 



DHEVdT • IRET(liJ*|l"VT(I»J*|) - RETt I *Jl"VYl I i J) 1 /L 

00130 

3H" 

C 



00131 

3b" 



IF ( | .E8.I 1 1 GO TO lo 

00133 

34" 



IF 11. EG. I2J GO TO 20 

00135 

37" 



P2HDX2" CSTL2 •( - 2 . • H < 1 . J > »M(t"l»Jl 1 

00I3A 

38* 



GO T ■ 100 

ODl 37 

3* • 


10 

P2.NDX2 • CSTIL *( CSTIL •( H U • 1 , J) -HI 1 , J 1 1 * DBl/K » 

001H0 

HO" 



GO TO 100 

001 HI 

HI" 


20 

P2HDX2 * CSTIL • ( -D02/K - CSTjL •( tt(I.J) - Hll-t.Jl 1 

001H1 

N2" 

c 



001*2 

83" 


100 

CONTINUE 

00 1 H 3 

HH • 



IF I3.EB.JI 1 GO TO 110 

001H5 

NS* 



IF <J'E0.J21 GO TO 120 

001H7 

NA* 



02KDTZ " CSTLZ •( Hl|iJ*tl «2."H(liJl • Hl|.J»ll 1 

00150 

N7" 



50 TO 200 

00151 

18" 


lie 

D2HDTZ " CSTIL «( CSTIL "1 H ( 1 , J" 1 1<1( I , J > 1 * DB3/* 1 

00152 

HS" 



GO TO ZOO 

00153 

50" 


120 

D7H0T2 * CSTIl, "1-pBH/R - fSTIL "1 Hll.JI * 1 

(JDlSH 

51* 


200 

CONTINUE 

00155 

SZ" 



RETURN 

0015A 

53" 



ENP 


ENO OF UNIYAC 1108 FORTRAN V COHWlLATJON 


□ •OUfiNOSTIC* MESSAGE I S 1 


GNaTOSOA 

6 N A T0S07 

sn a toso* 

GRaTOSD* 
SNAT0S10 
GNATG5 1 ) 
SNAT0&1Z 
GNATQ51 J 
fiNAT05lH 

sratosjs 
ONATOS l * 
GNAT0517 
GNaTOSI# 

gnatosi? 

GNATOSZO 

GRaTOSZI 

gN A T0522 

GNAT0523 

GNATC52R 

SNATD52S 

GNaTDSZ* 

GNATDSZ7 

GNAT0S28 

GNAT052? 

G»AToS30 

GNATDS3I 

GNAT05JZ 
GNaT0533 
GNAT0S3N 
GNAT0S3S 
GNAT053 A 

GNAT0537 

GNAT0S3# 

GN*T053* 

GNAT05N0 

GNATOSNl 

ghatosnz 



0 FOR HEATER 

UNIVAC HOB FORTRAN V LEVEL, 220* 0018 FSQ13H 

This compilation has done on 27 jan 7i at i<h33;nb 


GNAT05H3 


27 JAN 7 | 


In;33:nb, a? 


SUBROUTINE heater 

ENTRY POINT 00302Q 

STdRaGE USED 

(BLOCK. 

NAME, LENGTH) 

0001 

•code 

000022 

0000 

•data 

ooooio 

0002 

•blank 

oooooo 

0003 

dstate 

003)00 


external references (Block, nanej 

□ ODH NERR3* 


storage assignment for Variables iblock, type, relative location, namei 

0003 R 000000 OR 0003 R 000*20 ORE 0003 R OOtRHO ORU 0003 R 0022*0 o«V 0000 

0000 l 00000! J 


AO 

CTv 


00101 

!• 


SUBROUTINE HEaTER l OOHEAT. L3 1 

GNAT05HM 

00103 

2 • 


REAL L3 

GNAT05MS 

OotOH 

3» 


COMMON /DSTATE/ DRI2o,2D), 0Re(20,20I. ORUI2Q,201, DRVI20.20) 

GNaTOSM* 

□ Ot ON 

H • 

C 


GNAT05H7 

0010M 

S* 

c 

insert heat at SELECTED nodes 

GNaTDSVS 

00 1 OR 

*• 

c 


GNaTOSMV 

00105 

7 • 


l * 12 

GNATOSSD 

0010* 

8* 


J ■ 10 

GNAT05S1 

00107 

9 * 


DREUjJ) ■ ORE 1 1 ,-J t * D8NEAT / L3 

GN3TDS52 

00110 

10. 


[00 CONTINUE 

GNAT0553 

Oolll 

U» 


return 

GNaTOSSV 

00112 

12* 


END 

GNaTDSSS 


OOOQOQ 1 


END OF UNIVAC 1108 FORTRAN V C0H81l*T10N 


0 ‘DIAGNOSTIC* MESSAGE IS 1 



B FOR BCOUT, BCOUT 

UN i VAC I t OB FORTRAN V LEVEL 2204 00 1 B FS018H 
THIS COMPILATION has DONE ON 27 JAN 7t AT 19143119 


ON AT q SB 4 


SUBROUTINE bcooT 

ENTRY POINT 000111 

STORAGE USES 

'BLOCK p 

NaHEi LENGTH) 

0001 

♦COOE 

000122 

0000 

•data 

000017 

0002 

•BLANK 

OOOOOO 

0003 

DATA 

007 | 48 

0001 

dstate 

003)00 


EXTERNAL REFERENCES (Bl *K , nAhEJ 
OOOS NERR3I 


storage assignment for variables (block, ttpei relative location, nahej 

COOl 000015 2L 0003 R OOOOOO A 0003 000132 AVE <303 003000 CTL 

0001 R 000420 ORE 0001 R OOlHO QRU DOOM R 002240 PRV 0000 l OOOOOO t 

0003 R 000031 L 0003 p 000132 PBAR 0003 000021 PROF 0003 R 000741 R 

\D 0003 p 000131 RBAR 0003 R 001701 RE 0003 p 000134 REBAR 0001 r 00552H RU 

" J 0003 0001 11 STATE 0000 R 000002 UKALL 


oatoi 

1 • 

SUBROUTINE BCOUT IN, HOOT) 

GNAT0S57 

SO l 03 

2* 

REAL l 

GNaTDSSB 

□0101 

3» 

COMMON /DATA/ AI37O0) 

qnaTOssv 

OQIQS 

1* 

COMMON /OST £TE/ 0R(20i20), DRE!20 i201, 0RU(2Q,201. 0RV120i201 

GNAT0540 

00)04 

5* 

dimension AVE 1 1 0 1 , RATES ( 1400 > , STATE12000) 

GNATD541 

00)07 

*• 

0 1 MENS 1 ON CTL 1 20 ) t PROrUOl 

QNAT0542 

□ 0110 

7 • 

EQUIVALENCE 

SNAT0543 

□0110 

B* 

1 (Alll.CTLI 

GNAT0541 

□0110 

?• 

2, ( A ( 2 1 ) iPRQP I 

GNATDB4S 

00110 

I0» 

1, IAtfl),AVEl 

QNA T0S44 

00)10 

n* 

5, t * 1 1 o I J i state: | 

GNAT0S47 

00110 

12‘ 

4, (A12101 1 , RATES! 

GNAT0S48 

00111 

1 3 • 

EQUIVALENCE 1 PROP 141 ,L 1 

GNAT0549 

00112 

H* 

DIMENSION R ( 2Q i 20 1 

SKAT0S7D 

Da) |3 

15« 

DIMENSION RE 1 20 1 20 ) , RU(20,2D), RVUQ.20I 

GNAT0571 

001 11 

14» 

EQUIVALENCE 

GNAT0S72 

oom 

1 ?• 

2 ISTATeIIQII ,R1 

GNAT0S73 

O01 11 

1B» 

7, (RATES (101 1 , RE 1 

GNAT057 1 

00111 

19* 

8, ( RATE; I SOI 1 t PU ) 

GNAT0S7S 

001)1 

20* 

9, 1 RATES (12011 ,RV 1 

GNAT0574 

00115 

21* 

DOUBLE PRECISION PBAr, RBAR, rebar 

GNAT0S77 

^0116 

22* 

EQUIVALENCE IaVEU) iPBAR) i (AVE(3l ,RBaR> « UVE (S J , REBAR 1 

GNAT0S78 

aoi 17 

23* 

I « 20 

CNATDS>'9 

00120 

21* 

J ■ 10 

GNaTOSBO 


27 JAN 7 j 


l*l»33I‘t».*li2 


000* R OOOOOO OR 
oooo i DQOoai j 

0033 00*104*! R*TeS 

0003 r 00431*! RV 



00 121 

25* 


IF IN.E9.2} CO TO 2 

CNAT05B ) 

00123 

26* 

1 

CONTINUE 

CNAT05B2 

00121 

27* 


UlALt. * HOOT / IRt I tjl *RBARt /L/L 

G N ATO 5 B3 

00 125 

28* 


DRU l ; ■ J 1 ■ DRU(ltJ) - RU(J.J) • U* ALL /L 

GNAT05B1 

00|24 

29* 


DRV(I.J) • DRV (1 p J 1 * RV(I,JI • U* ALL /L 

GNAT05BS 

00122 

30* 


RETURN 

CHAT 05B4 

00130 

31* 

2 

CONTINUE 

, ,T05B7 

00131 

32* 


ORlJu’t - OR 1 1 1 3 1 * ( R 1 1 . JI*RbaR1 * URaLL /L 

CN*T05i8 

00132 

33* 


DrT- ( 1 t J ] • Dr£{ 1 , J| . (RE( I 1 Jl*REBAR*pBAR | • UftALL /L 

CNAT05B9 

00133 

31* 


return 

CNAT0S9Q 

00131 

35* 


END 

GNAT059 1 


END OF 

UNI VAC 

1104 FORTRAN V COMPILATION. 0 *OIaCNOSTIC* M£5f ACE { 5 ) 



vO 

CO 



B FOR TEMp, TEMr 

E’KIVac uoa FORTRAN V level 220* 001? 
t»is compilation sas done on 27 v u- T > 


6NAT0&V2 


3 ■ SO 


FUNCTION Temp ENTRY POINT OOII..J 

STORAGE USED (BLOCK, NAME t LENGTH) 

ODD) • CODE 000112 

0000 *DATA 000020 

0002 ‘BLANK OOOOOO 

0000 TEFCTN 00OQ7N 


External references (block, nahei 
DOOR NE«R3s 


storage assignment for variables (block, type, relative location, nahej 

0001 0000 1 A IOL 0001 000030 20L 0001 0000*1 SOL 0001 000070 MOL 

0000 R 000001 ElN 0000 ) 000002 ] 0000 I 0000D1 J S 0000 | 000003 u 

0000 R ODOOOT TEMP 

VO 

VO 


00101 

l« 

FUNCTION 

TEMP ( ERHQ, RHO 

1 




SNAT0593 

00103 

2* 

COMHON /TEFCTN/ T130J, E(30l 




GNAT0591 

0G|Q1 

3" 

DATA E / 

*59. 1*7 i»5l .2*3*- 

13*723 

31*915 t 

-2* *02/ 1*15.997 • 

■10*023, 

GNAT0S95 

00101 

H * 

1 -*,1*7, 

•2*0*8, 1,801, 

| A « V7B , 

22*8*5 i 

2* *302 I 31*00*1 

31.5531 

GNAT059* 

00 1 OR 

5* 

2 37,515, 

10*220* 12**87, 

15*00* , 

19. 311 * 

53*1001 57,273, 

*1*011, 

GNAT0597 

00101 

*» 

3 *1,457 , 

*8*225 i 71,731* 

75*178, 

76, *27, 

2*0i/ 


0NAT0S98 

0010* 

7 • 

DATA T / 

1*0*1 l*Oit 

200, i 

220« i 

210,1 2*0*i 

270. , 

GXATD599 

0010* 

8* 

1 275,, 

280*1 2t5*i 

2 90,, 

295.* 

3Q0 i i 3l0*i 

320* • 

GNATDAQO 

0010* 

9* 

2 330,i 

310., 350,, 

3*0,, 

380, , 

100., 120,, 

110,, 

gnato*gi 

0010* 

10* 

3 1*0* i 

180* * 500 t , 

520 • • 

SlO*i 

2*0./ 


GNaTOAOR 

00110 

11* 

DATA 1, 

N / 10* 28 / 





GNAT0A03 

00110 

12* 

C 






GNA T 0* 0 1 

00113 

13* 

ElN - ERHQ / RKO / ,2SQ5*AE5 




GNaTDADS 

OOIIH 

lit 

IF ( EIN 

*GT. E( 1 1 ) GO TO 

to 




GNATOAO* 

0011* 

IS* 

TEMP *T(ll 





GNATQA07 

0011? 

It* 

RETURN 






GNaTOAOB 

00120 

17* 

10 IF ( ElN 

• LT • EIN! ) GO TO 

20 




GNaT0*09 

00122 

1 S V 

TEMP ' TIN! 





GNATQA’O 

00123 

IV. 

RETURN 






gnato* 1 { 

00121 

20* 

20 CONTINUE 






GNATO* ] 2 

0012S 

21* 

IF I EIN 

• LT* Ed) 1 GO TO 

30 




ff*-„T0*l3 

00127 

22* 

IF ( EIN 

• GT* El 1*1) 1 GO 

TO 10 




gnatoah 

00 1 27 

23* 

C INTERPOLATE 

TEMPERATURE 





GNATO*[S 

00131 

21* 

TEMP « Till * (EIN-EUI) • 

(T(l*ll 

-T 1 I ) 1 

/ (El 1*1 )-E( 1 1 ) 


gnaTO* i * 

— 110132 

25* 

RETURN 






GNAT0AI7 

Ob 133 

2** 

30 1*1-1 






GNAT0A18 

00131 

27* 

GO TO 20 






GNAT0AI9 


27 jan 7| 


im;33:so>*Fq 


0003 R 00003* c 
0003 R OOOOOO I 



DO 1 3S 

2B * 

BO l*l*| 

GNAT0A2Q 

00136 

29* 

GO 10 20 

GNAT0621 

00137 

30* 

END 

GNAT0A22 


END or UNIVAC 11QB FORTRAN V CORE ILAT ION 


0 ^DIAGNOSTIC* MESSAGE I S 1 



TOT 


B FOR beta, beta 

UHIYAC HOB FORTRAN V LEVEL 220* 0018 FSOJBh 
This compilation has done on 27 jan 7i at ir;33,-S2 


GNA t D623 


22 JAN TI 


2 » 


FUNCTION BETA ENTRY POINT 00QQ7S 

STORAGE USED (SLOCK, NAME, LENGTH] 

0001 • CODE 000106 

0000 *DAT A 000016 

0002 -BLANK 000000 

0003 TEFCTN 000078 


EStCRNAL REFERENCES (BLOCK, NAME 1 
0008 NERR3S 


storage assignment for Variables c b lo c k. type* relative location, name) 

0001 00001 8 10L 0001 000080 I l ?G 0001 000032 20L 0001 OODOSQ 80L 0000 R OOOOOO BETa 

0003 fl 00003* E 0000 [ 000003 { 0000 ( QDOOOI M 0000 | 300002 Kl 0003 r OOOOOO T 


00101 

1 • 

FUNCTION BETA (TO I 

GNAT0A2R 

00103 

2* 

COMMON /TEFCTN/ Tt30I, E130I 

GNAT062S 

00108 

3* 

M • 28 

GNAT06Z* 

□ OIOS 

8* 

IF ( TO ,GT« TUI ) GO TD 10 

GN AT06Z7 

00107 

S* 

BETA - TUI 

GNAT062B 

001 10 

6 • 

RETURN 

G'UT0*2? 

QQl H 

t ~ 

10 IF I TO * LT * TtMl 1 GO TO 20 

GNAT063D 

00113 


BETA * T1M1 

GKATOii: 

DO] 18 

Vt 

RETURN 

GNAT0632 

001 18 

to- 

C FIND I SUCH THAT TUI *LE* TO *lE. TU*|J 

GNAT0633 

0011S 

il* 

20 HI » M-l 

SNAT063R 

□ Oil* 

12* 

00 30 1-|, HI 

GNAT063S 

00U1 

13* 

30 IF ( Tll-11 «GE» TO } GO TO 80 

GNAT0636 

00128 

IR- 

80 aETA - T ( 1 1 ♦ (Ell-ll-EUll • (T 0 -T(I 1 »/|T|I* 11 -TUI 1 

SHAT 0*3 7 

00I2S 

IS* 

RETURN 

GHATD63B 

30] 26 

16* 

END 

GIUT063? 


END OF UNIVAC 1 1 OB FORTRAN V COMPILATION, 


0 •DIAGNOSTIC* MESSAGE ( S ) 



102 


B FOR PRESS, PRE55 

UN I V AC l I OB FORTRAN V LEVEL 2204 0018 FSQ18H 

this compilation »as done on 27 jan 71 at 1*1:33:53 


SHAWANO 


27 j^n 7 l 


in:oj;S3.3&3 


SUBROUTINE press 

entry point ooooir 

storage used 

(BLOCK, 

NAME, LENGTH) 

0001 

•code 

000030 

0000 

•data 

□OODIO 

0002 

•blank 

000000 

0003 

DATA 

Do ? | am 


EXTERNAL REFERENCES {BLOCK, NAHE) 
ODOR NERR3S 


storage assignment f°R variables {block, type, relative location, nanei 

0003 R OOOOOD A 0003 000132 AVE 0003 0 000132 PBAR 0000 0 000000 PP 0003 00002m PROP 

0003 D 0001 BN RfiAR 0003 R OOOB32 RCONST 0003 R 000030 2 


□ 0101 

l* 

SUBROUTINE press { P, R, HM 1 

CNAT04T1 

00103 

2 • 

COHMON /DATA/ AI3730] 

GNAT04 12 

OOIOH 

3* 

double precision pbak, rear, pp, hh 

GNaTOA 13 

00105 

*)• 

DIMENSION PROP ( 1 D 1 , AVE ( 1 0 1 

onaTOamh 

00104 

5? 

EQUIVALENCE ( A (2 ] 1 .PROP) , ( A { 9 1 1 , A VE 1 

GNAT04MS 

00107 

4» 

EQUIVALENCE t A VE { 1 1 i PB AR 1 ■ IAVEC3) ,RBAR» 

CNaTDAMA 

ooito 

7* 

EQUIVALENCE tPR0P(5),Z|, (PRop(7l, RCONST 1 

ONAl 04R7 

001 1 1 

8« 

PP « 2* !RBAR*r) • RCONST • HH 

gnaTDasb 

00112 

V» 

P • PP • PBAR 

CNATOAMV 

00113 

10 * 

return 

GNAT04S0 

001 1 H 

{!• 

END 

GNaTOAS 1 


ENO OF ONIVAC life FORTRAN V COMPILATION. 0 .DIAGNOSTIC* MESSAGEC5J 



B fOR RTPRES, RTpR E S 

i/Miv^c i iob Fortran y l£v£e 2206 oojb fsoibh 
this compilation was done on 27 jan 71 at n;33;s*t 


ON AT0AS2 


function htpres EnTRT POINT 000210 

storage used ;block> na'E, length) 

0001 ‘CODE 0002IC 

ODDS .DAT A 0001 10 

0002 •BLANK 000000 


ExT £ RN*L references IBloCk, n*mE) 

0003 NEXP6S 
0001 EXP 
0005 NERH3S 


STORAGE ASSIGNMENT FOR VARIABLES IBLOCpc, TrpEi RELATIVE LOCATION. nA m E) 

OOOD ( 000011 IS 0000 R 00000*1 N 0000 R 000010 H25p2 0000 R 0000^2 dDt! 

□000 R ODODSO P 0000 R 0000*3 PHPC2 0000 R 0000*2 PmPC 2 S 0000 R 0000*1 PN2» 

0000 R 000056 P2 0000 R 000057 p3 0000 R 0000*0 PI 0000 H OOOOD3 R 

0000 R OODOQD RTPRES 0000 R 00005) T 0000 R 0C00Q2 TC OOOU R 000051 TMTC 

S-* 0000 R 000015 TNlO 0000 R 00001* TN2127 0000 R 00001? TrZ* 0000 R 000013 TNI 

p 0000 R 000012 TN* 


00 101 

J* 


FUNCTION RTPRES IRHO.TO) 




GNAT0653 

ooioi 

2* 

c 







0NAT0651 

00101 

3* 

c 

THIS ROUTINE CALCULATES STE*A«T'S EQUATION OF STaTe 

for oxtgengnaToass 

OOIOI 

1» 

c 







GNATG656 

00101 

5« 

c 

INPUT 

• 





6HAT06S7 

00)01 

6 * 

c 

RHO 

• 

DENSITY 

- 

LB / CU-FT 


6NAT0658 

OOIOI 

7 • 

c 

TO 


TEHPERATORE 

- 

DEGREES RANKJN 


5NaT0659 

OOIOI 

8* 

c 

OUTPUT 

- 





GNAT0660 

□ 0101 

V« 

c 

RTPRES 

• 

PRESSURE 

• 

PS I A 


GNaTD*61 

OOIOI 

10* 

c 

RTOPDT 

«* 

PARTIAL OF PRESSURE 

• 

PS I A / 0-R 


GNAT0662 

OOIOI 

1 1 • 

c 



WITH RESPECT TO 




G1ATC663 

OOIOI 

|2. 

c 



TEMPERATURE 




GNAT0661 

OOIOI 

13. 

c 

RTDPOR 

• 

PARTIAL OF PRESSURE 

• 

psia / ilb/cu«fti 


Gl*T0A65 

00101 

11. 

c 



WITH RESPECT TO 




GNAT0666 

001U1 

15. 

c 



OENSITT 




GNAT0667 

00101 

16. 

c 

CONSTANTS 

i . 




GNATD66B 

OOIOI 

17* 

c 

RHOC 

• 

CRITICAL density 


G.MOL f LITER 


GNaT0**9 

OOIOI 

ie. 

c 

TC 


CRITICAL TEHPE-RatPRF 

- 

DEGREES KELVIN 


GN a T 0670 

□ 0101 

IV* 

c 

R 

- 

GAS CONSTANT 


liter-atm / g-hOl-o 

KELT J N 

G"t A T067I 

□ 010 1 

20. 

c 







GNATD672 

00103 

21* 


OATa RHOC 

, 

TC , R / 




G1AT0673 

00103 



• 13.333. 

151.77, Q.0820S35 / 




GNAT0671 

00103 

23. 

c 

COEFFICIENTS FOR THE EOOatION 

OF STATE FOR UXTsEN 


GNAT0675 


27 JAN 7 | 


11 *33*51 *60* 


OQOO R 0C00 5 3 O0T2 
0000 R DUD0*1 POUT 
DOC J R ODODOt RHOC 
OOOD H 0000*5 TmTcZ 
0000 r 000011 TN 5 



00107 

2H* 


REAL NUB) , N25P2 

GNAT0674 

00110 

25* 


data n 

/ 

GNAT0677 

Oat io 

26* 


1 

3.3BJS90BE-3. -1 .31 A0A22E*0. "7 *3SS2852E*3 i 

GN A T 06 7 8 

oat id 

27* 


H 

I .920H907E*7 . -2.90260016*10 .-S. 701 OllH'ti 

GNAT0479 

Dot to 

26* 


7 

7.9A82238£*5i 6 .07O22SQE-3 , -2.7 l 0 l 964£*0 » 

GNAT0680 

DO 1 1 0 

29* 


0 

-3.S?91960E*1 , l .02ZO95AE-6, 1 , 90H5HS I E-H . 

SNAT0681 

□at to 

30* 


1 

1 .21708396-5. 2. H H25S9SE-3 ■ 1 . 736555 IE*2. 

CNAT0682 

ootto 

3 l * 


6 

3> 01 7528HE*5, -3.H952B52L*7. 8 , 8472H00E" I t 

GNaTD683 

ootio 

32* 


9 

-2 *6791 7A7e*2 • 1 .o547o7qE*5 . 5*6377 10 S E*3 * 
•1 • 12012?|t*0. 1 .H4829HVE*2. 9 . 9sa4e?2£-H i 

gnato* 8 h 

001 ID 

33* 


2 

GNAT0485 

001 to 

3H* 


S 

-0*00540 . -D* 157 . -0.35Q » 0.90 t 

GNATD486 

aoiu 

35* 


DATA TNG . TNHi TNS « TNIO. TN2H27 ■ TN2A / - 1 1 • H020232E-S , 

CNATQ6B7 

00112 

36* 


• 3 • 8^0981 N£*7 . *8 . 7078E03E* 10 i "7 • 1 8B3920E+1 t -4.99Z082HHE-H i 

GNAT06B6 

00112 

37* 


* 1 H 

t 

GNA T 06 8 9 

00112 

38* 

C 

convert the input fRoh lb/cu ft and degrees nankin to 

GNaTD490 

00112 

39* 

c 


G-MOL/LITER and DEGREES KELVIN 

GN*T069| 

D0121 

HO* 


P * RHO 

■ (500575 • 32.2 

GNAI0492 

00122 

HI* 


T ■ TO 

* .55555554 

GNAT0493 

00123 

H2* 


OOT l « 

1 * / T 

GNATQ69H 

D012H 

H3* 


OUT 2 * 

ODTl • QDTl 

GNAT069S 

00125 

HH* 


TKTC • 

T * TC 

GNAT069A 

OOUA 

HS* 


TNTCZ « 

THTC * THTC 

6NAT0697 

00127 

HA* 


P2 - 

P • P 

HNAT0495 

00130 

H7* 


P3 « 

P2 • P 

- V,TQ499 

00131 

HO* 


PH ■ 

P3 • P 

SNAT070C 

00132 

H9 • 


N2SP2 ■ 

' NUSI-P2 

GNAT070 l 

00133 

50* 


PN28 - 

P **N 1 28 ) 

GNAT0702 

0013N 

SI* 


PHPcaa 

» PN28 - RH0C**N 128) 

GNAT07Q3 

00135 

52* 


pftpra » 

' PHPC28 • PHPCzB 

SNAT070H 

00135 

53* 

c 



SNAT0TU5 

00135 

5H* 

c 

FA*. 

JVATR The EOUATIOH OF STATE 

GNaTETPA 

00135 

55* 

c 



GRAT0TQ7 

0013A 

56* 


pp«r * 

P*R*T 

aruTD7ci 

OOUA 

57* 


c 

* P2* <N ( I ) *T*N III* INI Jit (NlH l*N(Sl*OoT2)*OoT2)*OoT2)nNAT 0709 

00136 

58* 


c 

♦ P3*( In I A 1 *T*Nt7 1 )*T*NCa)*CN«9)*- 1 10 1 *O0T 1 1 «0DT l 1 

6N3 t 07 ' r 

00I3A 

59* 


c 

* PH*«nIUI*T*NI|2)*P*(nU31*nUHi*0DT|11 

ENATO' ) 

00134 

40* 


c 

* P3»ODTZ*E*ptN25pZl • l 

6*00712 

00136 

41 • 


c 

IN! IS)* [Nl I 6 l *Nl 1 7 1 *ODT 1 }*ODT> > * 

5**70713 

DOHA 

62* 


c 

1 iNlJCUlNl 19I*NC20)*00T1 J*ODTl 1 * 

GNAfO' 'H 

00136 

A3* 


c 

(NI2I )*IN(Z2I*n( 23)*QDTI l *00T l ) *p2 1 *P2 I 

GLA.T07 1 5 

00136 

SH* 


c 

« NI2H | •P*Pn2S*phPC28*EXP<N(26I *PHPC2*N(27 |*TNTCZ) 

S' .TC# ; A 

00136 

65* 

c 



..k* T07 1 / 

00137 

A6* 


RTPRES 

* POUT * | H .496 • 1 HH , 

GN,T07‘8 

00 I HO 

A7 • 


return 


„N£7071 9 

OQ | hi 

AS* 


END 


GNAI0723 


END OF UN l VAC UOB FORTRAN V COHBlLATION. 
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S FOR BULK, BULK 

UNlVAC I10B FORTRAN V LEVEL 2204 COIB F50 |Bh 
this compilation mas done on zi jan 7i at i*t;33:Si 


GNATD721 


27 Jan 7 1 


i*;3 


SUBROUTINE BULK ENTRY POINT D0022I 

STORAGE USED ] BLOCK , NAME, LENGTH) 

0001 »COOE 00023Q 

0000 * DATA 000044 

0002 »BLaNK 000000 

0003 DATA 007)44 


EXTERNAL REFERENCES IBLOCK, NAME) 

0004 TEMP 

0005 RTPfiES 

0006 NERR3S 


STORAGE 

ASSIGNHENT for 

VARIABLES 

(BLOCK, 

, TrPE, 

RELATIVE 

LOCATION, 

, NAME) 









ODD) 


000015 

I27G 

□ 001 


00003d 

1 34G 

0001 


000070 

1 Sq G 

0001 


DO 0 1 0 3 

1 55 G 

ooo 3 

« 

OOOBDU 

A 

0003 


300132 

AVE 

0003 


oooooo 

ctl 

oooo 

R 

000012 

opbar 

□ 000 

R 

0000)3 

ORBAR 

oooo 

R 

000014 

ORErar 

0003 

R 

001404 

H 

□ COD 

R 

000017 

hbeR 

0003 

R 

000044 

MBARl 

0003 

R 

000Q45 

KM A* 

0003 

R 

000043 

MMIN 

0000 

1 

000004 

I 

□ 000 

I 

ooooi | 

J 

oooo 

I 

000007 

Jl 

□ 000 

1 

000010 

J2 

0003 

R 

000031 

L 

0003 

H 

000034 

label 

0003 


000050 

LIMITS 

□ 003 

I 

000121 

NF 

0003 

1 

ooooso 

NG 

0003 

R 

000)22 

NODES 

0003 

I 

000074 

NS 

0003 

I 

000120 

NO 

0003 

R 

000144 

p 

0003 

0 

000132 

PBAR 

0003 

ft 

000044 

PCOL 

0003 


000024 

PROP 

oooo 

0 

000002 

PT 

□ 003 

R 

000744 

R 

0003 


004044 

rates 

OD03 

0 

□00134 

rbar 

0000 

R 

000015 

R8ARI 

0003 

R 

004704 

RE 

0003 

D 

000)34 

REBAR 

OOOO 

fl 

0030 1 4 

resari 

OOOO 

0 

□00004 

ret 

0000 

0 

000000 

RT 

0005 

R 

OOOOOO 

RTPHES 

□ 003 


000144 

state 

0004 

R 

OOOOOO 

temp 

0003 

ft 

300047 

«T 

0003 


000024 

WTH 

0003 


000030 

2 














□ 0101 

l* 

subroutine bvvk 

GNAT0722 

00103 

2* 

COMMON /DATA/ AI3700) 

SMAT0723 

001D4 

3 * 

REAL NODES, L, LABEL 

GNAT 07 2 4 

0D1Q5 

4 • 

EQUIVALENCE 

GNAT0725 

00105 

5* 

1 

(Alt), CTL) 

5 N A T0724 

00IQ5 

4* 

2, 

IA(2l 1 .PROP) 

GNAT0727 

00105 

7* 

3, 

IAI31 1 ( L*3EL1 

SNAT0728 

□ 0105 

B * 

3, 

1 A 14 I ) .LIMITS) 

GNAT072? 

00 I 05 

9* 

4, 

IAI?1 ) , AVE) 

GHATO 730 

□ OIOS 

10* 

5, 

(AI]Q)).STAtE) 

GHAT073I 

□ 0105 

11* 

* , 

(A (2|Q1 ) .RATES) 

GNAT 07 32 

00 10* 

12* 

DIM 

ENSION Pt2C,20l ,Rt 20,20 » ,H 120,20 ) , RE *20- 201 

GHAT0733 

00107 

13* 

EQUIVALENCE 

SNAT0734 

00107 

14* 

1 

(STATES II, P) 

GNAT0735 

-J30107 

IS* 

2, 

SSTATEC4D1 1 ,Rl 

GNAT0704 

00107 

14* 

3, 

(STATE (SOI ) ,H) 

GNAT0737 

□ 0137 

17* 

7, 

IRATES14011, RE 1 

GNAT0738 

00110 

|B* 

DIMENSION *VE I l 0 1 > HaTESU*QO>, STATE1200QJ 

GNAT0739 
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001 1 t 

19* 

double precision pbar, rbar 

i REBAR, HT, FT, RET 

GNAT0740 

00112 

20* 

EQUIVALENCE ( AVEI 1 1 iPBAR) . 

t AVE (3 1 ,RBAR ) * < AVE(S) iREBAR) 

GNA10741 

00(13 

21* 

dimension label ( lot 


GNAT0742 

out n 

22* 

equivalence 


GN A TD7 43 

DOUR 

23* 

6 (LA8ELI6) .HH1N1 


GNAT0744 

001 IN 

24* 

7. (LABEL 17 ) , HBAR 1 1 


GN*T074b 

0D1 1 4 

25* 

B, (LABELIBI ,HNAX> 


GNaTD746 

00| |4 

26* 

», ILABELI’ J ,PC0U 


GNAT0747 

001 I 4 

27* 

*. ILABELI 101 >«TI 


GNAT07H6 

00115 

2B* 

0 1 MENS ION LIMITS 1501 


GNAT0749 

001 1 6 

29 • 

DIMENSION NG { 20 ) , NS (201 ■ PRQP(ID) 

GNAT0750 

00117 

30* 

EQUIVALENCE (LIMITS (ITiNGl, 

(Limits t2i i ,nsi 

GNAT075I 

001|7 

31 • 

1 . (LIMITS! 41 | # NO J , 1L1MI TS ( 42) ,NF) , (LIMITS (431 iNOOESl 

GNATD752 

00120 

32* 

EQUIVALENCE (pRDP (111 ATM ) ■ 

(PR0P(5),ZI, (PROP ( 6 1 ,L 1 

GNAT075J 

Ooizo 

33* 

C 


GNAT0754 

00121 

34* 

PT * 0. 


GNATD755 

00122 

35* 

RT - a. 


GNAT07S6 

00123 

36* 

RET > 0. 


5NAT0757 

OOlZ* 

37* 

HN ! N * 1000, 


GNAT 0758 

00|25 

38* 

HMAX **1000, 


GNaT07S9 

00(26 

39* 

do lao 1 *no i nf 


GMAT076Q 

00(31 

40* 

Jl * t(G t t ) 


GMATQ76 | 

00132 

41* 

J2 » NS ( | ) 


6NATQ7A2 

00 | *3 

42* 

DO 100 J«Jl,J2 


GNaT0763 

00(34 

*.1* 

PT ■ PT * Ptl,Jl 


GNAT0764 

00137 

44* 

RT * RT » R( l ,J1 


GNATQ765 

00140 

45* 

RET » RET * RE ( 1 i 3 1 


GNAT0766 

00141 

46* 

100 CONTINUE 


GNAT0767 

00144 

'17* 

DPBAR * PT 7 NODES 


GNAT0768 

00145 

48* 

DRBAR ■ RT / NODES 


GNATQ769 

00146 

49* 

DREbAR • ret / NODES 


GNAT0770 

00147 

50* 

00 200 l*M0,NF 


GN*T077l 

00152 

51* 

Jl * Noll) 


GNAT0772 

00153 

52* 

J2 - NSi t 1 


GNAT0773 

00154 

53* 

DO 200 J*JI > J2 


GNATD774 

00157 

54* 

P( I i Jl * P( 1 iJI * d pbar 


GNAT0775 

00160 

55* 

H 1 1 ,JJ » R I I |J1 - ORBAR 


GNATQ774 

00161 

56* 

RE 1 1 , J } * RE ( J , J ] ■ DREBAR 


GNAT0777 

00162 

57* 

HMlN « AHINKhMIN, M(1»J) 1 


GNATQ77B 

00163 

58* 

MM AX ■ AmAXMHHAXi MC1.J) 1 


GMATQ779 

00164 

59* 

200 CONTINUE 


GMAT0780 

00167 

60* 

PBAR * PBAR * DPBAR 


GNATC7B 1 

00170 

61 * 

HBAR » RBAR • DRBAR 


GHATQ7B2 

00171 

62* 

REBAR ■ REBAR * DREBAR 


GMAT07B3 

00172 

63* 

RBAR1 • RBAR 


GNAT0784 

00173 

64* 

REBAR1 « REBAR 


SNATD7B5 

00174 

65* 

HBAR - TEMP ( REBAR ) , RBAR 1 

i 

GNATD706 

00175 

66* 

HBAR | « HBAR 


GNa73787 

00]76 

67* 

PCOL • RTPRES : RBAR 1 , HBAR 1 1 / 144. 

GNATD7B8 

00177 

68* 

AT • 32,2 • HBAR • NOOES * 

L**3 

GNAT07B9 

□ 0200 

69* 

return 


GNAT0790 

00201 

70* 

END 


GMAT079 1 


ENQ OF UN|VAC IlOS FORTRAN v COMPILATION 
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a FOR OUTPUT. OUTPUT 

Univac tics Fortran v level zzoa ooib fsoibh 
this compilation was done on 27 jan 71 at tmsisa 


GNAT0792 


27 J*N 71 


l*t;32:S8* NR 


SUBROUTINE 

OUTPUT 

ENTRY POINT 000110 

STORAGE USED (BLOCK. 

NAME, length! 

0001 

•CODE 

00D152 

0000 

• data 

000017 

0002 

•BLANK 

000000 

0003 

DATA 

GD7I&H 


EXTERNAL REFERENCES (SLOCK, NAME) 

DOOR 

TAPEPR 

0005 

TAPEWR 

0004 

D I SpLT 

0007 

tapeck 

D0|0 

NWPUS 

00 11 

NIOI* 

00 1 z 

N 1 02S 

00J3 

NERR3S 


STQRASE ASSIGNMENT FOR VARIABLES [BLOCK, TyPE, RELATIVE LOCATION, NAME) 


0001 

000020 

1 256 

0001 


0000*13 

SOL 

0000 


000022 

1 lOF 

0000 


000003 

»15F 

00u3 

R 

000000 

A 

0003 

000132 

AyE 

0003 


ooouoo 

CTL 

0003 


0010*1 

OR 

0003 

R 

001*01 

H 

0003 

0 

000134 

MBAr 

0003 

000015 

ITAPE 

0003 

I 

000020 

|UN1T» 

cooo 

I 

000001 

I* 

0003 

R 

000034 

LABEL 

0000 

1 

□00002 

LSTftT 

0003 

R 000 INN 

P 

0003 

0 

000132 

pSAR 

0003 

R 

000011 

pbarJ 

0003 


000021 

PROP 

0003 

R 

0007*1 

R 

0003 

001041 

RATES 

0003 

0 

000 1 31 

hbar 

0003 

R 

000012 

RBARl 

0003 


001701 

RE 

0003 


005521 

RU 

0003 

0Q4311 

RV 

0003 


000(11 

state 

□ 000 

R 

oocooo 

time 

0003 

R 

002121 

u 

0003 

R 

003211 

V 


00101 

1 • 

SUBROUTINE OUTPUT ( T, OT . OTCR, REF 1 , REFZi 

REF3 f 

GNAT0793 

00103 

z» 

COHHQN /DATYV A ( 3700 1 


GNAT07V1 

00101 

3» 

REAL LABEL 


GNAT0795 

00105 

1» 

E8UIVALENCE 


q(Ii\T079* 

00105 

5« 

l (All!, CTL) 


GNaT 0777 

00105 

4* 

2. ( A ( 2 l 1 i PROP) 


GNAT07F0 

00105 

7* 

T , ( A ( 31 1 , LABEL 1 


GNAT07?? 

00 105 

8* 

1, UdU.AVE] 


GNAT 0800 

DO 1 05 

V» 

5, (Alia; I , STATE) 


GNaTOSOI 

00105 

IQ* 

4, < A ( 2 1 01) i RATES ) 


GNaTOSOZ 

0Q1Q4 

11* 

DIMENSION AVEI1DI, RATES t 1 400) , STATEI2000) 


GNaT0BO3 

00107 

1 2* 

OlHENSION P(20 ,20) ,R( 20 .20) .HI 20,20 1 ,u( 20,201 

,V (20,201 

GNaTOSqI 

00110 

1 3* 

DIMENSION CR(20.201 t RU12Q.20;. RV<20,201, 

RE (20 . 20 ) 

GNaTOSDS 

aa til 

It* 

EOUIVALENCE 


gnaToeq. 

£30 1 i i 

15* 

t (STATE ( n.Pl 


S’.aT0807 

001 1 1 

14* 

2, (STATEIIOn .R) 


gnatobo* 

00111 

17* 

3, (STATECBOn ,H) 


gnatosdv 



HU I 


004 1 1 

ia* 

H, ( STATE 1 1 201 1 *0 ) 


GNAT0B10 

OOIII 

IT* 

5, ( ST AT E 1 1 *0 | 1 tV) 


gnatosi l 

OOlll 

20* 

4, (BATES 1 1 ) .OR > 


gn*:d8|Z 

OOIII 

2r 

7, 1 RATES ( *0 1 1 , RE 1 


GN aT OS 13 

OOIII 

22* 

8, (RATES 18011 i RU 1 


GNaTOSIH 

OOIII 

23* 

T i (RATES! 1201 | . RV ) 


gnatosis 

00112 

2H* 

0 1 HENS 1 ON LABEL 1 1 0 1 


GNaTOB 1 4 

001 12 

25* 

equivalence 


GNATQB | 7 

00113 

24* 

H (LABElIH) tPBARI 1 


GNATOB I 8 

00113 

27* 

S, (LABEL (5) ,RBAr| 1 


GNATOB 1 9 

001 M 

23* 

DOUBLE PRECISION FBar, rbar. 

HSAR 

GNAT0820 

001 is 

2?* 

EQUIVALENCE IaVEUIiPBaRI* (AVEI3) ,RBAR> i ( AVE (5 ) ,HB A») 

GNAT0821 

0011* 

30* 

dimension itapeiai 


GNAT0822 

00117 

31* 

equivalence ictli ihi , itape . 

(I TAPE! HI ,1UNITR) 

GN AT0B23 

00117 

32* 

c 


GNAT0B2H 

00120 

33* 

time ■ T 


GNaT0B25 

00121 

3H* 

PBAri » pBAR / |HH. 


GNATQ826 

00122 

35* 

ROAR| * RBAR • 32*2 


GNAT0B27 

00123 

3** 

RRITE (6,915) LABEL 


GNAT0828 

00131 

37* 

9 1 5 FORMAT tlKlt 2QX, S^HGENERAL 

NUMERICAL ANALYSIS of TRANSPORT, 

Q NaT0829 

00131 

38* 

1 10X. IHMF j hej vhiller // 

lOE 1 2 1 7 // ) 

GNaTO830 

00132 

39* 

IF <IUNITR<LE<01 GO TO 50 


GNAT0831 

00I3H 

HO* 

call TAREPR (O.IUNITR) 


GNAT0832 

00 | 35 

HI* 

CALL TAPERR < 0 , ION I TR. 3700 , A 

<lstati 

GNATD833 

0013# 

H2* 

SO CONTINUE 


GNAT083H 

00137 

H3* 

CALL DISPLT 1 P<REF1<<0049HH 

1 

GN aT0835 

00 1 HO 

HH* 

RRITE <4,9101 


GNAT0834 

001H2 

HS* 

910 FORMAT 1 INI 1 


GNAT0B37 

001H3 

Hi* 

CALL OISPLT 1 R.REF2.32.2 1 


GNAT0838 

001HH 

H7* 

RRITE <4,9101 


GNaT0839 

001 H4 

H8* 

CALL OISPLT < M.REE3, 1. 1 


GNAT DSHD 

DO | H7 

HV* 

RRITE <4,9101 


GNATD8HI 

00151 

SO* 

CALL OISPLT ( U, O-i 1, ) 


GNaTDSHZ 

001S2 

51* 

KRlTE (4,910) 


GNAT0BH3 

0Q15H 

52* 

CALL DISPLT ( V, 0,, 1, ) 


GNaTOBHH 

00155 

S3* 

IF t IUNITR.GT.O) CALL TAPECK 

(LSTAT) 

GNaTOSHS 

00157 

5H* 

RETURN 


GNaTOSH* 

001*0 

55* 

END 


GNaT03H7 


END OF UN1V*C Jio# FORTRAN V COMPILATION* 0 'DIAGNOSTIC* HE55AGE15I 
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a FOR OIF'- Y t 015PLT 

UNIVAC IluB FORTRAN V LEVEL 220* 00 1 B FSOIBH 
This COMPILATION WAS DONE on 27 JAN 71 AT |m;33!S9 


GNAT08H8 


27 JAN 7J 


iNiaoisv.Mo*) 


SUBROUTINE DISPLY ENTRY POINT 000173 
STORAGE USED (BLOCK. NAPE i LENGTH) 


0001 

•CODE 

000207 

□ 000 

•data 

oao&*3 

0002 

•.lank 

ooooao 

0003 

data 

0071 *H 


EXTERNAL REFERENCES (BLOCK, NANEl 


DOOR 

NWDU* 

000& 

NIDI* 

000* 

NJ02* 

0007 

NERR3S 


VARIABLES (BLOCK, TyPE, RELATIVE LOCATION, NAME) 


STORAGE ASSIGNMENT FOR 

0001 00001* | ) 2C 

0001 OOOISq 1S7G 

0000 ] 000*23 I 
0003 000050 LIMITS 

0003 I 000120 NO 


0001 00003* U 76 

0000 000*25 VOOF 

0000 I 000*21 11 

0003 1 000121 NF 
Ouui « SOOOVO TIME 


0001 OOOOH7 127S 
0000 000*31 »01F 

0000 I 000*22 12 
0003 1 000050 NG 
0000 R 000000 T 


0001 0000*T 1355 

0003 R 000000 A 
0000 I 000*20 J 
0003 000122 NODES 


000 1 000123 1$0G 

0003 000000 Ctl 

0000 1 000*2*1 A 
0003 I 00007*1 NS 


ooir 

1 • 

SUBROUTINE DJSPLT ( X, XBAR t CONST 1 

GHAT08MV 

00103 

2* 

COMMON /DATA/ A 1 3700 I 

GHaTOBSO 

ooiom 

3* 

EQUIVALENCE 

GNAT0B51 

0010*1 

*t» 

1 1 A ( 1 ) ,CTL 1 

GliA T0B52 

OOlOR 

5* 

3, Dim, LIMITS) 

GHATDB53 

00105 

*• 

EQUIVALENCE (AI33).TIME) 

GIIAT08SM 

□ Q10* 

J« 

DIMENSION LtMITS(SO) , HG ( 20 ) , NS120) 

ghatosss 

00107 

B* 

EQUIVALENCE (LIMITS 1 11 ,NG1 , (LIMITS (21 1 ,NSI 

GHAT0B5* 

00107 

?• 

] , (LIMITS INI ) ,N0) > (LIMITS I m2) ,«F) . (LI HITS (M3) .NODES) 

GNAT0BS7 

001 10 

10* 

DIMENSION X (20 • 20 ) ■ r(20.20) 

GI1AT0B5B 

00110 

11* 

c 

GIUTOBSV 

ootil 

12- 

00 100 JaNO.NF 

GNaTOIAO 

00 111 

13« 

11 - NGIJ) 

GNaTOB* 1 

001 is 

IN* 

12 • NSIJI 

GNAT0**2 

0011* 

IS* 

DO 200 1 *1 I i 1 2 

GHAT0843 

00121 

It* 

T(I.J) » t XII. J) > XBAR 1 • CONST 

GNATOBtM 

00122 

17* 

200 CONTINUE 

GHaT0S*5 

00I2H 

!*• 

too CONTINUE 

GNaTOB** 

^0012* 

IV* 

DO 300 K*NQ . Nf 

GHaT0«*7 

00131 

20* 

J • NF *1 -K 

gnatoc*b 

00132 

21* 

300 WRITE [4.9011 3. 1 T 1 1 . J) , 1 «N(j , i o 1 

GNAT0B4V 

001 M2 

22* 

WRITE (t,VOO| TIME 

GNAT0B70 



no 


00115 

23* 


IF 1 NF.LE.10 1 RETURN 

GHAT0871 

00117 

21* 


00 100 K*NO iNF 

GN*T0872 

00152 

25* 


J • NF *1 *K 

SNAT0B73 

00 153 

2** 


12 * NS | J ) 

CNATD871 

00151 

27* 

100 

WRITE (*.1011 J. IT< I tJ) t l-l 1 . 12) 

GNATQ875 

001*1 

28* 

10O 

FORMAT l 1 H t inT l M£ *. FJ2.3 J 

SHAT087* 

001*5 

21* 

101 

FORMAT 1 l H 1 13. 3 A t 10E12.5 1 

SNAT0877 

0014* 

30* 


RETURN 

SNAT0878 

001*7 

31* 


END 

SNATD871 


end of 

ON 1 VAC 

1108 FORTRAN V COMPILATION. 0 •DIAGNOSTIC* MESSAGE (S ) 




B FOR CR2TAP.CR2TAP 

ONI VAC 5108 FORTRAN V LEVEL 220* G0|8 F501IH 

this compilation was done on 27 jan 71 at im:3s;di 


GnatObbQ 


27 JlN 71 


M 


subroutine cr 2 tap entry point 000112 


STORAGE USED taLCCK. name* LENGTHl 

0001 * CODE 000 1 2*1 

0000 *DATA 0000S2 

0002 ‘BLANK 000000 


External references (block, namei 


0003 

NREWS 

DOOM 

NKDU* 

OOOS 

NIOIS 

000* 

NIOZ* 

0007 

NRDUS 

0010 

NERR 3 s 


STORAGE 

ASSIGNMENT for 

variables 

{BLOCK, 

TYPE, 

RELATIVE 

LOCATION, 

i na m Ej 







0001 

ooooo* 

lOOL 

0001 

0000 1 N 

1 1 1 G 

0001 

00QD2R 

ll*G 

0001 

000038 

122* 

0001 

DODO 1 ** 

J30S 

0001 

0000*1 

137G 

oooo 

000020 

2000F 

0000 

000022 

2100F 

0000 

000026 

ROOQF 

ODD) 

qooqTh 

»ODL 

0000 I 

000000 

DATA 

0000 1 

00001* 

I 

OOQQ 

1 000017 

IS 






00101 



00101 

£s 

C 

00101 

3. 

C 

00101 

*(• 

c 

00101 

&• 

c 

ooiot 

*• 

c 

00101 

7* 

c 

00 1 01 

8* 

c 

DO 101 

9* 

c 

00101 

10* 

c 

00103 

11* 


001 OR 

U* 


00 1 OR 

13* 

c 

00105 

1 R* 


0010* 

15* 


0010* 

1** 

c 

0010* 

17* 

c 

00107 

IS* 


OBI is 

19* 


00115 

20- 

c 

00120 

21* 


00120 

22* 

c 

'-so 12* 

23* 



• •• 


100 


SUBROUTINE CR2TAP (OpTERH < 1 NDaTA • I MAGES) SNATOBII 

THIS PROGRAM rill GENERATE A BCD CaRO IHaGE TaPE FOR USE aS aN s«aT08»2 
ALTERNATE INPUT DATA FILE. INPUT IS EJECTED ON UNIT FIVE AND AN G«aT0883 
•OPTERM * IS EXPECTED in card COLUMNS l.i at end OF Each OPERATION. GNaTDIBH 
ARGUMENT TYPE EXAHPLE DESCRIPTION GNaTDMS 

OPTeRM INTEGER DEKENO OPERATION TERMINATOR, Six HOLERITm CHARgNaTDSBA 

INPATA INTEGER ? ALTERNATE data input file DESIGNATOR GNaTOS87 

1 «lE. INDATA. LE.2TIXCLUDE S*7|25,26) qHaT0888 
IMAGES INTEGER 3? NUMBER OF CARD IMAGES PER HEADING GKaTOBSV 

normally lines per page less tad gnaTORVO 

INTEGER DATA . OPTERM GNaT089I 

DIMENSION DATA IIRI GNATD892 


REMIND THE ALTERNATE INPUT OAT* FILE 
REWIND INDATA 
CONTINUE 

PAGE EJECT PRELIMINARY TO WRITING THE RECCRD5 IN THE CARD TO TAPE 


GHAT0893 

GHaT0*9H 

GNAT089S 

GNaT089A 


OPERATION PLUS COLUMN INDICATORS 
WRITE U.HOOOI iltl-I.BI 
DO SOD I * 1 , IMAGES 
READ AN INPUT DATA CARD 


GNAT0S97 

GHATDB98 

GNAT0S99 

GNAT0900 


READ (S, 20001 DATA 

white the image on the alternate data input file 
write ( INDATA, 2009) DATA 


GNAT09D I 

GNaT0902 

GN*TD90J 


!3*t; ]. t 
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00)24 

29* 

C •*• 

■RITE OUTPUT TAPE 



0013*1 

25* 


■rite 1 4,2)00) ) , Data 



00139 

24* 

c 

IF operation terminator, cease 

PROCESSES 

00193 

27 • 


IF IPATAU ),£9.0PTERKI 50 TO 

900 


00195 

26* 

500 

CONTINUE 



00)97 

29* 


CO TO 100 



00)50 

30* 

900 

CONTINUE 



00150 

3|* 

C 

REWIND THE ALTERNATE 

INPUT DATA 

FILE 

00)5) 

32* 


REWIND INDATA 



00)52 

33* 


RETURN 



00153 

39* 

2000 

FORHAT ( 1 3 A 6 , a 2 ) 



00159 

35* 

2100 

FORHAT l)5,2H.*. 1JA6.A2. |H*1 



00155 

34* 

9000 

FORHATUOHl CARP C0L,)7,7|10| 



00154 

37* 


END 




GHATQ9Q9 
GMATQ905 
GMaT0?04 
GHAT0907 
GMATQ90# 
GMATD909 
GNATQ910 
GNATQ9L1 
CNAT0912 
B«aTD91 J 
GNaTO’19 

SMATD915 
SNATD91* 
GNATD9 17 


END 

CR2TAP 

OF UNIVAc 1)08 FORTRAN V COMPILATION, 
symbolic 

0 *DI ACROSTIC* MESSAGE 1 S ) 

23 DEC 70 2310512) 

□ 

01 95)554 

)9 

100 

IDELETEP) 

CR2TAP 

COPE RELOCATABLE 

23 DEC 70 23;09:2] 

1 

Q144Q494 

29 

1 

(OELETEPl 




0 

0)640676 

19 

19 
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B fOR ROTAPE, HDtAPE 

UNlVAC lloe FORTRAN V LEVEL 2204 OOJe FSOlgH 

this compilation *as ©one on n jan 71 at i m t 3 4 ; 02 


SNAyOVlS 


subroutine 

ROTAPE 

ENTRY POINT 0000*5 

STORAGE used (block, 

name, length) 

0001 

• CODE 

000047 

0000 

• data 

000014 

0002 

•blank 

000000 

0003 

DATA 

007144 


EXTERNAL REFERENCES (BLOCK, NAME ) 

0004 TAPER* 

0005 TAPEPS 
0004 TAPEPR 
0007 TAPER0 

0010 tapeck 

0011 NERR3S 


SToRASE ASSIGNMENT FOR VARIABLES {BLOCK, TyPE, RELATIVE LOCATION, NAME) 

0000 R 000000 A 0003 000000 CTL 0003 I 00001* (FILER 0003 

0003 | ODOQIt IUNITR 0000 ] 000001 JFJLER 0000 [ 000002 JKECR 0000 


00101 1* SUBROUTINE ROTAPE 

00103 2* CDNHON /DATA/ A(37QD] 

D0IQ4 3* DIMENSION CTLI20), I TAME (41 

00106 4* EfiUIvALENCE UUl.CTt,). (CTL ( 1 4) , I TAPE J 

00105 6* 1, ( ITAPE ( 1 ) t JUN I TR 1 , 1 1 TAPE 12 > , IF iLERl . UTAPEI3) ,IR£CRl 

00104 4* JUNITR * I UNI TR 

00107 7 • JF1LER - IFtLER 

00110 S« JHECR ■ 1RECR 

□Dill V* (ALL TAPER* (0, JUNITR) 

00112 10* call TAPEPS (0, JUNITR, JF1LER ,JR£CR*ll 

00113 II« CALL TAPEPR (O.JUNITrj 

00114 12* CALL TAPEHO (0 ,JUN 1 1 R ,3700 i A ,L5TAT 1 

001 Is 13* CALL TAPECK (LSTAT) 

0011* 14* CALL TAPER* IQ.JUNtTRi 

OD 11 7 |S* I UN | TR ■ JUNITR 

00120 I *• IF 1 LEM * JFILER 

00121 17» 1 RECR • JRECR 

00122 10* RETURN 

00123 IV* END 


END OF UNIVAC 1|0* FORTRAN V COMPILATION 


U ^DIAGNOSTIC* MESSAGE (51 


27 JAN 71 


14 *34 • 2*3*5 


000017 IRECR 00d3 00CQI5 ITAPE 

000000 JUNITR 0000 1 000003 LST a T 


GNaTOVI? 
GNAT0V20 
GNATQV2 l 
GNaT0?2Z 
GNAT0T23 
CNAT0T24 
SNATDV25 
5NAT0V24 
GNATDV27 
GNaTOVZS 
GNAT0V2? 
GNATDV30 
GNATOV 3 1 
GNAT0V32 
GNAT0733 
GNAT0434 
GNAT073S 
GNATQV34 
GNaT0»37 



8 FOR TaPEIO.TAPEJO 

ON l VaC 1)08 FORTRAN V LEVEL 22D6 0018 FSOISh 
This COMPILATION WAS done ON 27 JAN 71 at 1 H » 3R I 03 


GNATQV38 


SUgROUT INE TAPEPS ENTRY POINT 0003*2 

TAPERD ENTRT POINT OOORR I 

TY.PEWR ENTRY POINT 000527 

TAPERW ENTRT POINT 00311s 

TApEOF ENTRT POINT 000*43 

TAPEPR ENTRY POINT 000731 


tapeck entry point 000777 
STORAGE USED (BLOCK, NAME, LENGTH | 

0001 • CODE LO 1 005 

OOOD *DAT A 0001*12 

0002 'BLANK 000000 


EXTERNAL REFERENCES I BLOCK , NAHEJ 

0003 NTRAN 
OOOH GOFIL 
P 0006 QQREC 

J>- 0004 KILLER 

0007 NREWS 
00(0 NRBU* 

0011 N 1 0 1 S 

0012 N I 3 2 * 

0013 NkBUS 
00 l R NWEF* 

0015 NWOUS 

0014 NERR3S 


STORAGE 

assignment FOR 

VARIABLES 

(BLOCK, 

t t pe, 

RELATIVE 

LOCATION, 

i NAME) 


000] 

□00025 

10L 

0001 

000 325 

IDOL 

0001 

000253 

1 2qL 

0001 

000) 

000053 

ZQL 

0000 

OGOiOR 

200F 

0001 

000175 

22] G 

0001 

0001 

0000*3 

35L 

0001 

□ GOD**) 

ROL 

0001 

000D7R 

RSL 

0031 

0001 

00033*1 

7QL 

0001 

0001 MO 

SQL 

0001 

00033* 

“OL 

0000 

0009 1 

000075 

i* 

0000 1 

000000 

JFILE 

0000 

] 000035 

JREC 

□ 000 

0000 l 

000 I 0 i 

KREC 

0000 ] 

□00100 

MSTAT 

0000 

] 000 102 

NSTaT 

OJOQ 


00 100 I* C 

00100 2* C ITYPE • TTPE OF DATA TAPE 


27 JAN 7 


ir;or; 3 . sis 


000222 l*OL COO’- QGQ273 l»Ot 

000232 2R3G 001)1 005ER3 JOL 

000145 Sn L 0001 0U012R 55 t 

00JJD3 I 0000 J 00007*) j* 

000377 JSR 30 JO J 00037* SFJuE 

□00072 TTp£ 


5NAT0V39 

GMATOWO 
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00109 

3* 

C 


o « ntRan 


GNAT09H t 

00100 

*(• 

c 


1 • FORTRAN BINARY 


GNAT0912 

00100 

5* 

c 


1 UN IT - PHYSICAL UNIT FOR TAPE ASSIGNMENT 


GN*To9t3 

00100 

5* 

c 


IFlLE • FILES TO BE SKIPPED 


GN AT09tt 

00100 

7* 

c 


1REC - RECORDS TO BE SKIPPED 


GN A T09t5 

00100 

e* 

c 


I PORDS - NUMBER OF DATA WORDS 


GNATOtti 

00100 

?o 

c 


A - STORAGE area FOR DATA RONDS 


GN aTD9H7 

00100 

l 

c 


LSTaT - STATUS WORD FOR NTRAN READ/1RJTE 


GNAT099B 

00100 

11* 

c 


-1 . TRANSMISSION NOT COMPLETE 


GN*T09t9 

00100 

12* 

c 


-2 * eRd"°F”fILe for Re Ad * eho”Of”TAPe or jjRun-i 

FlLc 

t NAT0950 

00100 

13* 

c 


FOR WRITE 


4NAT0951 

00100 

It* 

c 


•3 . DEVICE ERROR 


GNAT0952 

00100 

15* 

c 


-1 ■ transmission aborted 


GN a T0953 

00100 

15* 

c 


H ■ NUMBER OF WOP.OS (IWOROSI TRANSMITTED WHEN 


GNAT095M 

00100 

17* 

c 


TRANSMISSION Is COMPLETE 


GNAT0955 

00100 

IB* 

c 


LSTAT • STATUS WORD FOR FORTRAN REAO/RRITE 


GNAT0954 

00100 

It* 

c 


N . NUMBER OF WORDS (1W0RDSI TRANSMITTED WHEN 


GN*T0bs7 

00100 

20* 

c 


transmission is complete 


GNAT095B 

00100 

21* 

c 




GNAT0959 

00100 

22* 

c 


msTat - status boro for fortran file skipping 


GNaT0940 

00100 

23* 

c 


0 * NORMAL COMPLETION AS REQUESTED 


GN*T094I 

ootoo 

2t* 

c 


1 * INAPPROPRIATE FORTRAN UNIT 


GNAT0942 

00 100 

25* 

c 


2 > NOT USED 


GNAT0943 

00100 

25* 

c 


3 * UNEXPECTED TERMINATION BY A PERMANENT READ 

ERROR 

GH*T0 ? *t 

□0100 

27* 

c 


NSTaT * status word for fortran record skipping 


GNAT0955 

0Q100 

28* 

c 


0 * NORMAL COMPLETION AS REQUESTED 


GNaT0944 

00100 

29* 

c 


1 * ihappropriate fortran UNIT 


GHAT 09 47 

OQIOO 

30* 

c 


2 - UNEXPECTED TERMINATION BT SKIPPING tHR»U«H 

AN 

GNATD94B 

ooioo 

31* 

c 


end-ofafile Or reaching THE beginnins of Tape 

GN*T0949 

Ooloo 

32* 

c 


3 . UNEXPECTED TERMINATION BT A PERMANENT READ 

ERROR 

GNaT0970 

00109 

33* 

c 




GN A T097I 

00101 

3t« 



SUBROUTINE TaPEPS UTTPEilONIT.lFlLEilRECl 


GN4T0772 

0010} 

35* 



0IMEN510N UFRE1291, JREC (29 ) 


GNAT0973 

00 1 OH 

35* 



INTEGER TTPEC2) 


GNAT0979 

00105 

37* 



DATA JFIlEi JSEC, / 29*0i 29*0 / 


GNAT3975 

00110 

38* 



DATA TTPE / 4m NTRA.5HF0RTRA / 


GNAT0974 

00112 

39* 



DATA 14 / * / 


GNAT0977 

00112 

to* 

c 




GUAT097B 

00112 

tl* 

c 


THIS ROUTINE "ILL POSITION EITHER AN NTRAN OR A FORTRAN 

GN4T0979 

00112 

H2* 

c 


BINARY TAFE TO The proper FILE AND/OR record* 


GNaT0980 

00112 

t3* 

c 




GNAT098I 

OQtlt 

1 1* 



IF IIFILE iEG > 0) GO TO 95 


GN*T09B2 

00115 

■IS* 



KFILE * IFILE 


GNATD9B3 

00117 

t5* 


5 

CONTINUE 


GNAT0989 

00120 

H7 • 



JSR * JFILEI IUNJTI * KFILE 


GNaT09S5 

00120 

tB* 

c 


TEST FOR TTPE Of TAPE 


GNATQ9B4 

00121 

1 9* 



IF UTTPE .EQ. 11 GO TO 30 


GNXT09B7 

00 121 

50* 

c 


NTRAN TAPE 


GNAT0988 

00123 

51* 



IF (JSR .GT, 0) GO TO 10 


GN A T09B9 

00125 

52* 



CALL NTRAN ( 1 UNIT • 1 0 J 


GNATQ990 

00125 

53* 



GO TO 35 


GNAT099I 

00127 

5t* 


10 

CONTINUE 


GN4T099Z 

00130 

55* 



CALL NTRAN (IUNITiBiKFILEM) 


GNAT0993 

00131 

54* 



CALL NTRAN (lUNlTiBtll 


GNAT099H 

-00132 

57* 



GO TO 90 


CNAT0995 

00133 

58 * 


30 

CONTINUE 


GNAT099A 

00133 

59* 

c 


BINARY TAPE 


GNAT0997 

0013* 

40* 



IF (JSR .ST. 01 CD TO 20 


GNAT099B 



13* 
137 
HO 
111 
IH2 
i HI 
llH5 
IIS* 
I l S7 
ISO 
IS! 

□ 0 1 52 
OOJ5S 

□ 015S 
001S3 
03 1 57 

□ OHO 
001*2 
001*2 
00]*3 

□ 0 1*3 
001*5 
001** 
001*7 
001*7 

00170 

00171 
00173 
00175 
0017* 
00177 

H 00200 

Ox Q 0201 
00201 
0O20S 
0020* 
0020* 
0020* 
0020 * 
0020* 
0020* 
00207 
00210 
00213 
0021 3 
0D2H 

00216 
0021 * 
002 |* 

00217 
00225 
0QZ2* 
0022 * 
00227 

-00227 

UJZ27 

00227 

00227 


*H 
* 2 * 
*3* 
*S« 
*5» 
**• 
*7* 
*a» 
*?« 
70* 
71. 
72* 
73* 
7S* 
75* 
7*. 
77 • 
78« 
79* 
80* 
S|» 
B2* 
E3= 
as* 
as* 
86* 
87* 
88 * 
8 ** 
fO* 
* 1 * 
*2* 
?3» 
?S* 
75* 
94* 
77* 
78. 
77" 
100 * 
let* 
102 . 
103* 
IOS* 
105* 
10 *. 
107* 
108* 
107. 
1 10 * 
111* 
112* 
113* 
1 1 S* 
115* 
1 I** 
117. 
1 I 8* 


remind (unit 

GO TO 35 
20 CONTINUE 

CALL WSFIL IIUN|T»KF ICC i NSTAT I 
|F 1HSTAT HE. 01 GO TO 70 
GO TO SO 
3S CONTINUE 
JSR « 0 
So CONTINUE 

JF1LEUUN1T1 - JSB 

JREC 1 1UN I T I ■ 0 

IF HREC .LT. 01 GO TO 70 

ss continue 

C POSITION tape to proper recorO 

IF (IREC .EO. 01 GO TO ¥D 
KREC > JRECUUNIT; . IREC 
IF UREC .LT. 01 GO TO 80 
JRECIIUNITI * KREC 
C TEST FOR TYPE OF T*PE 

IF lITTFE .EG. 11 GO TO 55 

C NTRAN TAPE 

CALL NTRAN 1 1 UNIT .7 . f REC I 
GO TO 70 
55 CONTINUE 
C BINARY TAPE 

CALL GOREC liUNiT.IREC.NSTAT) 

IF 1NST*T .Ed. 21 GO TO 70 

IF 1NSTAT HE, 01 GO T 0 ?□ 

GO TO 70 
80 CONTINUE 

JRECUUNIT) ■ 0 

kfile * *1 

IF IITYFEI 20,10.20 
C 

ENTRY TaPERO IITYPE.IUNIT.HORDS.a.LSTaTI 
INTEGER All) 

C 

C This ROUTINE MILL READ FROM EITHER A FORTRAN BINARY TaPE OR 

C an NTRaN TaPE and STORE I*ORdS AHOUNT of INFORMATION INTO 

C THE AP-RAY a, 

c 

JREttlUNlTl - JRECUUNIT] * I 
IF iJTYPEl 50,40,50 
*0 CONTINUE 

c read ntran tape and store data into a 

call NTRaN (IUNIT.Z.HOKoS.A.lSTaT) 

GO TO 70 
50 CONTINUE 

C READ BINARY TAPE AND STORE OaTa INTO A 

READ UUN1T1 lAi! l.l-l.HORDS) 

LSTAT • IWOROS 
GO TO 70 
C 

ENTRY TAPEAR IITTPE.IUN1T,IR0RDS,a.l5TaT| 

C 

c Th 1 5 ROUTINE * ILL *RjT£ EITHER A BLOC* OF INFORMATION HORqS 

c long on an ntran tape or a record of ifcrds on a fortran 

C BINARY tape 


GNAT0977 
GNAT 1000 
gnatiooi 
GN aT 1 002 
GNaT 1003 
GHaTIOOS 
GMaTIOOS 
GNaTIOD* 
GN*.TI007 
GNaTioOS 
GNaT 1007 
gnaT ioio 
GNATIOI 1 
GN AT J 0 12 
GNAT 1013 
GNAT1QH 
GNaTIOIS 
GNATIOI* 
GNAT 1017 
GNAT 1018 
GNaT 1017 
GNaT 1020 
GNAT L 02 I 
GNaT 1022 
S NaT13Z3 
GNaI 102S 
GNAT 1 025 
GNAT 1 02 * 
GNAT 1021 
6NaT|02, 
GNaT 1027 
GNAT 1 030 
GNAT 1031 
GNaTI 032 
GNaTI 033 
GNaTI C3S 
GNATI 035 
GN*TI03* 
GNaT 1037 
GNAT 1038 
GNAT |03? 
GNitlOsO 
GNaTIOSI 
GNAT1DS2 
GN.T 10S3 
gnaTIOss 
gnat 10S5 
GNaTIQS* 
GNAT10S7 
GN*T I0S8 
GNAT IOS? 
GNAT 1 05 0 
GNaTIOSI 
GNaT|052 
gnat 1053 
gnaUQSS 
gnaTioss 
GNaTIOS* 


00227 

1 19* 

C 



GNAT 1 057 

0O231 

120* 



jREc < 1 1 N I T 1 ■ JREC l ] UN ( T 1 * l 

gnatjo&b 

00231 

121* 

C 


TEST FOR TTFE OF TAPE 

GNAT 1 059 

00232 

122* 



IF IIYTPE1 l*0,MO,l*a 

GNAT 10*0 

00235 

123* 


I 40 

CONTINUE 

C*iATU*l 

00235 

124* 

c 


WRITE an ntran block 

GNaT,0*2 

0023* 

125* 



CALL NTRAN 1IHNIT,1,1«OROS,A,L5T*TI 

GNaTI0*3 

00237 

12** 



SO TO 40 

GNaT 10*4 

002*10 

127* 


1*0 

CONTINUE 

GNaT10*5 

002*10 

US* 

c 


write a fortran record 

GN A T|0*t 

0024) 

124* 



■HITe CZUNlT) IA<1),1»1,1«0 RdS! 

GN a T10*7 

00247 

130* 



LSTaT ■ J WORDS 

GNaU0*8 

QC2S0 

(31* 



CO TO 40 

6 NaTIU*4 

00250 

132* 

c 



GNAT 1070 

00251 

113* 



ENTRY TAPER* 1 I T YPE , I0N1 T J 

GNaT107I 

0025) 

134* 

c 



GNaT 1 072 

00251 

135* 

c 


This ROUTINE WILL REWIND EITHER AN NTRaN OR A FORTRAN 

GNaTI 073 

00251 

13** 

c 


binary tape. 

GNaT 1074 

00251 

137* 

c 



GNaT |075 

00253 

138* 



JFlLEllUNJTI » 0 

SNaT107* 

00254 

134* 



JREC IIUNITI * 0 

CHAT 1077 

00254 

140* 

c 


test for type OF tape 

GNAT 1078 

00255 

Ml* 



if (JTYPEI 120,110,120 

GNaT 1074 

002*0 

142* 


ila 

CONTINUE 

GNaTIOSO 

002*0 

143* 

c 


rewind The ntran tape 

GNaT 108 1 

002*1 

144* 



CALL NTRAN tIUNIT.IOI 

GNaTI 082 

002*2 

MB* 



CO TO 40 

GNAT 1083 

002*3 

M** 


Uo 

CONTINUE 

GNAT10B4 

002*3 

147* 

c 


REWINO FORTRAN BINARY TAPE 

GNaT 1085 

002*4 

MS* 



REWIND 1 UNIT 

GNaTIOB* 

002*5 

144* 



GO TO 40 

GNaT l 087 

00245 

ISO* 

c 



GNaTI 088 

002** 

151* 



ENTRY TAPEDF ilTYPE.lUNITl 

GNaT 1084 

002** 

152* 

c 



GNaTIOVO 

002** 

IS3* 

c 


this routine will write an end of file ha«k on eiTheR ah 

GNaT 1041 

002** 

154* 

c 


NTRAN tape or A FORTRAN BINART TaPE, 

SNaTIOTZ 

002** 

ISS* 

c 



GNaT 1043 

00270 

IS** 



JRECtlUNlT) * 0 

GNATI044 

00271 

157* 



JFILEIIUNITJ » JFILEllUNITl ♦ 1 

GNATIQ45 

00271 

!5»* 

c 


TEST FOR TYPE OF TAPE 

GNAT 104* 

00272 

154* 



IF t I TYPE | 140, ISO, MW 

GNAT 1097 

00275 

1*0* 


1*0 

CONTINUE 

GNAT1Q4B 

00275 

1*1* 

c 


WRITE AN EOF NARK OH THE NTRAN TaPE 

GNaTIO?? 

0027* 

1*2* 



CALL NTRAN ( 1 UN I T , 4 3 

GNaTI 100 

00277 

1*3* 



GO TO 40 

GNATI 101 

00300 

1*4* 


l"a 

CONTINUE 

GNATI 102 

00300 

1*5* 

c 


WRITE AN EOF HARK ON THE FORTRAN BINARY TAPE 

GNaTI 103 

00301 

!**• 



END FILE I UNI T 

GNaTI 104 

00302 

1*7* 



GO TO 40 

GNaT] 105 

00302 

H9=> 

r 



GNaTIIO* 

00303 

1*4* 



ENTRY TAPEPR 1 1 TYPE, IUN! T 1 

GNaTI 107 

00303 

170* 

c 



GNAT l 1 08 

□ 0303 

171* 

c 


THIS ROUTINE WILL PRINT ThE STATUS OF THE UNIT JUNlT, 

GNaT 1 1 09 

00303 

172«> 

c 



gnat i i i o 

00305 

173* 



KFILE ■ JFILEIIUNITJ * i 

snath I 1 

OQ30* 

174* 



KREC * JREC IIUNITI * l 

GNAT 1112 

00307 

175* 



WRITE u*, 2001 TYPE! ITYPE*1> , IUN1T, KFILE, KREC 

GNaT 1113 

-qoois 

17** 



50 TO 40 

GNATI 114 


118 


003(5 

177* 

c 





GNAT 1115 

□ OJJ* 

174* 



ENTRT TAPECX IJSTATl 



GNAT11I4 

00314 

t 79* 

c 





GNaT 1 1(7 

003(6 

lao* 

c 


THIS ROUTINE BILL TEST THE 

STATUS OF JSTAt FOR NTRAh 1/0 

gnaTH i» 

003(6 

1BJ. 

c 





GNATU19 

00320 

182* 


(00 

If IJSTAT .Efl. «|1 00 TO 100 



SNATU20 

00322 

(B3* 



(F USTATI 70,70.90 



SNaT (121 

00322 

t 0*t« 

c 


error has occurred 



GNAT1 122 

00325 

1 as* 


TO 

CONTINUE 



GNAT 1123 

00324 

IN* 



CALL KILLER 



GNAT1129 

00326 

187* 

c 





sNaTUZS 

00322 

188* 


»0 

CONTINUE 



GNATI124 

00330 

|R9* 



RETURN 



GNAT 1127 

00330 

170* 

c 





GNAT11Z4 

Oo33( 

191* 


200 

FORHaT ( 1 X a6 ,6hN UNI T 13 i 1 1 H IS 

ON 

F1LE13,7h SECORBIN) 

GNaTI 129 

00332 

192* 



END 



GNAT 1130 


end or 

D N 1 V A C 

1108 FORTRAN V COMPILATION, 

0 

‘DIAGNOSTIC* MESS AsE t S ] 




119 


a ELT 

WQPGS/VeR2,3 i 7Q03J2 | t752H 

4NAT 1131 

27 

JAR 

71 

l<t;3i : 

*,{07 

8 ElT 

tr*cer/cooe ,3 ,&aoi |5i *i9R5 


2 7 

JAN 

7 1 

it! j*i: 




120 


B XOT MAIN 


GHATI Jbl 


27 JaN 


li;31; *.S“3 


STARTING ADDRESS 011000 

CORE L IN ITS Ql ROOD 031371 IOOOOD 137523 1*3772 1*3777 


MAIN /CODE 

0 IDODPQ-IO013* 

1 0HQQ0-01521* 

NSTOPs/CODE 
I 015217-0152*1 

NlfcR* /CODE 

0 100137-10013'' 

1 0 152*5-015571 

2 I SOI “(0“ 1 0023 1 

NFKT* /CODE 

1 035572-01*157 

2 100232- 100215 

NCNyT*/CODE 
| 01*1*0-01**47 

2 1 0D2AA- 1 00331 


NFTV s /CODE 
1 01**70-01*712 

NOT) NS/COOE 

1 01*713-01730* 

2 100332-100371 

FpACKI/CODE 
l 017307-017352 

DEPTH /•••••• 

0 100375-100102 

NERRS /CODE 

0 150103-100512 

1 017353-017755 


N 1 0 I N*/C DDE 

1 0| 7754-020023 

2 100513-100573 
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NOUT* /CODE 

0 100571-100574 

1 02OD21-D2O471 
Z 100577-10041** 

NT»B* /CODE 
0 1 006 1 S- I 09713 


NBDCVl/cOOE 
0 10071*1-101130 


NLInPS/CODE 

0 101131-101137 

1 020475-022243 

2 101 |10-1Q|312 

H 1 N 1 NS/COOE 

1 022241-022135 

2 101313-101313 


OUTpllT/CODE 
□ 101311*101112 

I 022134-D22407 


D | SpLY /CODE 

0 101113-102275 

1 0224 | 0 - 0230|4 

BUUK /CC3E 

0 102274-102313 

1 023017-023214 

TEHP /CDOE 

0 102 . 311-102343 

1 023217-023340 

HEAtER/COOE 

0 102341-102373 
t 023341-023102 

0 JFF 2 /CODE 
D 102371-102125 
| 023103-023477 

ME*J*a /CODE 
D 103124-102170 

1 023700-021214 

BE OUT /CODE 

0 102171-102507 

1 021217-021370 

[•Iff /CODE 

0 102510-102511 

1 021371-025241 

V «EAN* /COOE 
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D l 025*1 2* 1 02* 1 0 
l 0252*5-025*5/ 


NQD»AU/-»*-»* 

0 102*11-113302 

RTPRES/COOE 

0 1 1 3303° 1 1 3*1 1 5 

1 0254*0-02*075 

e*p /code 

1 02*076-02*1** 

2 I 13H |*-| 13*135 

NEXp**/CODE 

t 026 j*7-02*2of 

2 I 13M3*-| 13*13* 


KXP AF*/ CODE 
t 02*205-02*3*3 
2 I 1 3H37 - 1 | 3H S3 

HXPAXS/CCOE 
| 02*3**1-02**10* 

2 I I3l*t*t-1 J3*f*l l t 

*L06 /root 

1 02*107-02*17*1 

2 I131l5-]1352| 


therm /code 

0 113522-113*12 

1 02**175-027517 

BJNSER/COCE 

0 113*13-113*35 

1 027520-027711 

TPCB /•«•••• 

0 113*34-123021 

SETA /CODE 

0 123022-123037 

1 027712-030017 

TEFCTR/-»»-.» 

0 1230*10-123133 

TAPE 1 0/CODE 

0 1231 3*1 -1 23275 

1 030020-03l02*t 

NFOUTX/COOE 

1 031025-0312*1* 

2 12327**123277 

nbuffs/cooe 
1 0312*17-031270 



2 | 23300- 1 2* 3 1 0 


NFJnP*/CODE 

1 03 1 27 1-03 1530 

2 12*311*12*312 

nr«nd*/cooe 
t 03 J 53 1 *03 1 623 


TRAcER/COOE 

0 I 2*3 | 3* 12*3)* 

1 03162**031635 

OaPoS / VER2 
) 03 ] 636*031 750 

2 1 2*3 1 5» | 2 *32 1 


NT**H /CODE 

0 12*322*1 2*322 

1 031751*033251 

2 12*323-12**52 

TU*BLS/COOE 
0 12**53-12**62 


TSCRHS/COOE 
0 12**63-12*535 


12*653 


T 5 »AP*/C 0 DE 
0 | 2 * 6 S*- 1 2 * 75 * 

TJNTLS/COOE 
0 12*755-125053 

ROTAPE/CDDE 
0 12505*. 125071 

| 033252-0333*0 

CR2T*P/C0DE 

0 125072-1251*3 

1 0333*1-033*6* 

NINPTS/CODE 

0 I 25 | **- | 2S l *5 

1 033*65*03*37* 

2 1251*6-125200 

OSTATE/****** 

0 125201*130300 


TRNSHT/*** • • • 

0 130301-130305 

C5TS /••••*• 

0 130306-13Q313 



124 


02 /•••••• 

0 UQ3H-J30323 

01 

0 I 3032H » | 3Q337 

PAT* 

0 1303H0-137523 


EKO Of ALLOCATION 1103 003?* 09099 



125 


CaRC_cOL . _ ..._1 z l 


1 .* 
z.* 

3 . • 

H.t 

S.‘ 

* aP°i-Lo 

* heater 

* scale 

* gy * - 

* 

• 

bottle CROSS 
NODES - 1 = 

FACTOR IS 2<*00 
2. E-a 

SECTION 
12. J * 

7.- 

5 input 



_ Si * 

PROP < *» 1 

* 1.2E-2 


9.» 

PROP < 8 1 « 

850* i 200. 


lO»* 

CTi-1 12i 

■ - 2 - E-a 


11. » 

CTLl 131 

= 2800. 


1 2 . * 

f Tl 15 t ■ 

I • H 


13. ♦ 

CTL(6I * 

b 2d * 6 


I 1 *.*. 

CTL12) * 

,5E~8| •£» -025 

IS, • 

CTL13I ■ 

.9 


1 6 * • 

CTLlHl * 

.QlZS 


|7,» 

ITAPECa) 

• lit.! 


1£- A 

CXJ.T 201 . 

5 2. _ 



!?.* *END 

nTeAH unit 1 IS DU fl L £ I record 


« 5 

- 20 NOoE DIAMETER 

10 


1 



GE«E«f t _A?jA l yS i S Of T»»iiS_PQHI p j K Ela*.|.L t Ef( 



•OOOQUOO 

•saoDoao-as .Qaoaoao 

•8500000+03 .*595085 

+07 .QOQOOUO 

• ooooeoo 

•QnoooaO 

•0000000 

•OODPQCO 




NTrAN 

unit 1 is 

On FILE 1 RECORD 1 











2fx 

- nfifinn 

. inama 

•QQEDD 

• nSQOn 

-fOPQCQ 

-i-DODDC 

•GDOOp 

•1*358+02 

. 1 *358+02 

. 1635B + Q2 



19 

•ooooo 

• oo odd 

* ODOCO 

• ooaon 

.00000 

• 1 *3*2 + 02 

, 1 63*2 + 02 

. 1 6362 + 02 

•1*3*2+02 

.163*2+02 



4* 

•00000 

•ooqoo 

•OOOOn 

• 1*3*7+02 

. 1*367+02 

. ,*3*7+02 

, 1 *367+02 

•16367+02 

. | *3*7 + 02 

, 1*3*7+02 



17 

•oodoo 

•00000 

• 1 *37(1+02 

• l*37(i + o2 

• 1*37(1 + 02 

.,*370+02 

• 1 *37q + o2 

• l*37o + o2 

• 1 *2-70+ o2 

•l*37 Q + 0 2 



1* 

-00000 

•oooon 

*1*375+02 

• 1 *375 + 02 

•1*375+02 

• , *375*02 

. 16375 + 02 

• 1 6375 + 02 

• 1*375 + 02 

• 1*37s+(j2 



15 

•00000 

• 1*379 + 02 

• 1*379 + 02 

• 1*379 + 02 

. 1*379+02 

. ; *37 9 + 02 

. 1 *379+02 

. 1 *379+02 

. 1*379+02 

. 1 6379+02 



L9 

• nnnnn 

• X*3JS3+02 

•4*38i+02 

_ »i*3 8.i±flz. 

i4*3 8j*H2 

11*381*1H_ . 

• ! *383 + 02 

• 1 *383 + 02 

• i *383+02 

. 1 *3Bl+n2 



13 

, j &3 8 7 + 0 2 

• 1 *38 7 + Q2 

• 1*387 + 02 

• 1 *387+02 

. 1*387+02 

. t 6387*02 

, 16387 + 02 

• 1 *387 + 02 

• 1 *387 + 02 

•1*387+02 



12 

■ I 639a*02 

• I*39o+tj2 

1 1 639(j+02 

• 1*39q+02 

• 1*390+02 

. 639 Q *02 

• 1 *39q + o2 

• 1*390 + 02 

, 1*39o + Q2 

• l*39 0 + 02 



1 1 

• 1 6395 + 02 

• 1*395 + 02 

* 1 *395*02 

• 1 *395 + 02 

■ 1*395 + 02 

• | *395*02 

• 1*395 + 02 

• 1 *395 + 02 

• 1*395 + 02 

* 1*395+02 



10 

• 16399*02 

t 1*399+02 

• | *399+02 

• | 6399+02 

•1*399+02 

• , *399*02 

■ 1*399 + 02 

• 1*399+02 

• 1*399 + 02 

• 16399 + Q2 



9 

• 1* 9 + 02 

• 1 *909 + 02 

• 1*909+02 

. 1*909 + 02 

. 1*909+02 

• t *908 + 02 

• 1*909 + 02 

. 1 *909*02 

• 1 *909 + 02 

. 1*909*02 



K— 

- i*-'* ft 7+0? 

• [ h *J p7 + ft? 

• i *9f>7 + n7 

. 1 Agn7 + n? 

. | *9q7 + Q 2 

• -j 6 -9-0.7 *02 

. 1 6907*02 

. 1*907+02 

. 1 *907+02 

• 1 1 907 + 02 



7 

•OODDD 

•3*91 ! +02 

• 1*9 1 1+02 

•1*91 1+02 

•1*91 1+02 

•1*81 1+02 

.16911+02 

• 1*811+02 

• 1*91 1 + 02 

• 1 *9 1 1 + 02 


* — 

A> 

-QOOOO 

• l-*9 t 2 

• | *91 A+02 

•1*916+02 

. 1 6916+02 

. , 681**02 

.1*916+02 

•1*914+02 

.1*916+02 

• 1*9J*+q2 



s 

•00000 

•00000 

• 1 *920 + 02 

• 1 *92o + 02 

. 1*920*02 

, , *82q + 02 

•1*92 0 +02 

• 1*920 + 02 

. 1*920+02 

•1*920+02 



9.-. 

•aoooo 

.00000 

• 1*823*02 

• 1 6923 + 02 

.16923+02 

. 1 *923*02 

. 1 *923 + 02 

• 1*923 + 02 

•16923+02 

• 1*923+02 



3 

•00000 

•00000 

*00000 

• 1 *928+02 

,1*928+02 

. | *928+02 

. 1*928+02 

•1*928+02 

. | *928+02 

• 1*928+02 



2 _ 

*00000- 

.cmgns 

• canon 

•nnnnn _ 

• OQQ-OO 

• j 3 ? *02 

• 1 *937 + 1,2 

• 1 *937 + 02 

• 1 *932 + 02 

« l*932*n2 



1 

•00000 

•00000 

-oaoao 

• 00000 

• 00000 

•OQOOn 

•OOOOO 

. 1*936+02 

.1*93**02 

• 1*93**02 


— 

4 1 we- *• 

-oaa 












2a 

. |63SB*02 

• I *308 + 02 

• | *356+02 










49 

* 16X62*02 

• 1*3*2 + 02 

* 1*3*2 + 02 

- 1*3*2 + 02 

. 1*362+02 








10 

• 1 6367*02 

• 1 *3*7 + 02 

• 1*3*7 + 02 

• 1*367 + 02 

• 1*3*7 + 02 

•163*7+02 

. 1*367 + 02 






1-7 

.-*-1.±Z7q.+ qZ.. 

- M *-47o + Q2 - 

m i*37n+n-2 

- 

. 1 *37n+n? 

, , 6 3 7(i -+0-2 

• 1 637 n*n2 

j-I*37fl*Q2 





Ifc 

. 16375*02 

• 1*375 + 02 

• 1*375 + 02 

• t *375+02 

•1*375+02 

■ ,6375*02 

. 16375+02 

•1*375+02 




- 

is 

. 1 6379+02 

• | *379+02 

i i *379 + 02 

• 1*379 + 02 

.1*379+02 

. ,437o*02 

.16379+02 

•1*379+02 

*i*i 7 *+oi 




19 

• 16383 + 02 

• 1*383 + 02 

• 1*383 + 02 

• 1*383' 

.1*383+02 

. 1 6383*02 

. 163B3 + 02 

• 1*383 + 02 

. 1*383+02 



— - - 

-43 

•• j 6387*02 

- 1*387 + 02 

• 1 *387+02 

• 1 *387* 

•1*387+02 

. 1*387*02 

.1*387+02 

1 1*387 + 02 

• 1 *387*02 

.1*387+02 



12 

. 16390+02 

• 1 *390+02 

• l *390+02 

• 1 *39q + q2 

• 1*390 + 02 

• i*39o*02 

• 1*39(3+02 

• | *39o+02 

• | *37q+02 

• 1*390+o2 



. . 


, - n. A n-7 

-*■ 1^3^5*02- 

- —46395+02 


* 1 A 1 9 r +r>? 

ft t 

• , Ai4c*n7 

* \ 



*T 











la 

•1*399+02 

• 1*399 + 02 

• 1 6399+02 

• 1 *399 + 02 

.1*399+02 

• 1 6399*02 

. 1*399+02 

. 1 *399+02 

. 1*399+02 

• 1*399 + 02 


— 

9 

,| 6809*02 

• 1*909*02 

- 1 *809 + 02 

• 1 *908 + 02 

, 1 *9q9+02 

•1*908+02 

.1*909+02 

•1*909*02 

.1*909+02 

.1*908+02 



e 

• 1 *907 *02 

• 1*907 + 02 

• 1*907 + 02 

• 1*907 + 02 

. 1 *9p7+02 

• 1 6907*02 

. 1 6907+02 

• 1*907 + 02 

. 1 *907 + 02 

• 6907+02 



7 

-1*91 1 +Q2 

- 1*91 l +02 

• 1*9 1 1*02 

- 1*91 ;+07 

•1*91 1+02 

•1*911 *02 

.1*911+02 

. 1 *911+02 

.16811+02 




* 

•|*91 6+Q2 

• 1*916 + 02 

• 1 *9 | 6 + 0? 

•1*916+02 

•1*916+02 

• ,*9 16*02 

. 1*914+02 

• l *9 l 6 + 02 

• I 69 1 6*02 



— 

5- 

. 1*820*02 

— - H*H20±02 "-1 *82-0+02 

. (-*820+02 

.1*920+02- 

*820+02 

,4*92(1+412- 

^l*‘i2HtQ2 





H 

. 1*923+02 

•1*923+02 

• 1*923 + 02 

• 1 *923 + 02 

. 1 *823 + 02 

• ,*923*02 

. 1*923+02 

. 1 *923 + 02 





3 
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